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. Abstract 

We review the cosmological implications of the flat directions of the Mini- 
mally Supersymmetric Standard Model (MSSM). We describe how field conden- 
sates are created along the flat directions because of inflationary fluctuations. 
The post-inflationary dynamical evolution of the field condensate can charge 
up the condensate with i? or L in a process known as Affleck-Dine baryogen- 
esis. Condensate fluctuations can give rise to both adiabatic and isocurvature 
density perturbations and could be observable in future cosmic microwave exper- 
iments. In many cases the condensate is however not the state of lowest energy 
but fragments, with many interesting cosmological consequences. Fragmenta- 
tion is triggered by inflation-induced perturbations and the condensate lumps 
will eventually form non-top ological solitons, known as Q-balls. Their proper- 



ties depend on how supersymmetry breaking is transmitted to the MSSM; if by 
gravity, then the Q-balls are semi-stable but long-hved and can be the source 
of all the baryons and LSP dark matter; if by gauge interactions, the Q-balls 

can be absolutely stable and form dark matter that can be searched for directly. 
We also discuss some cosmological applications of generic flat directions and Q- 
balls in the context of self-interacting dark matter, inflatonic solitons and extra 
dimensions. 
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1 Introduction 



The interplay between particle physics and cosmology plays an increasing role in un- 
derstanding the physics beyond the Standard Model (SM) [1] and the early Universe 
before the era of Big Bang Nucleosynthesis (BBN) [2, 3]. On both fronts we currently 
lack hard data. Above the electroweak scale E ~ 0(100) GeV, the particle content 
is largely unknown, while beyond the BBN scale T ~ 0{1) MeV, there is no direct 
information about the thermal history of the Universe. However, there are some ob- 
servational hints, as well as a number of theoretical considerations, which seem to be 
pointing towards a wealth of new physics both at small distances and in the very early 
Universe. Perhaps most importantly, new data is expected soon from accelerator ex- 
periments such as LHC and from cosmological measurements carried out by satellites 
such as MAP [4] and Planck [5]. 

In cosmology the recent observations of the cosmic microwave background (CMB) 
radiation, which has a temperature ~ 2.728 ± 0.004 K [6], have given rise to an era 
of precision cosmology. The Cosmic Background Explorer (COBE) satellite [7] first 
detected in a full-sky map a temperature perturbation of one part in 10^ at scales 
larger than 7 degrees [8] . The irregularities are present at a scale larger than the size 
of the horizon at the time when the microwave photons were generated and cannot be 
explained within the traditional hot Big Bang model [9]. The recent balloon experi- 
ments BOOMERANG [10] and MAXIMA [11], together with the ground-based DASI 
[12] experiment have established the existence of the first few acoustic peaks in the 
positions predicted by cosmic inflation [13, 14, 15, 9]. Inflation, a period of exponential 
expansion in the very early Universe, is a direct link to physics at energy scales that 
will not be accessible to Earth-bound experiments for any foreseeable future. Inflation 
could occur because a slowly rolling scalar fleld, the inflaton, dynamically gives rise 
to an epoch dominated by a false vacuum. During inflation quantum fluctuation are 
imprinted on space-time as energy perturbations which then are stretched outside the 
causal horizon. These primordial fluctuations eventually re-enter our horizon, whence 
their form can be extracted from the CMB (for a review, see [16, 9]). 
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Inflation can be considered as a model for the origin of matter since all matter arises 
from the vacuum energy stored in the infiaton field. However the present models do 
not give clear predictions as to what sort of matter there is to be found in the Universe. 
Prom observations we know that baryons constitute about 3% of the total mass [3] , 
whereas relic diffuse cosmic ray background virtually excludes any domains of anti- 
baryons in the visible Universe [17]. Almost 30% of the total energy density is in non- 
luminous, non-baryonic dark matter [18]. Its origin and nature is unknown, although 
various simulations of large scale structure formation suggest that there must be at 
least some cold dark matter (CDM), comprising of particles with negligible velocity, 
although there may also be a component of hot dark matter (HDM), comprising of 
particles with relativistic velocities [19]. The rest of the energy density is in the form 
of dark energy [20, 21]. 

The striking asymmetry in the baryonic matter has existed at least since the time 
of BBN and plays an important role in providing the right abundances for the light 
elements. The present Helium i^He), Deuterium [D) and Lithium abundances suggest 
a baryon density and an asymmetry relative to photon density of order 10~^° [3]. Such 
an asymmetry is larger by a factor of 10^ than what it should have been by merely 
assuming a initially baryon symmetric hot Big Bang [22] . Therefore baryon asymmetry 
must have been created dynamically in the early Universe. 

The origin of baryon asymmetry and dark matter bring cosmology and particle 
physics together. Within SM all the three Sakharov conditions for baryogenesis [23] 
are in principle met; there is baryon number violation, C and CP violation, and 
an out-of-equilibrium environment during a first-order electroweak phase transition. 
However, it has turned out that within SM the electroweak phase transition is not 
strong enough [24, 25, 26, 27], and therefore the existence of baryons requires new 
physics. Regarding HDM, light neutrinos are a possible candidate [19, 28], but there 
is no candidate for CDM in the SM. HDM alone cannot lead a successful structure 
formation because of HDM free streaming length [29, 19]. Therefore one must resort 
to physics beyond the SM also to find a candidate for CDM [18]. 

The tangible evidence for small but non- vanishing neutrino masses as indicated by 
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the neutrino oscillations observed by the Super-Kamiokande [30] and SNO collabora- 
tions [31] is definitely another indication for new physics beyond the SM. The main 
sources of neutrino mass could be cither Dirac or Majorana. A Dirac neutrino would 
require a large fine tuning in the Yukawa sector (one part in 10^^) while a Majorana 
mass would appear to require a scale much above the electroweak scale together with 
an extension of the SM gauge group SU{3)c x SU{2)l x U{1)y- In the Majorana case 
the lightness of the neutrino could be explained via the see-saw mechanism [32, 33]. 

A theoretical conundrum is that the mass scale of SM is ~ O(IOO) GeV, much lower 
that the scale of gravity Mp = (SttGjv)"^/^ = 2.436 x 10^^ GeV, and not protected 
from quantum corrections. The most popular remedy is of course supersymmetry (for 
a review, see [34, 35, 37]), despite the fact that so far supersymmetry has evaded 
all observations [38]. The minimal supersymmetric extension of the SM is called the 
MSSM. Supersymmetry must be broken at a scale ~ 0{1) TeV, presumably in some 
hidden sector from which breaking is transmitted to the MSSM, e.g., by gravitational 
[34, 35] or gauge interactions [39] . 

In the MSSM the number of degrees of freedom are increased by virtue of the 
supersymmetric counterparts of the SM bosons and fcrmions. One of them, known 
as the lightest supersymmetric particle (LSP), could be absolutely stable with a mass 
of the order of supersymmetry breaking scale. LSP would be a natural candidate for 
CDM (see e.g. [18]). In addition, because of the larger parameter space, electroweak 
baryogenesis in MSSM in principle has a much better chance to succeed. However, 
there are a number of important constraints, and lately Higgs searches at LEP have 
narrowed down the parameter space to the point where it has all but disappeared [40, 
41, 42]. 

Electroweak baryogenesis within MSSM thus appears to be heading towards deep 
trouble. Moreover, although MSSM can provide CDM, there is no connection between 
dark matter and electroweak baryogenesis. On the other hand, by virtue of supersym- 
metry, MSSM has the intriguing feature that there are directions in the field space 
which have virtually no potential. They are usually known as flat directions, which 
are made up of squarks and sleptons and therefore carry baryon number and/or lepton 
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number. The MSSM flat directions have been all classified [43] . 

Because it does not cost anything in energy, during infiation squarks and sleptons 
arc free to fluctuate along the flat directions and form scalar condensates. Because 
inflation smoothes out all gradients, only the homogeneous condensate mode survives. 
However, like any massless scalar field, the condensate is subject to infiaton-induced 
zero point fiuctuations which impart a small, and in infiation models a calculable, 
spectrum of perturbations on the condensate. After infiation the dynamical evolution 
of the condensate can charge the condensate up with a large baryon or lepton number, 
which can then released into the Universe when the condensate decays, as was first 
discussed by Affleck and Dine [44] . 

The potential along the MSSM flat direction is not completely flat because of su- 
persymmetry breaking. In addition to the usual soft supersymmetry breaking, the 
non-zero energy density of the early Universe also breaks supersymmetry, in par- 
ticular during inflation when the Hubble expansion dominates over any low energy 
supersymmetry breaking scale [45, 46]. Flatness can also be spoiled by higher-order 
non-renormalizable terms, and the details of the condensate dynamics depend on these. 

In most cases, the MSSM condensate along a flat direction is however not the state 
of lowest energy. The condensate typically has a negative pressure, which causes the 
inflation-induced perturbations to grow. Because of this the condensate fragments, 
usually when the Hubble scale equals the supersymmetry breaking scale, into lumps 
of condensate matter which eventually settle down to non-topological solitons dubbed 
as Q-balls by Coleman [47] . Q-balls carry a global charge, which in the case of MSSM 
is either B or L. 

The properties of Q-balls depend on supersymmetry breaking. If transmitted to 
MSSM by gravity, the Q-balls turn out to be only semistable but nevertheless long- 
hved compared with the time scales of the very early Universe [48] . When they decay, 
they may provide not only the baryonic matter but also dark matter LSPs [49]. If 
supersymmetry breaking is transmitted from the hidden sector to MSSM by gauge 
interactions, the resulting Q-balls would be stable and could exist at present as a 
form of dark matter [50] . In this case one can make direct searches for their existence 
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[51]. In both cases there is a prediction for the relation between the baryon and dark 
matter densities. Moreover, the condensate perturbations are inherited by the Q-balls, 
and can thus be a source of both isocurvature and adiabatic density perturbations 
[52, 53, 54]. 

This review is organized as follows. In Section 2, we recapitulate some basic cos- 
mology, and in particular baryogenesis. We briefly discuss various popular schemes 
of baryogenesis and describe the original Affleck-Dine baryogenesis. In Section 3, we 
present some background material for inflation, mainly concentrating on supersymmet- 
ric models. Quantum fluctuations and reheating are also discussed. In Section 4, we 
present the MSSM flat directions and discuss their properties. Various contributions to 
the flat direction potential in the early Universe are listed. Low energy supersymmetry 
breaking schemes, such as gravity and gauge mediation, are also discussed. In Section 
5, we discuss the dynamical properties of fiat directions and the running of the flat 
direction potential due to gauge and Yukawa interactions. Leptogenesis along LH^ 
flat direction, and the condensate evaporation in a thermal bath, is also described. 
We discuss fragmentation of the condensates for both gravity and gauge mediated 
supersymmetry breaking and present the relevant numerical studies. In Section 6, 
Q-ball properties are presented in detail. We describe various types of Q-balls, their 
interactions and their behavior at finite temperature. We discuss surface evaporation, 
diffusion, and dissociation of charge from Q-balls in a thermal bath. In Section 7, 
we focus on the cosmological consequences of Q-balls. We consider Q-ball baryogen- 
esis and non-thermal dark matter generation through charge evaporation for different 
types of Q-balls. We discuss Q-balls as self-interacting dark matter and present ex- 
perimental and astrophysical constraints on stable Q-balls. In Section 8, we briefly 
survey beyond-the-MSSM-condensates by considering inflatonic Q-balls and Affleck- 
Dine mechanism without MSSM flat directions. We also describe solitosynthesis, a 
process of accumulating large Q-balls in a charge asymmetric Universe. 
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2 Baryogenesis 



2.1 Baryon asymmetric Universe 

There are only insignificant amounts of anti-particles within the solar system. Cosmic 
ray showers contain ~ 10~^ anti-protons for each proton [55] , but the anti-protons are 
by-products of the interaction of the primary beam with the interstellar dust medium. 
This strongly suggest that galaxies and intergalactic medium is made up of matter 
rather than anti-matter, and if there were any anti-matter, the abundance has to be 
smaller than one part in 10^. The absence of annihilation radiation from the Virgo 
cluster shows that little anti-matter is to be found within a 20 Mpc sphere, and the 
relic diffuse cosmic ray background virtually excludes domains of anti-matter in the 
visible Universe [17]. 

The best present estimation for the baryon density comes from BBN [56] combined 
with the CMB experiments and it is given by [57] 

0.010 < Qb/i^ < 0.022 , (1) 

where flh = Pb/ Pc defines the fractional baryon density pb with respect to the critical 
energy density of the Universe: pc — 1.88 /i^ x 10~^^ g cm~^. The observational uncer- 
tainties in the present value of the Hubble constant; Hq — 100 h km • s~^ • Mpc~^ 
(/i/3000)Mpc~^ are encoded in h. Various considerations such as Hubble Space Tele- 
scope observations and type la supernova data suggest that h = 0.70 [58]. However, 
from the age of the globular cluster which comes out to be 11 Gyr, h seems to take 
lower value of about 0.5 [59]. The present convention is to take 0.5 < h < 0.8. In 
terms of the baryon and photon number densities we may write 

77 = ^^^^ = 2.68 X 10-U/i' , (2) 

where rib is the baryon number density and is for anti-baryons. The photon number 
density is given by = (2C(3)/7r^)T'^. Observations of the deuterium abundance in 
quasar absorption lines suggest [60] 

4(3) X 10"^° <rj< 7(10) x 10"^° . (3) 
11 



The conservative bounds are in parenthesis. 

Often in the hterature the baryon asymmetry is given in relation to the entropy 
density s = 1.8(?*n^, where g^, measures the effective number of relativistic species 
which itself a function of temperature. At the present time ~ 3.36, while during 
BBN g^ 10.11, rising up to 106.75 at T » 100 GeV. In the presence of supersym- 
metry at T » 100 GeV, the number of effective relativistic species are doubled to 
213.30. 

The baryon asymmetry B, defined as the difference of baryon and anti-baryon 
number densities relative to the entropy density, is bounded by 

5.7(4.3) X 10"^^ <B= ""'"""^ < 9.9(14) x 10"^^ , (4) 

s 

where the numbers in parenthesis are conservative bounds [60]. If at the beginning 
T) = 0, then the origin of this small number can not be understood in a CPT invariant 
Universe by a mere thermal decoupling of nucleons and anti-nucleons at T ~ 20 MeV. 
The resulting asymmetry would be too small by at least nine orders of magnitude, see 
[22]. 

2.2 Thermal history of the Universe 
2.2.1 Expanding Universe 

The hot Big Bang cosmology assumes that the Universe is spatially homogeneous and 
isotropic and can be described by the Priedmann-Robertson- Walker (FRW) metric 



ds^ = g^^dx^dx"" = dr - a\t) 



+ r\de' + sin' ed(t)' 



(5) 



1-Xr2 

where a{t) is the scale factor that determines the expansion or the contraction of the 
Universe; the constant K defines the spatial geometry. If X = 0, the Universe is 
flat and has Euclidean geometry, otherwise there is a spatial curvature corresponding 
either to a closed eUiptic [K = +1) or an open hyperbolic {K = —1) geometries. The 
value of K cannot however fix the global topology; for instance, an Euclidean topology 
can be flat and inflnite R^, or a surface of a 3-torus T^. However, topology has other 
observable consequences, e.g., for the pattern of CMB temperature fluctuations [61]. 
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There are two characteristic scales corresponding to the homogeneous and isotropic 
Universe: the curvature scale rcurv = o-{t)\K\~^/'^, and the Hubble scale 



H 



-1 



At). 

where dot denotes derivative w.r.t. t. The Hubble time is denoted by 

rf dt 



(6) 



m ^ . (7) 



^Hub / rj, — 1 

The behavior of the scale factor depends on the energy momentum tensor of the 
Universe. For a perfect fluid 

= -PQixi^ + {p + p)Ufj,u^ , (8) 

where p is the energy density and p is the pressure of a fluid and the four velocity 
Ufj, = dx^/ds. For the FRW metric and for the perfect fluid the equations of motion 
gives the Friedmann-Lemaitre equation 

3M| a(t)2' ^""^ 
also known as the Hubble equation. The acceleration equation is given by 

and the conservation of the energy momentum tensor Tf^^ — gives 

^ = -3pa^ (11) 
da 

Note that pa^ is constantly decreasing in an expanding Universe for a positive pressure. 

The early Universe is believed to have been radiation dominated with p — p/3 and 
a{t) oc t^^'^, followed by a matter dominated era with p — and a{t) oc i^/^. The early 
Universe might also have had an era of acceleration, known as the inflationary phase, 
which could have happened only if 

d>0 ^ p + 3p<0. (12) 



A geometric way of defining inflation is [9] 

d{H-ya{t)) 
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which states that the Hubble length as measured in comoving coordinates decreases 
during inflation. We will use this particular deflnition of inflation very often while 
discussing the number of e-foldings and density perturbations. 
The Hubble expansion rate is related to the temperature by 



where is the total number of relativistic degrees of freedom and it is given by 

9.{T) = Y.9^(^y + lJ:9^(^y. (15) 



. . sTj 8 ^ VTx 

i=b 1=1 

Here Ti denotes the effective temperature of species i, which decouples at a temperature 

During the radiation era when H = (l/2t), one flnds 

±«2.42,,'/^(i^)\ (16) 

2.2.2 Entropy 

An ideal gas of particles respects the Fermi-Dirac or Bose-Einstein distributions 

/,(p,/i,T) = [exp((£;, - f,,)/T) If 1]-^ , (17) 

where Ef — |p|^ + m^, jii represents the chemical potential of the species i, — /+ cor- 
responds to Bose/Fermi statistics. The value of ji is equal and opposite for particles 
and anti-particles. Therefore in the early Universe a flnite net chemical potential is 
proportional to the particle anti-particle asymmetry. The bound on charge asymmetry 
relative to the photon number density is severe, less than one part in 10^^ at temper- 
atures close to BBN [62], while baryon asymmetry is comparatively larger, but still 
small enough for //g, ~ to be an excellent approximation. Neutrinos may however 
carry a net B — L charge which need not be vanishingly small at early times, although 
a large enough neutrino chemical potential can affect nucleosynthesis, for example, 
see [28]. 

The number density n, energy density p, and pressure p can be expressed in terms 
of temperature, and Qi is the number of internal degrees of freedom [63] 

"^^^^ ^ (If/ /^(P'^'^)^'p = , 
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P.(T) = (^/ E,^p)f.(p,^,,T)d'p = ^T'Jf (=F) , 



lfM^r "Y.Zir\ , (19) 



where 

For a relativistic case with <S 1, 

= 2C(3), = = for bosons, 

= /f(+) = 4'^^(+) = g^, forfermions, (20) 

where ( denotes the Riemann Zeta function and ^(3) = 1.202. Thus the energy density 
of radiation reads 

P^-^9*T'. (21) 
For non- relativistic particles with Xi ^ 1, one obtains for both bosons and fermions 

nnr{T) = ^ = aA^) e-^/\ p„r = 0. (22) 

If the chemical potential is non-zero, the exponential Eq. (22) also includes a factor 

The entropy density is defined as 

= (23) 

where d{sa^) = is a thermodynamically conserved quantity. The decoupling tem- 
perature can be expressed as [64] 

where S is the total entropy and Sa the entropy in the degrees of freedom that have 
decoupled at To- 
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2.2.3 Nucleosynthesis 

According to BBN (for reviews see [3, 57]) the light elements ^H, ^He, ^He, and Li 
have been synthesized during the first few hundred seconds. The abundances depend 
on the baryon-to-photon ratio 

-J-^. (25) 

All the relevant physical processes take place essentially in the range from a few MeV 
~ 0.1 sec down to 60 — 70 KeV ~ 10^ sec. During this period only photons, pairs, 
and the three neutrino flavors contribute signiflcantly to the energy density. Any 

additional energy density may be parameterized in terms of the effective number of 
light neutrino species N^, so that 

= 10.75+^(iV,-3). (26) 

Nucleosynthesis starts off with a freezing out of the weak interaction between neutron 
and proton at To Ri 0.8 MeV. Free neutrons keep decaying until deuterium begins to 
form through n+p d+j. Deuterium synthesis is over by ^ 0.086 MeV (assuming 
?7 = 5 X 10~^°). At To, neutron abundance has been depleted to X„(t£)) = n/{n+p) f» 
0.122. All the surviving neutrons are now captured through n-\- D ^ i^H,^ He), and 
subsequently by virtue of the process i^H,^ He) + n — >^ He, which has a binding 
energy of 28.3 MeV. The total mass fraction of primordial hehum, which is denoted 
by Fp(^ife), is given by 

YY>{^He) ^ 2X^(Td) = 0.245 . (27) 

Adopting the experimentally allowed range of 0.22 < Fp < 0.26, one can constraint 
that number of light neutrino species by [56] 

AT, < 4 . (28) 

The four LEP experiments combined give the best fit as [65] 

= 2.994 ± 0.12 . (29) 

Nucleosynthesis also constrains many non-conventional ideas, for instance alternative 
theories of gravity such as scalar-tensor theories [66] . 
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Besides ^He, D and ^He are produced at the level of 10 ^, and ^Li at the level 
of 10~^°. The theoretical prediction has some slight problems in fitting the observed 
'^He and ^Li abundances. Both seem to indicate 1.7 x 10"^'^ < f] < 4.7 x 10~^°, 
corresponding to 0.006 < Ofe/i^ <- o.017 with a central value n^h^ = 0.009 [57]. The 
abundance ratio D/H is comparable with '^He and "^Li abundances at the 2a level in 
the range 4.7 x 10"^° < rj < 6.2 x 10"^°, which corresponds to 0.017 < Q^/i^ < 0.023. 
The likelihood analysis which includes all the three elements (D,^ife,and '^Li) yields 
[57] 

4.7 X 10"^° < ?7 < 6.2 X 10"^° , 0.017 < Qbh'^ < 0.023 . (30) 

Despite the uncertainties there appears to be a general concordance between theo- 
retical BBN predictions and observations, which is now being bolstered by the CMB 
data from several different experiments. The results from the ground based DASI ex- 
periment indicates Q^/i^ = 0.0221^:^^3 [12], while the results from the BOOMERANG 
balloon-borne experiment imply Q^/i^ = 0.02llo:oo3 [10]- MAXIMA, another balloon 
experiment, quotes a somewhat larger value flhh^ — 0.0325 ± 0.006 [67]. 

2.3 Requirements for baryogenesis 

As pointed out by Sakharov [23], baryogenesis requires three ingredients: (1) baryon 
number non- conservation, (2) C and CP violation, and (3) out-of-equilibrium condi- 
tion. All these three conditions are believed to be met in the very early Universe^ . 

2.3.1 Non- conservation of baryonic charge 

In the SM, baryon number B is violated by non-perturbative instanton processes 
[70, 71]. At the quantum level both baryon number current and the lepton number 

current J£ are not conserved because of chiral anomalies [72]. However the anoma- 
^ There have been attempts [68, 69] for baryogenesis via a repulsive interaction between baryons 
and anti-baryons which would lead to their spatial separation before thermal decoupling of nucleons 
and anti-nucleons. However at such early times the causal horizon contained only a very small fraction 
of the solar mass so that the asymmetry could not be smooth at distances greater than the galactic 
size. 
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lous divergences of and Jj^ come with an equal amplitude and an opposite sign. 
Therefore B — L remains conserved, while B + L may change via processes which 
interpolate between the multiple non-Abelian vacua oi SU{2). The probability for the 
B + L violating transition is however exponentially suppressed [70, 71]. As was first 
pointed out by Manton [73], at high temperatures the situation is different, so that 
when T » Mw, baryon violating transitions are in fact copious (see Sect. 2.4.2). 

In addition to baryogenesis, B violation also leads to proton decay in GUTs. For 
instance, the dimension 6 operator (QQQL)/ A generates observable proton decay 
unless A > 10^^ GeV. In the MSSM the bound is A > lO^^ GeV because the decay 
can take place via a dimension 5 operator. In the MSSM superpotential there are also 
terms which can lead to AL = 1 and AB — 1. Similarly there are other processes 
such as neutron-anti-neutron oscillations in SM and in supersymmetric theories which 
lead to AB = 2 and AB = 1 transitions [74]. These operators are constrained by the 
measurements of the proton lifetime, which yield the bound Tp > 10^^ years [65]. 

2.3.2 C and CP violation 

Weak interactions ensures maximum C violation while neutral Kaon is an example 
of CP violation in the quark sector which has a relative strength ~ 10~^ [65]. CP 
violation could also expected to be found in the neutrino sector. Beyond the SM there 
are many sources for CP violation. An example is the axion proposed for solving the 
strong CP problem in QCD [75]. Quantum fluctuations of light scalars in the early 
Universe, in particular during inflation, can create different domains of various C and 
CP phases. C and CP can also be spontaneously broken during a phase transition, 
so that domains of broken phases form with different CP-charges [76]. 

2.3.3 Depeirture from thermal equilibrium 

Departure from a thermal equilibrium cannot be achieved by mere particle physics 
considerations but is coupled to the dynamical evolution of the Universe. If S- violation 
processes are in thermal equilibrium, the inverse processes will wash out the pre- 
existing asymmetry (Anb)o [77]. This is a consequence of (S'-matrix unitarity and 
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CPT-theorem [78]. However there are several ways of obtaining an out-of-equilibrium 
process in the early Universe. 

• Out-of-equilibrium decay or scattering: 

The Universe in a thermal equilibrium can not produce any asymmetry, rather 
it tries to equilibrate any pre-existing asymmetry. If the scattering rate T < H, 
the process can take place out-of-equilibrium. Such a situation is appropriate 
for e.g. GUT baryogenesis [78, 79]. 

• Phase transitions: 

Phase transitions are ubiquitous in the early Universe. The transition could be 
of first, or of second (or of still higher) order. First order transitions proceed by 
barrier penetration and subsequent bubble nucleation resulting in a temporary 
departure from equilibrium. Second order phase transitions have no barrier be- 
tween the symmetric and the broken phase. They are continuous and equilibrium 
is maintained throughout the transition. 

Prime examples of first order phase transitions in the early Universe are the 
QCD and electroweak phase transitions. The nature and details of QCD phase 
transition is still very much an open debate [80, 25, 81], and although a mech- 
anism for baryogenesis during QCD phase transition has been proposed [82], 
much more effort has been devoted to the electroweak phase transition [24, 26] 
(see Sect. 2.4.2). 

• Non-adiabatic motion of a scalar field: 

Any complex scalar field carries C and CP, but the symmetries can be broken 
by terms in the scalar potential. This can lead to a non-trivial trajectory of 
a complex scalar field in the phase space. If a coherent scalar field is trapped 
in a local minimum of the potential and if the shape of the potential changes 
to become a maximum, then the field may not have enough time to readjust 
with the potential and may experience completely non-adiabatic motion. This 
is similar to a second order phase transition but it is the non-adiabatic classical 
motion which prevails over the quantum fluctuations, and therefore, departure 
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from equilibrium can be achieved. If the field condensate carries a global charge 
such as the baryon number, the motion can charge up the condensate. This is 
the basis for the Affleck-Dine baryogenesis [44] (see Sect. 2.5). 

2.3.4 Sphalerons 

In the SM B + L is very weakly violated in the vacuum [70] . At finite temperatures 
violation is large [73, 83, 84, 85, 86] by virtue of the sphaleron configurations, which 
mediate transitions between degenerate gauge vacua with different Chern-Simons num- 
bers related to the net change oi B + L. Thermal scattering produces sphalerons which 
in effect decay in B + L non-conserving ways below 10^^ GeV [87], and thus can ex- 
ponentially wash away B + L asymmetry. Sphalerons and associated electroweak 
baryogenesis has been reviewed in [88, 24, 26, 89, 90, 91]. Let us here just give a brief 
summary of the main ingredients. 

• Chiral anomalies 

An anomaly means that a classical current conservation no longer holds at the 
quantum level; an example is the chiral anomaly [72] . In the SM there is classical 
conservation of the baryon and lepton number currents Jg and J|^, but because 
of chiral anomaly the currents are not conserved. Instead [70], 

MS Stt ^ ' 87r ^9 9 36/ ^' 

« = -^N,WrW,,^ + ^N,{l-]^F'^'^F^,, (31) 

where Ng is the number of generations, a2 and cci {Winv and F^^,) are respectively 
the SU{2) and U{1) gauge couplings (field strengths), and the various numbers 
inside the brackets correspond to the squares of the hypercharges multiplied by 
the number of states. Note that while at the quantum level B + L is violated, 
B — L is still conserved. 

• Gauge theory vacua 

The vacuum structure of the gauge theories is very rich [71, 92] . In case of SU (2) , 
the vacua are classified by their homotopy class {Q„(r)}, characterized by the 
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winding number n which labels the so called ^-vacua [71, 92]. A gauge invariant 
quantity is the difference in the winding number (Chern-Simons number) 

Ncs = n+-n^^^l d^xW^Wa^. ■ (32) 

In the electroweak sector the field density WW is related to the divergence of 
B + L current. Therefore a change in B + L reflects a change in the vacuum 
conflguration and is determined by the difference in the winding number 

A{B + L)^J d^'xd^rs+L = -^Ng J d'^xWrWa^u = -2NgNcs ■ (33) 

For three generations of SM Icptons and quarks the minimal violation is A{B + 
L) = 6. Note that the proton decay p e~^7i^ requires A{B + L) = 2, so that 
despite violation, proton decay is completely forbidden in the SM. 

The probability amplitude for tunneling from an n vacuum at i ^ — oo to an 

n + Ncs vacuum at t — >• +oo can be estimated by the WKB method [70] 

mc.W~exp(z^^)~10— <'^. (34) 

Therefore, as advertised, the baryon number violation rate is totally negligible 
in the SM at zero temperature, but as argued in a seminal paper by Kuzmin, 
Rubakov and Shaposhnikov [93], at flnite temperatures the situation is com- 
pletely different. 

• Thermal tunneling 

The sphaleron is a fleld conflguration sitting at the top of the potential barrier 
between two vacua with different Chcrn-Simons numbers and can be reached 
simply because of thermal fluctuations [93]. Neglecting U{1)y, the zero temper- 
ature sphaleron solution was flrst found by Manton and Klinkhamer [73, 83]. 
At flnite temperature the energy obeys an approximate scaling law [85, 86] 

E,,,{T) = E,,,{omT))/{my- 

E,,,(T)^^^^^B(X/g,), (35) 

where m^(T) = {l/2)g2{^{T)) is the mass of the W-boson and the function 
B has a weak dependence on X/ g2, where A is the quartic self coupling of the 
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Higgs. Below the critical temperature of the electroweak phase transition, the 
sphaleron rate is exponentially suppressed [94]: 

r ~ 2.8 X lO^^T^ ' e-^-''/^ . (36) 

where n is the functional determinant which can take the values 10~^ < < 10~^ 
[95] . Above the critical temperature the rate is however unsuppressed. Requiring 
that the Chern-Simons number changes at most by ANqs ~ 1, one can estimate 
from Eq. (32) that ANcs ~ g^^hWl ~ 1 ^ ~ Therefore a typical 

energy of the sphaleron configuration is given by 

Esph ~ l%Himr ~ . (37) 

i/2 ''sph 

At temperatures greater than the critical temperature there is no Boltzmann 
suppression, so that the thermal energy oc T > Egph- This determines the size 
of the sphaleron as 

hph > ■ (38) 
This is exactly the infrared cut-off generated by the magnetic mass of order 
~ 5f|T. Therefore, based on this coherence length scale one can estimate the 
baryon number violation per volume ~ Z^^^, and per unit time ~ I sph- On 
dimensional grounds the transition probability would then be given by 



''sph'' 



K{a2TY . (39) 



where k is a constant which incorporates various uncertainties. However, the 
process is inherently non-perturbative, and it has been argued that damping of 
the magnetic field in a plasma suppresses the sphaleron rate by an extra power 
of a2 [96], with the consequence that Tgph ~ Q;|T^. Lattice simulations with hard 
thermal loops also give Tsph ~ 0{lQ)alT^ [97] 

2.3.5 Washing out S + L 

In the early Universe the transitions ANcs = +1 and ANcs = — 1 are equally proba- 
ble. The ratio of rates for the two transitions is given by 

£fp^^g-AF/r^ (40) 

sph — 
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where AF is the free energy difference between the two vacua. Because of a finite 
B + L density, there is a net chemical potential Hb+l- Therefore 



AF ~ f,l_,^T' + 0{T') ^ ^ + 0{T') . (41) 
One then obtains [87] 

^ = r.,, ^ - r.,, _ r. N.^UB^L . (42) 

It then follows that an exponential depletion of ub+l due to sphaleron transitions 
remains active as long as 

^>H ^ T<at^^lO''GeY. (43) 
^ 9* 

This result is important because it suggests that below T = 10^^ GeV, the sphaleron 

transitions can wash out any B + L asymmetry being produced earlier in a time scale 

r ~ {T^/NgTsph}- This seems to wreck GUT baryogenesis based on B — L conserving 

groups such as the minimal SU{5). 

2.4 Alternatives for baryogenesis 

There are several scenarios for baryogenesis (for reviews on baryogenesis, see [78, 90, 
91]), the main contenders being GUT baryogenesis, electroweak baryogenesis, lepto- 
genesis, and baryogenesis through the decay of a field condensate, or Affleck-Dine 
baryogenesis. Here we give a brief description of these various alternatives. 

2.4.1 GUT-beiryogenesis 

This was the first concrete attempt of model building on baryogenesis which incorpo- 
rates out-of-equilibrium decays of heavy GUT gauge bosons X, y — > qq, and X,Y ^ ql 
(see e.g., [78, 79, 88, 98]). The decay rate of the gauge boson goes as Fx ~ axMx, 
where Mx is the mass of the gauge boson and Ox^ is the GUT gauge couphng. As- 
suming that the Universe was in thermal equilibrium at the GUT scale, the decay 
temperature is given by 

Tn^9:'/'a'^\MxMpy/\ (44) 
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which is smaller than the gauge boson mass. Thus, at T T^, one expects nx ^ 

Tlx ^ 

and hence the net baryon density is proportional to the photon number 
density Ub = ABrij. However below Td the gauge boson abundances decrease and 
eventually they go out-of-equilibrium. The net entropy generated due to their decay 
heats up the Universe with a temperature which we denote here by Trh- Let us naively 
assume that the energy density of the Universe at Td is dominated by the X bosons 
with px ~ MxUx, and their decay products lead to radiation with an energy density 
p = {n'^/ 30) g^T^f^, where ~ C(IOO) for T > Mqut- Equating the expressions for 
the two energy densities one obtains 

Therefore the net baryon number comes out to be 

^ ^ ^ ASrzx ^ 3 Tw,^^ (46) 
s g^n^ AMx 

Trh is determined from the relation r| ^ H'^{Td) ~ (n^ / 90) g^T^/M^. Thus, 

Uncertainties in C and CP violation are now hidden in AB, but can be tuned to yield 
total B ~ 10^^° in many models. 

Above we have tacitly assumed that the Universe is in thermal equilibrium when 
T > Mx- This might not be true, since for 2 <-> 2 processes the scattering rate is 
given by r ~ a^T, which becomes smaller than H at sufficiently high temperatures. 
Elastic 2 — > 2 processes maintain thermal contact typically only up to a maximum 
temperature ~ 10^'' GeV, while chemical equilibrium is lost already at T ~ 10^^ GeV 
[99, 100]. It has been argued that QCD gas, which becomes asymptotically free at 
high temperatures, never reaches a chemical equihbrium above ~ 10^^ GeV [101]. 
In supergravity the maximum temperature of the thermal bath should not exceed 
10^0 GeV [102] (see Sect. 3.6.3). 
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2.4.2 Electroweak baryogenesis 

A popular baryogenesis candidate is based on the electroweak phase transition, during 
which one can in principle meet all the Sakharov conditions. There is the sphaleron- 
induced baryon number violation above the critical temperature, various sources of CP 
violation, and an out-of-equilibrium environment if the phase transition is of the first 
order. In that case bubbles of broken SU{2) x U{1)y are nucleated into a symmetric 
background with a Higgs field profile that changes through the bubble wall [93, 103, 24]. 

There are two possible mechanisms which work in a different regime; local and non- 
local baryogenesis. In the local case both CP violation and baryon number violation 
takes place near the bubble wall. This requires the velocity of the bubble wall to be 
greater than the speed of the sound in the plasma [104, 105], and the electroweak 
phase transition to be strongly first order with thin bubble walls. 

The second alternative, where the bubble wall velocity speed is small compared to 
the sound speed in the plasma, appears to be more realistic. In this mechanism the 
fermions, mainly the top quark and the tau-lepton, undergo CP violating interactions 
with the bubble wall, which results in a difference in the reflection and the transmission 
probabilities for the left and right chiral fermions. The net outcome is an overall chiral 
flux into the unbroken phase from the broken phase. The flux is then converted into 
baryons via sphaleron transitions inside the unbroken phase. The interactions are 
taking place in a thermal equilibrium except for the sphaleron transitions, the rate of 
which is slower than the rate at which the biibble sweeps the space. 

One great difficulty with the electroweak baryogenesis is the smallness of CP vio- 
lation in the SM. It has been pointed out that an additional Higgs doublet [106, 107, 
105, 108] would provide an extra source for CP violation in the Higgs sector. However, 
the situation is much improved in the MSSM where there are two Higgs doublets 
and Hd, and two important sources of CP violation [109]. The Higgses couple to the 
charginos and neutralinos at one loop level leading to a CP violating contribution. 
There is also a new source of CP violation in the mass matrix of the top squarks 
which can give rise to considerable CP violation [110]. 

Bubble nucleation depends on the thermal tunneling rate, and the expansion rate 
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of the Universe. The tunnehng rate has to overcome the expansion rate in order to 
have a successful phase transition via bubble nucleation at a given critical temperature 
Tc > Tt > Tq. The actual value of the baryon asymmetry produced at the electroweak 
baryogenesis is still an open debate [111, 40, 112, 113, 42], but in general it is hard 
to generate a large baryon asymmetry. For Tc ~ 100 GeV, N — 3, a2 — 0.033, and 
B{\/g2) ~ 1.87, one obtains the condition [114, 87, 115, 26] 



which implies [115] 



EsphiTc) ^^^^^ 



T 

-'- c 



+ 91og(10)+log(K). (48) 



Esph{Tc) ^ f„„ 1A-1 



> 45 for K = 10"' , (49) 

c 

> 37 for K = 10"^ . (50) 



The standard bound is often taken to be that of Eq. (49). In terms of the Higgs field 
value at Tc, one then obtains from Eq. (35) 

for the above values of q;2, B. Then the bounds in Eqs. (49,50) translate to 

^>l-3 (1), (52) 

where the number in parenthesis is for Eq. (50). Eq. (52) respectively, implies that 
the phase transition should be strongly first order in order that sphalerons do not 
wash away all the produced baryon asymmetry. This result is the main constraint on 
electroweak baryogenesis. 

In order to save the SM electroweak baryogenesis some attempts such as matter 
induced effects have also been invoked; exciting the SM gauge degrees of freedom in a 
time varying Higgs background [116, 117], or via dynamical scalar field which couples 
to the SM fields [118]. 

Lattice studies suggest that in the SM the phase transition is strongly first order 
only below Higgs mass itlh ~ 72 GeV [27, 119, 120, 121]. Above this scale the 
transition is just a cross-over. Such a Higgs mass is clearly excluded by the LEP 
measurements [65], thus excluding electroweak baryogenesis within the SM. 
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2.4.3 Electroweak baryogenesis in MSSM 

In the MSSM the ratio ^{Tc)/Tc can increase considerably. The MSSM Higgs sector 
at finite temperature has been considered in [122, 123, 124, 125], for lattice studies see 
[126, 127, 128]. In the MSSM the right handed stop Ir couples to the Higgs with a large 
Yukawa coupling. This leads to a strong first order phase transition [123, 124, 125]. 
The LEP precision tests then indicate that the lightest left handed stop should be 
heavy heay^ with mq > 500 GeV. This implies the for lightest right handed stop mass 

mf ^ml + 0.15M| cos(2/?) + ( 1 - ^ ) , (53) 

where At — At — /i/ tan(/3) is the stop mixing parameter, and /i is the soft-SUSY 
breaking mass parameter for the right-handed stop. The coefficient P of the cubic 
term PTH^ in the effective potential reads 

and can be at least one order of magnitude larger than (isM- In principle this modifi- 
cation can give rise to a strong enough first order phase transition. 

The sphaleron bound implies Higgs and stop masses in the range [41, 42] 

110 GeV <mH< 115 GeV , and 105 GeV < mj^ < 165 GeV . (55) 

The present LEP constraint on the Higgs mass is mn > 115 GeV [129]. Hence, even 
an MSSM-based electroweak baryogenesis may be at the verge of being ruled out. The 
definitive test of the MSSM based electroweak baryogenesis will obviously come from 
the Higgs and the stop searches at the LHC and the Tevatron [40, 130, 41, 111, 42]. 

2.4.4 Leptogenesis 

Even if B+L is completely erased by the sphaleron transitions, a net baryon asymmetry 
in the Universe can still be generated from a non- vanishing B — L [131], even if there 
were no baryon number violating interactions. Lepton number violation alone can 
produce baryon asymmetry B ~ —L [132], a process which is known as leptogenesis 
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(for a recent review [133], and references therein). For lepton number violation one 
however has to go beyond the SM. 

A popular example is 5*0(10) GUT model, which can either be broken into SU{5) 
and then subsequently to the SM, or into the SM gauge group directly. The most 
attractive aspect of SO{10) is that it is left-right symmetric (for details, see [134, 135]), 
and has a natural foundation for the see-saw mechanism [32, 33] as it incorporates a 
singlet right-handed neutrino Nr with a mass Mr. A lepton number violation appears 
when the Majorana right handed neutrino decays into the SM lepton doublet and Higgs 
doublet, and their CP conjugate state through 

Nr^<!> + 1, + (56) 

There also exist AL = 0, and AL = 2 processes mediated by the right handed neutrino 
through 

' ' ^^^^ 

which are dimension 5 operators [136, 137]. (There are other processes involving t- 
quarks which may also be important [138, 139]). CP asymmetry is generated through 
the interference between tree level and one-loop diagrams. 

The total baryon asymmetry and total lepton asymmetry can be found in terms of 
the chemical potentials as [140] 

B = I^(2/Xgi + Huui + l^dRi) , L = ^{2nii + HeRi) , (58) 

i i 

where i denotes three leptonic generations. The Yukawa interactions establish an 
equilibrium between the different generations {nu — and /igi — /ig, etc.), and one 

obtains expressions for B and L in terms of the number of colors iV = 3, and the 
number of charged Higgs fields Nh 

AN UN-' + gNNn 

B — — ^ — — 7rT7 — l^i , (59) 

3 ^ 6N + 3Nh 



together with a relationship between B and B — L [140] 

/ 8N + ANh 
\22N+13Nh 



B=(Z^'^" ](B-L). (60) 
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A similar expression has also been found in [141, 88, 133], although there seems to be 
small of order one differences. The baryon asymmetry based on the decays of right 
handed neutrinos in a thermal bath has been computed in [138, 142, 143, 144, 145]. 
In a recent analysis [145], it was pointed out that the baryogenesis scale is tightly 
constrained together with with the heavy right handed neutrino mass Tg ~ Mi^r — 
O(10^°) GeV, with an upper bound on the hght neutrino masses J2i f^i < eV. The 
current bound on the right handed neutrino mass is around Mr ~ C(10^^) GeV for 
light neutrino masses miu ~ ~ m-^^ ~ 0(0.1) eV. 

High scale leptogenesis is ruled out in a supersymmetric theory because of the grav- 
itino problem (see Sect. 3.7.1). However, if the masses of the right handed neutrinos 
are such that the mass splitting is comparable to their decay widths, it is possible to 
obtain an enhancement in the CP phase of order one [144], and possibly a low scale 
thermal leptogenesis [146]. Otherwise, one could resort to non-thermal leptogenesis 
[147, 136, 137, 148], or, to the scattering process discussed in [149], or to sneutrino 
driven leptogenesis [150, 137]. 

2.4.5 Baryogenesis through field condensate decay 

Scalar condensates may have formed in the course of the evolution of the early Uni- 
verse. In particular, during inflation all scalar flelds are subject to fluctuations driven 
by the non-zero inflaton energy density so that flelds with very shallow potentials 
may easily take non-zero values. An example is the MSSM, where for the squark 
and slepton fields there are several directions in the field space where the potential 
vanishes completely [46, 43]. These directions are called (perturbatively) flat. Field 
fluctuations along such flat directions will soon be smoothed out by inflation, which ef- 
fectively stretches out any gradients, and only the zero mode, or the scalar condensate, 
remains. This mechanism is quite general and applicable to any order parameter with 
flat enough potential. Baryogenesis can then be achieved by the decay of a condensate 
that carries baryonic charge, as was flrst pointed out by Affleck and Dine [44]. As we 
will discuss, the flat direction condensate can get dynamically charged with a large B 
and/or L by virtue of CP- violating self-couphngs. 
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Baryogenesis from MSSM flat directions has the virtue that it only requires two 
already quite popular paradigms: supcrsymmctry and inflation. In the old version 
[44] baryons were produced by a direct decay of the condensate, to be discussed in 
Sect. 2.5.2. It was however pointed out first by Kusenko and Shaposhnikov [151] in the 
case of gauge mediated supersymmetry breaking, and then by Enqvist and McDonald 
in [48] in the case of gravity mediated supersymmetry breaking, that the MSSM flat 
direction condensate most often is not stable but fragments and eventually forms non- 
topological solitons called Q-balls [47]. These issues will be dealt in Sects. 6 and 7. 

2.5 Old Affleck-Dine baryogenesis 
2.5.1 Classical motion of the order parameter 

In the original Affleck-Dine baryogenesis [44] it was assumed that the order parameter 
along the fiat direction is displaced from the origin because of inflationary fluctuations. 
Because of inflation, only the long wave-length model of the order parameter will 
survive so that a spatially constant condensate field is formed along the fiat direction. 
This we shall sometimes call the Affleck-Dine (AD) field. In an expanding Universe 
the coherent AD field (j) obeys the usual equation of motion, 

dV 

(j) + 3H(f>+ — = , (61) 

where H is the Hubble parameter. 

To follow the time evolution of the AD field, let us consider a toy model with the 
potential 

V{cj>) = m'\cf>\' + A(0^ + r') + ^ + ■ ■ ■ ■ (62) 

Although this potential is unrealistic in that it does not take correctly into account 
of supersymmetry breaking induced by the non-zero cosmological constant of the in- 
fiationary era, it nevertheless captures the main features of the initial cosmological 
evolution of the AD field. 

The theory Eq. (62) has a partially conserved current — with 

d^^f = dM*d^'<P - id^'<P*<P) = iX{<P*^ - (t>^) . (63) 
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The current is conserved for small 0. The role of the higher order term |0|^ is just 
to stabilize the potential. In the toy model Eq. (62), we identify the baryon number 
density Ub with jo- The model also has a CP invariance under which 0* but 
which is violated by the initial conditions, which are taken to be 



4> = i4>o, = 0, 



(64) 



where 0o is real. Writing 
e.g., [152]) 



'^R + i4>i one finds the coupled equations of motion (see 



+ 2,H(t>R + 



M2 
M2 



^R = 4A(30,0|-0f) 



(65) 



In a matter dominated Universe H = 2/(3t), so that for large times t ^ m ^ the 
motion is damped and Eq. (65) has then oscillatory solutions of the form 



(f)k — — sinimt + (5fe) , k — I,R , 
mt 



(66) 



where the amplitudes and the phases Sk depend on the parameters m, A, M, and 
the initial conditions Eq. (64). For large times the baryon number is then found to be 



Ub = 2(0/0ij - (l)R(j)i) = — ^ sin((5j - Sr) . 



(67) 



If 00 -C mM, as was tacitly assumed by Affleck and Dine [44], one may disregard 
the higher-order terms. In that case one obtains [45] Ai = 0o and Ar = aR\(f)l/m'^, 
where qr — 0.85 is determined numerically. Likewise, numerically one finds that 
5i-Sr^ 1.54. Thus, 



riB = 



1.7A04 



and the generated baryon number per particle is 

mnB 1.7A0O 



R 



P4> 



(68) 



(69) 



Eq. (69) is true for matter dominated Universe; for radiation dominated Universe one 
obtains a similar result, but the numerical prefactor 1.7 should be replaced by —1.3. 



31 



If 4)1 > mM, one finds Ai = ai{mMy/'^ and Ar = aR{M^ /mY''^ witfi o/ = 
0.94, qr = -2.86, 6i = Oil, and Sr = -0.41. It then follows from Eq. (67), that 

Thus, the baryon generation mechanism is remarkably robust. The initial conditions 
do not matter, nor the actual expansion rate of the Universe. The baryon number 
generated per 0-particle is always large and, with A ~ m^/(0)^, typically ub ^ 1- 
Although these conclusions were derived in a toy model, similar results hold true also 
for the MSSM flat directions. 

Thus, to summarize, along a flat direction where squarks and sleptons have non- 
zero expectation values, evolution of the AD field condensate, starting from a CP 
violating initial value, will dynamically generate large baryon number density and 
charge the condensate with B and/or L. 

2.5.2 Condensate decay 

To provide the Universe with the observed baryon to entropy ratio, ub/s 10"^'^, 
the AD condensate must eventually transform itself into ordinary quarks. Originally 
[44] , it was thought that this could happen via the decay of the AD field components 
(squarks and sleptons) to ordinary quarks and leptons. The AD condensate can be 
thought of as a coherent state of 0-particles where cf) — 0oe*"** and |0o| 3> m. When 
supersymmctry breaking is switched on, the AD field starts to oscillate about the old 
vacuum (0) ^ m. Writing 4> = {4>) + 4>', one observes that all the fields to which the 
excitations 0' couple are heavy with masses 0{{(f))). The field (f)' itself has a mass of 
the order of supersymmctry breaking, 0{m). Therefore, 0' can decay to fight fields 
only through loop diagrams involving heavy fields, with an effective coupling of the 
typs {g"^ / {4'))4>''^dilj\ where is a light fermion and g some coupling constant. The 
decay rate is thus, [44] 
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r~9*^. (71) 

Because of the oscillations of the AD field, the Universe will eventually become 
dominated by the energy density in the oscillations, ~ m^0^, so that H ~ plp"^ /Mp. 
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The AD field will decay when F ~ if, or ~ (m^Mp)^/^. This implies a reheating 



Depending on A, and the size of the initial fluctuation 0o of the AD condensate, 
ub/ s can be either small or large. Therefore determining the initial value is of utmost 
importance [153]. This requires us to consider theories of inflation in more detail, which 
will be done in the next Section. Following that, we shall discuss the disappearance of 
the AD condensate by fragmentation into (quasi)stable lumps of condensate matter, 
whose state of lowest energy is a spherical non-topological sohton, a Q-ball [47]. 



temperature T^h — s^^^ ~ p^"^ while the baryon number density is Ub = Rp^/m, 
where R is given in Eq. (69). Therefore one finally obtains 




(72) 
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3 Field fluctuations during inflation 

Apart from explaining the initial condition for the hot Big Bang model, the flatness 
problem, and the horizon problem, cosmological inflation [154, 13, 14] is one of the 
most favored candidate for the origin of structure in the Universe (for reviews on 
inflation, see [15, 155]). There are many models of inflation, but by far the simplest 
is one in which inflation is generated by the large energy density of a scalar field. 
The scalar field driven inflation not only explains the homogeneity and the flatness 
problems but also the observed scale invariance of the density perturbations. 

Inflation based on a scalar field theory is described by the following Lagrangian: 

1 

C=^R+-d,ct>d^ct>-V{(t>), (73) 
where R is the curvature scalar. The energy-momentum tensor reads 

T^. = d^(t>d,4> - \gMd'<P - 9,.uV{(l)) (74) 
so that the energy density and the pressure are given by 

p = r„„ = ^4,' + ^l^^{v4,)' + v{4,), (75) 

One of the initial conditions for inflation is that there must be a homogeneous patch of 
the Universe which is bigger than the size of the Hubble horizon [156] (also supported 
by numerical studies, see [157]). However such a stringent condition can be evaded in 
a chaotic inflation beginning at the Planck scale [158, 159, 15, 169]. More comphcated 
situation can be obtained if there are several fields that participate in inflation; the 
classic example is assisted inflation [161, 162]. 

3.1 Fluctuation spectrum in de Sitter space 

The plane wave solution of a massive scalar field 0(x, t) in a spatially flat Robertson- 
Walker metric can be decomposed into Fourier modes by 

^ (2^^ / '^'^ ('^^^^)"'''^ + ^-'O ■ 
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Solving the Klein-Gordon equation for the scalar field in a conformal metric: ds^ — 
g^ydx^dx^ = a^{T, x){dT'^ — dx'^), the mode function can be given by [163, 164, 165, 166] 

r = -H-'e-"' , and = ^ - ^ , (78) 

where m is the mass of the scalar field, H^^^ and H^^^ are the Hankel functions and 
Ci, C2 are constants. The readers might be tempted to take the limit r <^ 0, in order 
to match the above solution with the plane wave solution in a Minkowski background. 
However this leads to a quasi static de Sitter solution [167]. More technically, it has 
been shown that using a point splitting rcgularization scheme, it is possible to obtain 
a Bunch-Davies vacuum for a de Sitter background which actually corresponds to 
taking ci = 0, and C2 — 1- A simple but intuitive way has been developed in [168], 
where it has been argued that during a de Sitter phase, the main contribution to the 
two point correlation function comes from the long wavelength modes; /c|r| -C 1 or 
k <^ Hexp{Ht). Therefore the two point function is defined by an infrared cutoff 
which is determined by the Hubble expansion [168] 

The result of the integration yields [166, 165, 167, 168] an indefinite increase in the 
variance with time 

(0^)-^^- (80) 

This result can also be obtained by considering the Brownian motion of the scalar field 
[160]. 

For a massive field with m -C i7, and u ^ 3/2, one does not obtain an indefinite 
growth of the variance of the long wavelength fiuctuations, but [166, 165, 167, 168] 

<*^)= 8^(1 (81) 

In the limiting case when m — > if, the variance goes as (0^) ^ H^. In the limit 
m ^ H, the variance goes as (0^) ~ (i7^/127r^m) [167]. Only in a massless case (0^) 
can be treated as a homogeneous background field with a long wavelength mode. This 
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result plays an important role for the rest of this review as it implies that in a de 
Sitter phase any scalar field, including the AD condensate, are subject to quantum 
fluctuations. 



3.2 Slow roll inflation 

A completely flat potential can render inflation eternal, provided the energy density 
stored in the flat direction dominates. The inflaton direction is however not completely 
flat but has a potential with some slope. An inflationary phase is obtained while 

« ^V(^U (82) 

3i/0 ^ (83) 

where prime denotes derivative with respect to 0. In the above the approximations 
are: 0^ < V^(0), and 4> < ^'(0); which lead to the slow roll conditions (see e.g. 
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<4>) = -y[v) ^^^^ 

|r/(0)| = M2^«1. (85) 

Note that e is positive by deflnition. 

These conditions are necessary but not sufficient for inflation. They only constrain 
the shape of the potential but not the velocity of the field 4>. Therefore a tacit as- 
sumption behind the success of the slow roll conditions is that the inflaton fleld should 
not have a large initial velocity. 

Inflation comes to an end when the slow roll conditions are violated, e ~ 1, and 
?7 ~ 1. However, there are certain models where this need not be true, for instance in 
hybrid inflation models [169], where inflation comes to an end via a phase transition, 
or in oscillatory models of inflation where slow roll conditions are satisfled only on 
average [170]. 

One of the salient features of the slow roll inflation is that there exists a late time 
attractor behavior. This means that during inflation the evolution of a scalar fleld at 
a given fleld value has to be independent of the initial conditions. Therefore slow roll 
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inflation should provide an attractor behavior which at late times leads to an identical 
field evolution in the phase space irrespective of the initial conditions [171]. In fact 
the slow roll solution does not give an exact attractor solution to the full equation 
of motion but is nevertheless a fairly good approximation [171]. A similar statement 
has been proven for multi-field exponential potentials without slow roll conditions 
(i.e., assisted infiation) [161]. The attractor behavior of the infiaton leads to powerful 
predictions which can be distinguished from other candidates of galaxy formation [19]. 
The standard definition of the number of e-foldings is given by 

N ^In'-^ = r-'Hdt^'f ld<^, (86) 
a{t) Jt M^J^^„,V' ^' ^ ^ 

where 0end is defined by e{(pend) ^ 1) provided inflation comes to an end via a violation 
of the slow roll conditions. The number of e-foldings can be related to the Hubble 
crossing mode k — akH^ by comparing with the present Hubble length aoi?o- The 
final result is [9] 

Nik)^62 - In-— - In , + In -^^ - - In ^ (87) 

"0^0 Vk Vend Prh 

where the subscripts end {rh) refer to the end of inflation (onset of reheating). The 
details of the thermal history of the Universe determine the precise number of e- 
foldings, but for most practical purposes it is sufficient to assume that N{k) 50, 
keeping all the uncertainties such as the scale of infiation and the end of infiation 
within a margin of 10 e-foldings. A significant modification can take place only if 
there is an epoch of late inflation such as thermal inflation [172], or in theories with a 
low quantum gravity scale [173]. 

3.3 Primordial density perturbations 

Initially, the theory of cosmological perturbations has been developed in the context 
of FRW cosmology [175], and for models of inflation in [176, 177, 178, 179, 180]. For 
a complete review on this topic, see [16]. For a real single scalar field there arise 
only adiabatic density perturbations. In case of several fluctuating flelds there will in 
general also be isocurvature perturbations. We briefly describe the two perturbations 
and their observational differences. 
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3.3.1 Adiabatic perturbations and the Sachs- Wolfe effect 

Let us consider small inhomogeneities 0(x, t) — 4>{t) + d4>{x, t) such that Sep <S 4>. Per- 
turbations in matter densities automatically induce perturbations in the background 
metric, but the separation between the background metric and a perturbed one is not 
unique. One needs to choose a gauge. A simple choice would be to fix the observer to 
the unperturbed matter particles, where the observer will detect a velocity of matter 
field falling under gravity; this is known as the Newtonian or the longitudinal gauge 
because the observer in the Newtonian gravity limit measures the gravitational poten- 
tial well where matter is falling in and clumping. The induced metric can be written 
as 



[1 + 2^)dT^ - (1 - 2^)Sikdx'dx''\ , 

where $ has a complete analogue of Newtonian gravitational potential. In the case 
when the spatial part of the energy momentum tensor is diagonal, i.e. 5Tj = Sj, it 
follows that $ = ^' [16]. Right at the time of horizon crossing one finds a solution for 
S4> as 

im')-^^ (89) 



2k^ 

where denotes the instance of horizon crossing. Correspondingly, we can also define 
a power spectrum 



27r2 
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(90) 



k=aH 



Note that the phase of 5(j)k can be arbitrary, and therefore, infiation has generated a 
Gaussian perturbation. 

In the limit /c ^ 0, one can find an exact solution for the long wavelength inhomo- 
geneities k -C aH [181, 16], which reads 

/I /•* \' H 
^ ci - / adt'] +C2— , (91) 
\a JO / a 

^ = '-Lfadt'-c], (92) 

(j) a \ Jo ) 

where the dot denotes derivative with respect to physical time. The growing solutions 
are proportional to Ci, the decaying proportional to C2. Concentrating upon the grow- 
ing solution, it is possible to obtain a leading order term in an expansion with the help 
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of the slow roll conditions: 

f - 

Note that at the end of inflation, which is indicated by a = 0, or equivalently by 
H — —H'^, one obtains a constant Newtonian potential ~ Ci. This is perhaps the 
most significant result for a single field perturbation. 

In a long wavelength hmit one obtains a constant of motion ( [180, 182, 16] defined 

as 

C=o + w = -. 95 

6 1 + W p 

If the equation of state for matter remains constant there is a simple relationship which 
connects the metric perturbations at two different times [180, 182, 16] 

Mtf) - TT^TT— 4r^^^(^0 ■ (96) 

The comoving curvature perturbation [183] reads in the longitudinal gauge [16] for 
the slow roll inflation as 

= - ^ [H-^^k + ^k) ■ (97) 

For CMB and structure formation we need to know the metric perturbation during the 
matter dominated era when the metric perturbation is ^{tf) ~ (3/5)ci. Substituting 
the value of ci from Eq. (94), we obtain 



(98) 

k=aH 



In a similar way it is also possible to show that the comoving curvature perturbations 
is given by 

(99) 



k=aH 

where denotes the field perturbation on a spatially flat hypersurfaces, because on a 
comoving hypersurface 5(p — 0, by definition. Therefore, on flat hypersurfaces 

S(f)k^^St, (100) 
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where 5t is the time displacement going from flat to comoving hypersurfaces [19, 9]. 
As a result 

nk = H5t. (101) 

Note that during matter dominated era the curvature perturbation and the metric 
perturbations are related to each other 

$fe = -^7^fe. (102) 
o 

In the matter dominated era the photon sees this potential well created by the 
primordial fluctuation and the redshift in the emitted photon is given by 

^--^k. (103) 

At the same time, the proper time scale inside the fluctuation becomes slower by an 
amount 5t/t — Therefore, for the scale factor a oc t^/^, decouphng occurs earlier 
with 

a 6 t 6 

By virtue of T oc this results in a temperature which is hotter by 

+ (105) 

This is the celebrated Sachs- Wolfe effect [184] , which we shall revisit when discussing 
isocurvature fluctuations. 

3.3.2 Spectrum of adiabatic perturbations 

Now, one can immediately calculate the spectrum of the metric perturbations. For a 
critical density Universe 



P 



k 

where — > — A;^, in the Fourier domain. Therefore, with the help of Eqs. (90,98), 
one obtains 
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where the right hand side can be evaluated at the time of horizon exit k — aH. In fact 
the above expression can also be expressed in terms of curvature perturbations [19, 9] 

h = (108) 

and following Eq. (97), we obtain 51 = A/25VTi{k) = {A/2b){H / 4)Y{H /inf , exactly 
the same expression as in Eq. (107). With the help of the slow roll equation 3if0 = 
—V, and the critical density formula 3H^Mp = V, one obtains 



(109) 



where we have used the slow roll parameter e = {M'^/2){y' /VY . The COBE satellite 
measured the CMB anisotropy and fixes the normalization of 5^{k) on a very large 
scale. For a critical density Universe, if we assume that the primordial spectrum can 
be approximated by a power law and ignoring gravitational waves: 



k 



(n-l)/2 



5^{k) = 1.91 X 10-' , (110) 

\ l^pivot ) 

where n is the spectral index and k^i^ot — 7.5aoi?o is the scale at which the normal- 
ization is independent of the spectral index. 
The spectral index n(A;) is defined as 

^ (ilnP$ N 

This definition is equivalent to the power law behavior if niki) is fairly a constant 
quantity over a range of k of interest. The power spectrum can then be written as 

7'$(A;) oc A;"-^ (112) 

If n = 1, the spectrum is flat and known as Harrison-Zeldovich spectrum [185]. For 
n 7^ 1, the spectrum is tilted and n > 1 is known as blue spectrum. In terms of the 
slow roll parameters, one can write [171] 

2er,-Ae\ ^ -2ev + e St ^ '^^^ + Ve + (113) 



dlak dink dink 

where 



' =^P V'2 ' ^ =^P Y3 ■ (114) 

41 



Thus one finds[186] 

n-l^-6e + 2r]. (115) 

Slow roll inflation requires that e <S 1, 177I <^ 1, and therefore naturally predicts 
small variation in the spectral index within A In A; ^ 1. The recent Boomerang data 
suggest [10] 

|n-l|<0.1. (116) 



The rate of change in 77 is also very small, and can be estimated in a similar way [187] 

= -16677 + 246^ + 2e . (117) 

dink 



It is possible to extend the calculation of metric perturbation beyond the slow roll 
approximation basing on a formahsm similar to that developed in [188, 189, 190, 191]. 

3.3.3 Gravitational waves 

Gravitational waves are linearized tensor perturbations of the metric and do not couple 
to the energy momentum tensor. Therefore, they do not give rise a gravitational 
instability, but carry the underlying geometric structure of the space-time. The first 

calculation of the gravitational wave production was made in [192], and the topic 
has been considered by many authors [193]. For reviews on gravitational waves, see 
[16, 194]. 

The gravitational wave perturbations are described by a line element ds^ + Sds'^, 
where 

ds^ = a^{T){dT'^ — dx'dxi) , 5ds^ = —a^{T)hijdx'dx^ . (US) 

The gauge invariant and conformally invariant 3-tensor hij is symmetric, traceless 
S'^^hij — 0, and divergenceless Vihij — (Vj is a covariant derivative). Massless spin 2 

gravitons have two degrees of freedom and as a result are also transverse. This means 
that in a Fourier domain the gravitational wave has a form 

= h+e^ + h^e^^ . (119) 
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For the Einstein gravity, the gravitational wave equation of motion follows that of a 
massless Klein Gordon equation [16] . Especially, for a flat Universe 

^5 + 3^^+(^)^} = 0' (120) 

As any massless field, the gravitational waves also feel the quantum fluctuations in an 
expanding background. The spectrum mimics that of Eq. (90) 



(121) 

k=aH 



Note that the spectrum has a Planck mass suppression, which suggests that the ampli- 
tude of the gravitational waves is smaller compared to that of the adiabatic perturba- 
tions. Therefore it is usually assumed that their contribution to the CMB anisotropy 
is small. The corresponding spectral index can be calculated as [186] 

din A; ~ ^ ' 

Note that the spectral index is negative. 

3.4 Multi- field perturbations 

In multi-field inflation models contributions to the density perturbations come from 
all the flelds. However unlike in a single scalar case, in the multi-fleld case there 
might not be a unique late time trajectory corresponding to all the flelds. This is 
true in particular for those flelds that are effectively massless during inflation, such 
as the MSSM flat direction flelds. Therefore, in these cases scalar perturbations will 
depend on the field trajectories and thus on the choice of initial conditions, with an 
ensuing loss of predictivity. In a very few cases it is possible to obtain a late time 
attractor behavior of all the flelds; an example is assisted inflation [161]. Let us here 
nevertheless assume that there is an underlying unique late time trajectory resulting 
in a simple expression for the amplitude of the density perturbations and the spectral 
index [195, 155]. 
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3.4.1 Adiabatic and isocurvature conditions 

There are only two kinds of perturbations that can be generated. The first one is 
the adiabatic perturbation discussed previously; it is a perturbation along the late 
time classical trajectories of the scalar fields during infiation. When the primordial 
perturbations enter our horizon they perturb the matter density with a generic adia- 
batic condition, which is satisfied when the density contrast of the individual species 
is related to the total density contrast 5^ 




(123) 



where h stands for baryons, c for cold dark matter, 7 for photons and v for neutrinos. 

The other type is the isocurvature perturbation. During infiation this can be viewed 
as a perturbation orthogonal to the unique late time classical trajectory. Therefore, if 
there were N fiuctuating scalar fields during infiation, there would he N — 1 degrees 
of freedom which would contribute to the isocurvature perturbation. 

The isocurvature condition is known as 5p = 0: the sum total of all the energy 
contrasts must be zero. The most general density perturbations is then given by a 
linear combination of an adiabatic and an isocurvature density perturbations. 

3.4.2 Adiabatic perturbations due to multi-field 

In a comoving gauge Eq. (97) with TZ — —HS4>/4> holds good even for multi-field 
infiation models, provided we identify each field component of along the slow roll di- 
rection. There also exists a relationship between the comoving curvature perturbations 
and the number of e-foldings N [181, 196, 195, 155] 

dN 

n^5N=—5(t>a, (124) 

0(Pa 

where N is measured by a comoving observer while passing from fiat hypersurface 
(which defines Scf)) to the comoving hypersurface (which determines TZ) [195, 197]. 
The repeated indices are summed over and the subscript a denotes a component of 
the infiaton. A more intuitive derivation has been given in [155, 9]. 
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If again one assumes that the perturbations in 5(f)a have random phases with an 
amphtude {H/2nY, one obtains 

V dN dN 



^ 757r2 Ml 



(125) 



p ' 



For a single component dN/dcf) = {M^'^V/V), and then Eq. (125) reduces to Eq. (109). 
By using slow roll equations we can again define the spectral index 

where V^a = dV/d(t)a, and similarly N^^ = dN/d(f)a. For a single component we recover 
Eq. (115) from Eq. (126). These results prove useful in constraining the AD potential 
by cosmological density perturbations, as will be discussed in Sect. 5.3. 

3.4.3 Isocurvature perturbations and CMB 

One may of course simply assume a purely isocurvature initial condition. For any 
species the entropy perturbation is defined by 

5. = ^-^, (127) 

Thus, if initially there is a radiation bath with a common radiation density contrast 
Sr, a baryon-density contrast Sb = 35r/4, and a CDM density contrast Sc, then 

^ ^ _ 3^^ ^ p,6p, {?>/A)pM ^ PL±SmPlSp^ ^ S,, (128) 

4 PrPc PrPc 



where we have used the isocurvature condition 5pr + Spc — 0, and the last equality 
holds in a radiation dominated Universe. 

However a pure isocurvature perturbation gives five times larger contribution to 
the Sachs- Wolfe effect compared to the adiabatic case [198, 19, 9]. This result can be 
derived very easily in a matter dominated era with an isocurvature condition 5pc — 
—Spr, which gives a contribution TZk — (1/3)5"^. Therefore from Eqs. (102,105), we 
obtain AT^/T = —3^/15. There is an additional contribution from radiation because 
we are in a matter dominated era, see Eq. (128), S* ~ 5c = "(3/4)5^.. The sum total 
isocurvature perturbation AT^/T = —S/15 — S/3 = —63/ 15, where 3 is measured on 
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the last scattering surface. The Sachs- Wolfe effect for isocurvature perturbations fixes 
the slope of the perturbations, rather than the amplitude [199, 200]. Present CMB 
data rules out pure isocurvature perturbation spectrum [201, 202], although a mixture 
of adiabatic and isocurvature perturbations remains a possibility [201, 203, 204, 205]. 
In the latter case it has been argued that the adiabatic and isocurvature perturbations 
might naturally turn out to be correlated [206, 207]. The most general power spectrum 
is not a single function but a 5 x 5 matrix, which contains all possible adiabatic and 
isocurvature perturbations together with their cross-correlations. As discussed in [207], 
resolving the perturbation spectrum in all its generality would be an observational 
challenge that probably would have to wait for the determination of the polarization 
spectrum by the Planck Surveyor Mission. 

It is sometimes useful to consider the ratio a of the total power spectra, defined as 
Vtot = Vad + Visa [208, 204], where a is defined as [204] 



a 











k=aH 


7 



(129) 



where 81^ is the perturbation in the photon energy density due to isocurvature pertur- 
bations and 5" is the perturbation due to adiabatic perturbations. 

3.5 Inflation models 

A detailed account on infiation model building can be found in many reviews [15, 98, 
155, 9]. Here we briefiy recall some of the popular models with a particular emphasis on 
supersymmetric infiation. First we recapitulate some aspects of non-supersymmetric 
models. 

3.5.1 Non-supersymmetric inflation 

The very first attempt to build an infiation model was made in [154], where one loop 
quantum correction to the energy momentum tensor due to the space-time curvature 
were taken into account, resulting in terms of higher order in curvature invariants. 
Such corrections to the Einstein equation admit a de Sitter solution [209], which was 
presented in [154, 210]. Infiation in Einstein gravity with an additional E? term was 
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considered in [211] (for a discussion of inflation in pure i?^ gravity, see [212]). Such 
a theory is conformally equivalent to a theory with a canonical gravity [213] with 
a scalar field having a potential term. A similar situation arises in theories with a 
variable Planck mass, i.e., in scalar tensor theories [66]. Inflation in these models has 
been studied extensively [214]. 

The simplest single field inflation model is arguably chaotic inflation [158, 215] 
with a generic potential 

y = — ^</'" , (130) 

where a is a positive even integer. In chaotic inflation slow roll takes place for ^ Mp, 
and the two slow roll parameters are given by [9] 

6=y^, 77 = a(a-l)— . (131) 

Inflation ends when e = 1, or, aMp. The cosmological scales leave the horizon 



when (f) = \/2NaMp, and the spectral indices for scalar and tensor perturbations turn 

out to be [9] 

2 + a 3.1a 
n = 1 , r = . 132 

The amplitude of the density perturbations, if normalized at the COBE scale, yields 
the constraint A ~ 4 x 10~^^. 

An exponential potential, such as might arise in string theories and theories with 
extra dimensions, 

yW = y„exp(-^-^). (133) 

would give rise to a power law a{t) oc for the scale factor, so that inflation occurs 
when p > 1. Multiple exponentials with differing slopes give rise to what has been 
dubbed as assisted inflation [161]. 

3.5.2 F-term inflation 

In four dimensions the N — 1 supersymmetric potential receives two contributions: 
one from the F-term, which is related to the chiral supermultiplets, and the second 
from the D-term, which contains the gauge interactions. For a detailed discussion of 
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supersymmetric inflation we refer to the review by Lyth and Riotto [155]. Here we 
give a brief resume of the two types of inflation. 

Historically, supersymmetric inflation was first introduced to cure some of the prob- 
lems associated with the fine tuning of new inflation [216], but since then utilizing 
supersymmetry as a tool for inflation has gained in popularity (we describe supersym- 
metry in Sect. 4.2, and for supergravity, see Sect. 4.5.2.). The F-term potential can 
be derived from the superpotential W 

y(0, 0*) = F«F, , F,^-(^^^ , (134) 

where for renormalizable interactions W has a mass dimension three. 

Supersymmetry is broken whenever |Fp 7^ 0. A simple working example is to 
consider the superpotential 

W - XS{(I)^ - 4) (135) 

which is invariant under a global R symmetry with the superflelds S and N carrying 

respectively the R charges 1 and 0. The scalar components of these superflelds can be 
written in the form 

where we have used i?-transformation in order to make S real. The potential follows 
from Eq. (134): 

V = - \'<pM - 4>l) + ^{<Pi + 4>lf + AV2(0? + <Pl) . (137) 

The supersymmetric vacuum is located at cr = 0, 0i — (t)o, and 02 — 0. Note that 
the potential has a flat direction along a-axis when a > (Jinst = 0o- When a < 
(Tinst-i the mass squared of 0i becomes negative and suggests a phase transition along 
the 01 direction. When this happens cr, 0i, and 02 begin to oscillate around their 
supersymmetry preserving vacua. If (pi — (f)2 — 0, the height of the potential is 
given by = A^0O) and as a consequence one obtains a period of inflation. This is the 
simplest example of a flat direction giving rise to an inflation potential, and it is known 
as the hybrid model, flrst described in a non-supersymmetric context in [215, 160] and 
in a supersymmetric context in [217]. 
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In order to have a graceful exit from inflation one requires a slope for the flat 
direction such that u can roll down and approach Omst- The flatness of the potential 
can be lifted in two ways: by radiative corrections [218], or by the low energy soft 
supersymmetry breaking effects. 

Due to supersymmetry breaking the fermions obtain a mass of the order {d'^Wj 90^) = 

while the two complex scalars receive a mass squared }^S'^ ± A^^q- The one- loop 

radiative correction to the potential is given by [219] 

1 A/f2 
*V=^E(-)''M/ln-^. (138) 

where /j denotes the number of fermions, is the fermion mass squared, and M the 
cut-off or the renormalization scale. The summation should be taken over all helicity 
states i. In the present example the effective potential along the fiat direction is given 
by 

y = AVj(l + ^^ln^). (139) 

where C is a constant essentially counting the states running in the loops. If the loop 
correction dominates over the tree level potential, there is a period of inflation which 
typically ends when 

The COBE normalization sets the scale of inflation to 

^ ci/^x X 10^^ GeV (141) 

while the spectral index is given by 

Therefore, depending on the coupling constant A and the number of e-foldings N, it 
is possible to have a wide range of inflation energy scales which all provide a spectral 
index n ~ 0.96 - 0.98. 

Soft supersymmetry breaking contributions induce ~ 0{TcV). One could also 
imagine that the mass of a appears dynamically if a has couplings to bosons and 
fermions; these may induce a typical running mass cx: a'^\n{a/M) [220, 155, 221]. 
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3.5.3 D-term inflation 

In the above discussion we have neglected the gauge contribution. The D-term — 
—ga{(f)*T"'(j)) gives rise to a scalar potential (see [37, 34]) 

n0>*) = ^E^"^a (143) 

where (T")i ^ satisfy [T",r''] = i/"^'=T'= {f'"' is the structure constant). 

The simplest realization of D-term inflation reproduces the hybrid potential with 
three chiral superfields, S, 0+, and 0_ with (non-anomalous) C/(l) charges 0, — 1 
[222]. The superpotential can be written as 

W = A,S(/)+(/)_ . (144) 

The scalar potential then reads [222] 

V = x'\s\' + 10-r) + A^i0+0-r + y (i<^+r - i<^-r + ef , (145) 

where g is the gauge coupling and ^ is the Fayet-Iliopoulos D-term. Note that the 
potential allows unique supersymmetry preserving vacua with a broken gauge symme- 
try S — — 0, and 0_ = ^. By virtue of the coupling, when \S\ > Sinst — gC/\ 
the fields — > 0, and therefore inflation occurs because of the Fayet-Iliopoulos 
D-term V = The slope along the infiaton direction S can be generated by the 

one-loop contribution and reads 

Inflation ends when slow roll condition breaks down for S ^ ((//27r^/2)Mp, and the 
predictions for the inflationary parameters are similar to the previous discussion. D- 
term inflation based on an anomalous U{1) symmetry (which could appear in string 
theory [223]) is no different. 

Hybrid inflation is successful but has also problems that are related to the initial 
conditions. In [224], it was pointed out that hybrid inflation requires an extremely 
homogeneous field configuration for the fields orthogonal to the infiaton. In our ex- 
ample the orthogonal fields to the infiaton must be set to zero with a high accuracy 
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over a region much larger than the initial size of the horizon. It is possible to solve 
this impasse by having a pre-infiationary matter dominated phase when the field or- 
thogonal to the inflaton direction oscillates and decays into lighter degrees of freedom, 
gradually settling down to the bottom of its potential [224]. 

3.5.4 Supergravity corrections 

When the field values are close to the Planck scale, supergravity (SUGRA) effects be- 
come important and may ruin the fiatness of the infiaton potential. The soft breaking 
mass of the scalar fields are typically [225, 217, 226, 45, 46] 

"^'on-^-om^ (147) 

Once the infiaton gains a mass ~ H, the field simply rolls down to the minimum of 
the potential and inflation stops. Indeed, in SUGRA the slow roll parameter 

N^M^^-%^-0(l), (148) 

where rrvsuGRA ~ "^lc/5y + ^susy/^M^) ~ rn^susY + 0{1)H'^. Note that the latter 
contribution dominates in an expanding Universe and violates the slow roll condition. 
For field values smaller than Planck scale it is always possible to obtain e <^ 1, but in 
supergravity t] can never be made less than one for a single chiral field with a minimal 
kinetic term. This is known as the r] problem in SUGRA models of inflation [217]. 

When there are more than one chiral superfields, it might be possible to cancel the 
dominant 0{1)H correction to the infiaton mass by choosing an appropriate Kahler 
term [226, 217] (see also discussion in [155]). In hybrid infiation models derived from an 
F-term the dominant 0{1)H correction in the mass term can be canceled if |iV| = 
exactly, which however seems to lead to an initial condition problem, as discussed 
above. The fact that the superpotential is hnear in S in Eqs. (135,144) guarantees 
the cancellation of the dominant contribution in the mass term for a minimal Kahler 
term ~ \S\'^. For non- minimal Kahler potential such as K = \S\'^ + I3\S\'^/M'^ + 
one obtains (d'^K/dSdS*)^^ 1 — A(3\S\'^/Mp. These contributions again lead to 
a problematic (3 x 0{1)H contribution to the infiaton mass unless the value of the 
unknown constant /3 is suppressed. 
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In [222], it was shown that the rj problem does not appear for D-term inflation 
even for the minimal Kahler potential because the main contribution to the inflation 
potential does not come from the vcv of the inflaton field alone, but from the Fayet- 
Iliopoulos term which belongs to the D-sector of the potential. Based on this fact 
many D-term inflation models have been written down [227, 228]. Therefore, hybrid 
inflation, whether realized as an effective potential coming from F-sector or from D- 
sector, appears to be among the most promising models for supersymmetric inflation. 

3.6 Reheating of the Universe 
3.6.1 Perturbative inflaton decay 

Traditionally reheating has been assumed to be a consequence of the perturbative de- 
cay of the inflaton [229, 198, 22]. After the end of inflation, when H < m^^ the inflaton 
field oscillates about the minimum of the potential. Averaging over one oscillation re- 
sults in [230] pressureless equation of state where (p) = (0^/2 — V^cj))) vanishes^, so 
that the energy density redshifts as during matter domination with — pi{ai/af 
(subscript i denotes the quantities right after the end of inflation). If F^ represents 
the decay width of the inflaton to a pair of fermions, then the inflaton decays when 



H{a) = Y (l/3Mp)pj(ai/a)^/^ ~ F^. When the inflaton decays, it releases its energy 
into the thermal bath of relativistic particles whose energy density is determined by 
the reheat temperature T^^, given by 



However the inflaton might not decay instantaneously. In such a case there might 
already exist a thermal plasma of some relativistic species at a temperature higher 
than the reheat temperature already before the end of reheating [22]. If the inflaton 
decays with a rate F<^, then the instantaneous plasma temperature is found to be [22] 





(149) 




(150) 



^This will be discussed in a more detail in Sect. 5.8. 
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where denotes the effective relativistic degrees of freedom in the plasma. The 
temperatures reaches its maximum T^ax soon after the infiaton field starts oscillating 
around the minimum. Once the maximum temperature is reached, then ~ a~^/^, 
and T ~ a~^/^ until reheating and thermalization is completely over [231, 232, 233, 
234, 235]. 

The process of thermalization has two aspects; achieving kinetic equilibrium, and 
achieving chemical equilibrium. Kinetic equilibrium can be reached by 2 ^ 2 scatter- 
ing and annihilation. For chemical equilibrium one requires particle number changing 
interactions such as 2 ^ 3 processes. In [232], soft processes which allow for small 
momentum transfer with a larger cross-section have been advocated for chemical equi- 
libration, while in [231], hard processes have been invoked. Therefore, depending on 
the interactions, thermahzation time scale could be short, such as in the case of soft 
scattering processes, or, it could be long compared to the Hubble time if only hard 
processes are operative. Recently it has been argued [233, 234, 235], that thermaliza- 
tion time scale can be as long as the time it takes for the infiaton decay products with 
typical energies C(m^) to lose the energy ~ (m^ — Tru)- The main conclusion is that 
inelastic scattering interactions 2^3 can thermalize the Universe faster compared to 
elastic interaction 2-^2. Inelastic interactions can achieve the kinetic and chemical 
equilibrium both, and therefore, F^^^ could be considered as the true thermalization 
time scale. In [235], the authors have studied thermalization in QCD based approach 
with emphasis upon late thermalization and hadronization. 

3.6.2 Non-perturbative infiaton decay 

Much effort has lately been devoted to non-perturbative effects which are essentially 
non-thermal. These may lead to a rapid transfer of the infiaton energy to other degrees 
of freedom by the process known as preheating. The requirement is that the infiaton 
quanta couple to other (essentially masslcss) fields x through e.g. terms like (f^x^. 
The quantum modes of x may then be excited during the infiaton oscillations via a 
parametric resonance. Preheating has been treated both analytically [236, 237, 238, 
239, 240, 241, 242, 243], and on lattice [244]. 
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Like bosons, fermions can also be excited through preheating [245, 246, 136, 247]. 
In fact it has been argued that fermionic preheating is perhaps more effective than 
bosonic preheating [246, 136, 248]. However note that supersymmetry is effectively 
broken during the inflaton oscillations [247, 249, 250] . As a consequence, one naturally 
expects corrections to the inflaton potential during the oscillations [251]. Therefore 
in most supersymmetric models of inflation preheating might not turn out to be very 
relevant. 

3.6.3 Gravitino and infiatino problems 

The reheat temperature should certainly be above the BBN temperature T > 0{1) MeV, 
but there also exists an upper bound from gravitino overproduction. In supergravity 
the supcrpartner of the graviton is a spin-3/2 gravitino, which gets a mass from the 
super-Higgs mechanism [252] when supersymmetry is spontaneously broken. Typically 
supergravity is broken in a hidden sector by some non-perturbative dynamics. Super- 
symmetry breaking is then mediated via gravitational (or possibly other) couplings to 
the observable sector in such a way that sfermions and gauginos get masses of order 
electroweak scale [253, 34]. In addition, the gravitino also gets a mass which in the 
simplest gravity mediated models is of order 1 TeV [254] (see Sect. 4.4). 

If the gravitino is not the hghtest supersymmetric particle (LSP), it will decay. 
Gravitino has two hehcity states ±3/2 and ±1/2. The latter one is mainly the goldstino 
mode which is eaten by the super-Higgs mechanism. The goldstino coupling strength is 
inversely proportional to the momentum, so that at low energies the gravitino coupling 
is mainly dictated by the goldstino mode [255]. At temperatures much above the 
sparticle masses, it is the massless ±3/2 mode that governs the gravitino interactions. 
The helicity ±3/2 mode can decay into gauge bosons and gauginos through a dimension 
5-operator with a lifetime 

AMI , ^ 

r3/2^A,A~-^- (151) 

"^3/2 

Typically r ~ 10^ — 10^ s for a gravitino mass in the range 10 TeV < m3/2 < 100 GeV. 

Although the gravitino interactions with matter are suppressed by the Planck mass, 
they can be generated in great abundances very close to the Planck scale [256]. In- 
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flation would dilute their number density [257], but during reheating they would be 
regenerated though scattering of gauge and gaugino quanta, with adverse consequences 
[102, 258, 259, 260, 261, 262, 263, 264]. The resulting gravitino abundance has been 
estimated to be [102] 



""'I' - 2.4 X 10-1^ ( -4^) \l - 0.018 In ( 

Vio^GevyL Vio^Gevy 



Trh \ 



where s defines the entropy density and T^h denotes the reheating temperature of the 
Universe. The abundance Eq. (152) could well be increased by an order of magnitude 
if gravitino interactions with other chiral multiplets arc included [260]. In [261] it was 
argued that at finite temperatures gravitino overproduction could be enhanced, but 
the calculation was criticized in [262, 263]; for a recent discussion on this topic, see 
[264]. 

Since the gravitino is a late decaying particle, BBN yields a restriction on the 
reheat temperature [258, 2]. For instance, gravitino decay products can enhance the 
abundance of D +^ He due to photo fission of '^He which implies [2] 

^ < (10-14-10-11)^ 
s 

Trh < (10^ - 101°) GeV , 100 GeV < mg/a < 10 TeV . (153) 
The constraint on the reheating temperature is [2] 

for m3/2 < 1.6 TeV. 

In gauge mediated supersymmetry breaking scenarios, to be discussed in Sect. 4.4., 
the gravitino can have a very light mass ~ 10-^ GcV [265] and can be a hot dark matter 
candidate [266]. Small (or large) gravitino masses can also be obtained in SUGRA 
models with non-minimal Kahler terms, such as the no-scale model [267] . In anomaly 
mediation the gravitino mass is large with 7713/2 ~ mgoft/ot ^ rusoft [268]. In general, 
if the gravitino is not LSP and heavier than 10 TeV, it decays before nucleosynthesis 
and thus does not cause any cosmological problems [269]. 

Gravitinos could also be produced by non-perturbative processes, as was first de- 
scribed in [247], where the formalism for exciting the helicity ±3/2 component of the 
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gravitino was developed. Later the production of the hehcity ±1/2 state, which for a 
single chiral multiplet is the supcrpartner of the inflaton known as inflatino, has been 
studied by several authors [249, 250, 248, 270, 271]. The decay channels of the inflatino 
have been first discussed in [272, 273]. It has been suggested [272] and also explicitly 
shown [271] that in realistic models with several chiral multiplets, the helicity ±1/2 
gravitino production is not a problem for nucleosynthesis as long as the inflationary 
scale is sufficiently higher than the scale of supcrsymmetry breaking in the hidden 
sector and the two sectors are gravitationally coupled. A very late decay of inflatino 
could however be possible, as argued in [273, 274]. In [274], it was argued that if 
the inflatino and gravitino were not LSP, then late off-shell inflatino and gravitino 
mediated decays of heavy relics could be significant. 
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4 Flat directions 



4.1 Degenerate vacua 

At the level of renormalizable terms, supersymmetric field theories generically have 
infinitely degenerate vacua. This is a consequence of the supersymmetry and the 
gauge symmetries (and discrete symmetries such as i?-parity) of the Lagrangian, which 
allow for certain types of interaction terms only. Therefore, in general there are a 
number of directions in the space of scalar fields, collectively called the moduli space, 
where the scalar potential is identically zero. In low energy supersymmetric theories 
such classical degeneracy is accidental and is protected from perturbative quantum 
corrections by a non-renormalization theorem [275] . In principle the degeneracies could 
be lifted by non-perturbative effects. However such effects are likely to be suppressed 
exponentially and thus unimportant because all the couplings of low energy theories 
are typically weak even at relatively large vevs. Therefore in the supersymmetric limit 
when Mp oo, the potential for the flat direction always vanishes. 

In the MSSM the moduli fields are quark, lepton and Higgs chiral fields. In string 
theories there are often additional moduli fields associated with the conformal field 
theory degrees of freedom and world sheet discreet i?-symmetries [276]. The moduh 
space of string theory can also be lifted by a soft supersymmetry breaking masses of the 
order of the gravitino mass 777,3/2. Since the moduli interactions with others fields are 
usually Planck mass suppressed, the string moduli are also a cause for worry because 
they may decay after nucleosynthesis. This problem has been dubbed as the moduli 
problem [225, 277]. However, the MSSM fiat directions are made up of condensates of 
squarks, Higgses, and sleptons, and can evaporate much before nucleosynthesis. 

However there is an effective potential for the fiat direction condensate fields which 
arises as a result of supersymmetry breaking terms and higher dimensional operators 
in the superpotential. In this sense the MSSM fiat directions are only approximately 
fiat at vevs larger than the supersymmetry breaking scale. 
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4.2 MSSM and its potential 

Let us remind the reader that the matter fields of MSSM are chiral superfields $ = 
(f) + \/29'^ + 99F, which describe a scalar (f), a fermion and a scalar auxihary field 
F. In addition to the usual quark and lepton superfields, MSSM has two Higgs fields, 
Hu and Hd- Two Higgses are needed because H'^ ^ which in the Standard Model gives 
masses to the ti-quarks, is forbidden in the superpotential. 
The superpotential for the MSSM is given by [34] 

Wmssm = XuQHuU + XdQHdd + XeLHdC + iiH^Hd , (155) 

where H^, Hd,Q, L,u,d,e in Eq. (155) are chiral superfields, and the dimensionless 
Yukawa couplings A„, A^, Ag are 3x3 matrices in the family space. We have sup- 
pressed the gauge and family indices. Unbarred fields are SU{2) doublets, barred 
fields SU (2) singlets. The last term is the term, which is a supersymmetric version 
of the SM Higgs boson mass. Terms proportional to H*Hu or H^Hd arc forbidden in 
the superpotential, since Wmssm must be analytic in the chiral fields. and Hd are 
required not only because they give masses to all the quarks and leptons, but also for 
the cancellation of gauge anomalies. The Yukawa matrices determine the masses and 
CKM mixing angles of the ordinary quarks and leptons through the neutral compo- 
nents of Hu = {H^,H^) and Hd = {H^Hj). Since the top quark, bottom quark and 
tau lepton are the heaviest fermions in the SM, we assume that only the (3, 3) element 
of the matrices A„, A^, Ag are important. In this limit only the third family and the 
Higgs fields contribute to the MSSM superpotential. 

The SUSY scalar potential V is the sum of the F- and D-terms and reads 

V ^Y.\Fi? + \Y.9lD^D'^ (156) 

i ^ a 

where 

p^^ dWMSSM j^a^^ij.a^_ (^57) 

d(j)i 

Here we have assumed that 4>i transforms under a gauge group G with the generators 
of the Lie algebra given by T". 
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The ijl term provides masses to the Higgsinos 

C D -pi{H^Hj - HlH'a) + c.c , (158) 
and contributes to the Higgs {massY terms in the scalar potential through 

-£ D y D HWny + \h:\' + + 1^^ 1') ■ (i59) 

Note that Eq. (159) is positive definite. Therefore, it cannot lead to electroweak 
symmetry breaking without including supersymmetry breaking {massY soft terms for 
the Higgs fields, which can be negative. Hence, should almost cancel the negative 
soft {massY term in order to allow for a Higgs vev of order ~ 174 GeV. That the 
two different sources of masses should be precisely of same order is a puzzle for which 
many solutions has been suggested [278, 279, 280, 281]. 

Note also that Eq. (155) is the minimal superpotential because we have not included 
terms which are gauge invariant and analytic in the chiral superfields but which violate 
either baryon number B or Icpton number L. The most general gauge invariant and 
renormalizable superpotential would not only include Eq. (155), but also the terms 

Wal=i = ^A^^'=L,L,e, + A'^^'%Q,-4 + Ai''L,i7^, (160) 

Wab=i = ]^>^"^%djdk, (161) 

where i = 1, 2, 3 represents the family indices. The chiral supermultiplets carry baryon 
number assignments S = +1/3 for Qj, 5 = —1/3 for ui, di, and B — for all others. 
The total lepton number assignments are L — +1 for Lj, L — —1 for e^, and L = for 
all the others. The terms in Eq. (160) violate lepton number by one unit, while those 
in Eq. (161) violate baryon number by one unit. 

Unless A' and A" terms are very much suppressed, one would obtain rapid proton 
decay which violates both B and L by one unit. Many other processes also give rise to 
violation in baryon and lepton number (for a review, see [74]). Therefore, there must 
be a symmetry forbidding the terms in Eqs. (160,161), while allowing for the terms in 
Eq. (155). The symmetry is known as i?-parity [255], which is a discrete parity defined 
for each particle as 
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with Pr — +1 for the SM particles and the Higgs bosons, while Pr — —1 for all the 
sleptons, squarks, gauginos, and Higgsinos. Here s is spin of the particle. Without 
the product (—1)^*, the expression is known as matter parity [282], and denoted by 
Pm = (— l)^'-'^"^-'. The quantity (—1)^* is equal to 1 whenever conservation of angular 
momentum holds at a given vertex. In this case matter parity and i?-parity are 
equivalent. If i?-parity is conserved then there will be no mixing between the sparticles 
and the ones which have Pr — +1. This completely forbids potentially dangerous 
terms in Eqs. (160,161). 

Matter parity is actually a discrete subgroup of the continuous U{1)b-l group. 
Therefore, if a gauged U{1)b-l is broken by scalar vevs which carry even integer val- 
ues of 3(S — L), then Pm survives as an exactly conserved discrete remnant [283]. 
Besides forbidding B and L violation from the renormalizable interactions, i?-parity 
has interesting phenomcnological and cosmological consequences. The lightest sparti- 
cle with Pr — —1, the LSP, must be absolutely stable. If electrically neutral, the LSP 
is a natural candidate for non-baryonic dark matter [284, 285]. It may be possible to 
produce LSPs in a next generation collider experiments. 

4.2.1 F-and D-renormalizable flat directions of MSSM 

For a general supersymmetric model with N chiral superfields Xi, it is possible to 
find out the directions where the potential Eq. (156) vanishes identically by solving 
simultaneously 

dW 

= X^T^X = , Fx. = = . (163) 

oXi 

Field configurations obeying Eq. (163) are called respectively D-fiat and F-flat. 

D-fiat directions are parameterized by gauge invariant monomials of the chiral 
superfields. A powerful tool for finding the flat directions has been developed in 
[286, 287, 288, 46, 289, 43], where the correspondence between gauge invariance and 
flat directions has been employed. The conflguration space of the scalar flelds of the 
MSSM contains 49 complex dimensions (18 for Qi, 9 each for Ui and Jj, 6 for Lj, 3 for 
ej, and 2 each for and i?^), out of which there are 12 real D-term constraints (8 for 
SU (3)c 3 for SU (2)^^, and 1 for t/(l)y), which leaves a total of 37 complex dimensions 
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[46, 43]. The trick is to construct gauge invariant monomials forming SU{3)c singlets 
and then using them as building blocks to generate SU{3)c x SU (2)^, and subsequently 
the whole SU{3)c x SU{2)l x U{1)y invariant polynomials [46, 43]. However these 
invariant monomials give only the D-flat directions. For F-flat directions, one must 
solve explicitly the constraint equations Fx^ — 0. 

A single flat direction necessarily carries a global U{1) quantum number, which 
corresponds to an invariance of the effective Lagrangian for the order parameter (f) 
under phase rotation e*^0. In the MSSM the global U (1) symmetry is B — L. For 
example, the Li/^-direction (see below) has B — L — —1. 

A flat direction can be represented by a composite gauge invariant operator, X^,, 
formed from the product of k chiral superfields $j making up the flat direction: 
Xm = ^1^2 ■ ■ ■ ^m- The scalar component of the superfield is related to the 
order parameter through Xm = cq 



Lm 




if„ = -= I ° , L ^ 







(164) 



4.2.2 An example of F-and D-flat direction 

The flat directions in the MSSM are tabulated in Table 1. An example of a D-and 
F-flat direction is provided by 

/ 

where is a complex fleld parameterizing the flat direction, or the order parameter, 
or the AD field. All the other fields are set to zero. In terms of the composite gauge 
invariant operators, we would write X^ = LHu {m — 2). 

From Eq. (164) one clearly obtains F^^ = XuQu + /iHd = F£ = XdHdC = for 
all (f). However there exists a non-zero F-component given by F^^ — iiH^- Since /i 
can not be much larger than the electroweak scale Mw ~ TeV, this contribution 
is of the same order as the soft supersymmetry breaking masses, which are going to 
lift the degeneracy. Therefore, following [46], one may nevertheless consider LHu to 
correspond to a F-flat direction. 

The relevant D-terms read 

^W) = ^«^3i^« + L^rsL = _ ^ . (165) 
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B-L 




B-L 


HuHd 







-1 


udd 


-1 


QLd 


-1 


LLe 


-1 


QQud 





QQQL 





QLue 





uude 





QQQQu 


1 


QQuue 


1 


LLddd 


-3 


uuuee 


1 


QLQLdd 


-2 


QQLLdd 


-2 


uudddd, 


-2 


QQQQdLL 


-1 


QLQLQLe 


-1 


QLuQQdd 


-1 


uuuddde 


-1 



Table 1: Renormalizable F and D fiat directions in tiie MSSM 

Therefore the LH^ direction is also D-flat. 

The only other direction involving the Higgs fields and thus soft terms of the order 
of /X is HuHfi- The rest are purely leptonic, such as LLe, or baryonic, such as udd, or 
mixtures of leptons and baryons, such as QLd. These combinations give rise to several 
independent flat directions that can be obtained by permuting the flavor indices. For 
instance, LLe contains the directions L^L'^e^, L'2.L^e\, and LiL^e2. 

Along a flat direction gauge symmetries get broken, with the gauge supermultiplets 
gaining mass by super- Higgs mechanism with irtg = g{4>). Several chiral supermulti- 
plets typically become massive by virtue of Yukawa couplings in the superpotential; for 
example, in the LH^ direction one finds the mass terms VFmass = ^u{<f>)Qu + ^e{<P) Hae. 

Of course, there may simultaneously exist several flat directions. For the purpose 
of AD mechanism it is the lowest dimensional operator which determines the baryonic 
charge of the eventual condensate. In what follows we will therefore mostly consider 
a single flat direction. 
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4.3 Lifting the flat direction 

Vacuum degeneracy along a flat direction can be broken in two ways: by supersym- 
metry breaking, or by higher order non-renormahzable operators appearing in the 
effective low energy theory. Let us first consider the latter option. Supersymmetry 
breaking will then be discussed in more detail in Sects. 4.4 and 4.5. 

4.3.1 Lifting by non-renormalizable operators 

Non-renormalizable superpotential terms in the MSSM can be viewed as effective 
terms that arise after one integrates out fields with very large mass scales appearing 
in a more fundamental (say, string) theory. Here wc do not concern ourselves with the 
possible restrictions on the effective terms due to discrete symmetries present in the 
fundamental theory, but assume that all operators consistent with symmetries may 
arise. Thus in terms of the invariant operators X^, one can have terms of the type 
[45, 46] 

where the dimensionality of the effective scalar operator d — mk, and h is a coupling 
constant which could be complex with \h\ ~ 0{1). Here M is some large mass, 
typically of the order of the Planck mass or the string scale (in the hctcrotic case 
M ~ Mqut)- The lowest value of A; is 1 or 2, depending on whether the flat direction 
is even or odd under i?-parity. 

A second type of term lifting the flat direction would be of the form [45, 46] 

where ■0 is not contained in Xm- The superpotential term Eq. (167) spoils F-flatness 
through F,tp ^ 0. An example is provided by the direction uiU2U3eie2, which is lifted 
by the non-renormalizable term W = [h' / M)uiU2d2ei. This superpotential term gives 
a non-zero contribution = {h'/M)uiU2ei ~ {h'/M)(f)^ along the flat direction. 

Assuming minimal kinetic terms, both types discussed above in Eqs. (166,167) 
yield a generic non-renormalizable potential contribution that can be written as 
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where we have defined the couphng |A|^ = + By virtue of an accidental 

i?-symmetry under which (f) has a charge R = 2/d, the potential Eq. (168) conserves 
the U{1) symmetry carried by the flat direction, in spite of the fact that at the su- 
perpotential level it is violated, see Eqs. (166,167). The symmetry can be violated if 
there are multiple fiat directions, or by higher order operator contributions. However 
it turns out [46] that the B — L violating terms are always subdominant. This is 
of importance for baryogenesis considerations, where the necessary B — L violation 
should therefore arise from other sources. 

The process of finding all the possible non-renormalizable superpotential contri- 
butions lifting a particular fiat direction is similar to finding the D-fiat directions 
discussed in Sect. 4.2.1. All the non-renormalizable operators can be generated from 
SM gauge monomials with i?-parity constraint which allows only even number of odd 
matter parity fields (Q, L, u, d, e) to be present in each superpotential term. At each 
dimension d, the various F = constraints are separately imposed in order to construct 
the basis for monomials. 

As an example, consider fiat directions involving the Higgs fields such as HuHa 
and LHu directions. Even though they are already lifted by the fi term, since fi is 
of the order of supersymmetry breaking scale, for cosmological purposes they can be 
considered fiat, as was discussed in Sect. 4.2.2. At the d — A level the superpotential 
reads 

W, D ^{HuHd? + ^{L,H^){L,H^) . (169) 

Let us assume A, Ajj 7^ 0. Note that Fh^ — constraint implies XH^{HuHd) — 0, 
which acts as a basis for the monomials. An additional constraint can be obtained 
by contracting Fh^ = by Ca/sH^, which forms the polynomial H^Hd = in the 
same monomial basis. Similarly the constraint F//^ = 0, along with the contraction 
yields X'^{LiHy){LjHy) = 0. This implies that LiHu = for all i. Therefore the 
two monomials LH^ and H^Hd can be lifted by 0? = 4 terms in the superpotential 
Eq. (169). 

The other renormalizable fiat directions are LLE, uud, QdL, QQQL, QuQd, uude 
and QuLe, dddLL, uuuee, QuQue, QQQQu, uudQdQd, and {QQQ)4^LLLe. These are 
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lifted primarily by the superpotential terms which involve either Hu or if d is odd, 
or those which contain neither Hu nor Hd. The complete list of superpotential terms 
which lift the flat directions can be found in [43]. It was shown that all the MSSM 
flat directions are lifted by c? = 4, 5, 6, 7, 9 terms in the superpotential. The unique 
flat directions involving {Q, u, e) is lifted by d — 9, {L, d) by d — 7, and (L, d, e) by 
d — 5. The flat directions involving (L, e), {u, d) and (L, d, e) are all lifted by d = 6 
terms in the superpotential, while the rest of the flat directions are lifted already by 
d = A superpotential terms. 

4.3.2 Lifting by soft supersymmetry breaking 

Vacuum degeneracy will also be lifted by supersymmetry breaking, as will be discussed 
in more detail in Sects. 4.4 and 4.5. It is induced by the soft terms, which in the 
simplest case read 

where the supersymmetry breaking mass mo and A are typically of the order of the 
gravitino mass mz/2- An additional soft source for supersymmetry breaking are the 
gaugino masses nig. The A-term in Eq. (170) violates the C/(l) carried by the flat 
direction and thus provides the necessary source for B—L violation in AD baryogenesis. 
In general, the coupling A is complex and has an associated phase 6\. Writing cf) — 
1 01 exp(i^), one obtains a potential proportional to Q.o'$,{Q\~\'n&) in the angular direction. 
This has n discrete minima for the phase of 0, at each of which is broken. 

4.4 Supersymmetry breaking in the MSSM 

In the MSSM there are several proposals for supersymmetry breaking, which we shall 
discuss below. However most of the time it is not important to know the exact mech- 
anism of low energy supersymmetry breaking. This ignorance of the origin of super- 
symmetry breaking can always be hidden by simply writing down explicitly the soft 
breaking terms with arbitrary couplings. 



dMd-^ 



+ h.c. 



(170) 
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4.4.1 Soft supersymmetry breaking Lagrangian 

The most general soft supersymmetry breaking terms in the MSSM Lagrangian can 
be written as (see e.g. [35]) 

Csoft = -\ {M^X'X' + c.c.) - (m2)}0^>, - {^hj4>4j + ^a^^'Vi</',-0fe + cc.) , (171) 

where Mx is the common gaugino mass (m^)^ are 3x3 matrices determining the 
masses for squarks and sleptons, denoted as mQ, m?, m|, m|, m|; 6^ is the mass term 
for the combination Hy^H^, and finally, a^^^ are complex 3x3 matrices in the family 
space which yield the 74-terms a^, a^, ae- There are a total of 105 new entries in the 
MSSM Lagrangian which have no counterpart in the SM. However the arbitrariness 
in the parameters can be partly removed by the experimental constraints on flavor 
changing neutral currents (FCNC) and CP violation [290]. In order to avoid FCNC 
and excessive CP violation, the squark and slepton (mass)^ matrices are often taken 
to be flavor blind, so that the squark and slepton mixing angles can be rotated away. 
Similarly, one may assume that the 0^ couplings are proportional to the Yukawa cou- 
pling matrix, so that a„ = AuqX^, aa — AdoXd, and = ^eo^e- Large CP violating 
effects can be avoided if the soft parameters do not involve new CP phases in addi- 
tion to the SM CKM phases. One can also fix ji parameter and h to be real by an 
appropriate phase rotation of if„ and H^i. 

There are a number of possibilities for the origin of supersymmetry breaking. Fayet- 
Iliopoulos mechanism [291] provides supersymmetry breaking by virtue of a non-zero 
D-term but requires a U{1) symmetry. However, this mechanism does not work in the 
MSSM because some of the squarks and sleptons will get non-zero vevs which may 
break color, electromagnet ism, and/or lepton number without breaking supersymme- 
try. Therefore the contribution from the Fayet-Iliopoulos term should be negligible at 
low scales. 

There are models of supersymmetry breaking by F-terms, known as O'Raifeartaigh 
models [292], where the idea is to pick a set of chiral supermultiplets $i D (0i,'0iFj) 
and a supcrpotential W in such a way that Fj = —5W/ 5(f)* = have no simultaneous 
solution. The model requires a linear gauge singlet superfield in the supcrpotential. 
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Such singlet chiral supermultiplet is not present in the MSSM. The scale of supersym- 
metry breaking has to be set by hand. 

The only mechanism of supersymmetry breaking where the breaking scale is not 
introduced either at the level of superpotential or in the gauge sector is through dy- 
namical supersymmetry breaking [293, 294] . In these models a small supersymmetry 
breaking scale arises by dimensional transmutation. It is customary to treat the su- 
persymmetry breaking sector as a hidden sector which has no direct couplings to the 
visible sector represented by the chiral supermultiplets of the MSSM. The only allowed 
interactions are those which mediate the supersymmetry breaking in the hidden sector 
to the visible sector. 

The main contenders are gravity mediated supersymmetry breaking, which is as- 
sociated with new physics which includes gravity at the string scale or at the Planck 
scale [253, 34], and gauge mediated supersymmetry breaking, which is transmit- 
ted to the visible sector by the ordinary electroweak and QCD gauge interactions 
[295, 265, 296, 297]. There are other variants of supersymmetry breaking based upon 
ideas on gravity and gauge mediation with some extensions, such as dynamical super- 
symmetry breaking (see [298], and references therein), and anomaly mediation (see 
[268, 299], and references therein), which we do not consider here. 

4.4.2 Gravity mediated supersymmetry breaking 

Let us assume that supersymmetry is broken by the vev (F) ^ and is communicated 
to the MSSM by gravity. On dimensional grounds, the soft terms in the visible sector 
should then be of the order [34] 

(172) 

Note that rasojt — > as Mp — > cxo. In order to obtain a phenomenologically acceptable 
soft supersymmetry mass msojt ~ O(IOO) GeV, one therefore requires the scale of 



supersymmetry breaking in the hidden sector to be (F) ~ 10^° — 10^^ GeV. 

Another possibility is that the supersymmetry is broken via gaugino condensate 
(0|A"A^]0) = 5"^A^ 7^ 0, where A is the condensation scale [300, 34]. If the composite 
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field A^A*" belongs to the (F) ~ A^/Mp-term, then again on dimensional grounds one 
would expect the soft supersymmetry mass contribution to be [34] 

^soft ~ ■ (173) 

In this case the nature of supersymmetry breaking is dynamical and the scale is given 
by A ~ 10^=^ GeV. 

The supcrgravity Lagrangian must contain the non-rcnormalizable terms which 
communicate between the hidden and the observable sectors. For the cases where the 
kinetic terms for the chiral and gauge fields are minimal, one obtains the following soft 
terms [34] 

The gauginos get a common mass Mi = M2 = M3 = mi/2, the squark and slepton 
masses are mg = m\ — m| = m\ — m\ — mg, and for the Higgses = m|^^ = rriQ. 
The A-terms arc proportional to the Yukawa couplings while h = Bq/i. 

Some particular models of gravity mediated supersymmetry breaking give more 
detailed estimates of the soft supersymmetry terms. They include: Dilaton dominated 
models [301], which arise in a particular limit of superstring theories, which have 
ml — vp^i2-, and mi/2 = — = ^f^^m■i|2^■, Polonyi models [302], where m,Q — ^3/2, 

= (3 — V^)ms/2, and mi/2 = 0(1713/2); and No-scale models [303], which also arise 
in the low energy limit of superstrings and in which the gravitino mass is undetermined 
at the tree level while the at the string scale mi/2 ^ mo, Aq, 1713/2- 

The predictions for the mass spectrum and other observable can be found renormal- 
ization group (RG) equations; these will be described in connection with the dynamical 
evolution of the AD field. Therefore, a generic fiat direction in gravity mediated super- 
symmetry breaking has two important components: the soft supersymmetry breaking 
terms, and the RG induced logarithmic dependence of the vev. 

4.4.3 Gauge mediated supersymmetry breEiking 

In gauge mediated supersymmetry breaking one employs a heavy messenger sector 
which couples directly to the supersymmetry breaking sector but indirectly to the 
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observable sector via standard model gauge interactions only [295, 304]. As a result 
the soft terms in the MSSM arise through ordinary gauge interactions. There will still 
be gravitational communication, but it is a weak effect. 

The simplest example is a messenger sector with a pair of SU{2) doublet chiral 
fields I, I and a pair of SU{3) triplet fields g, q, which couple to a singlet field z with 
Yukawa couplings A2, A3, respectively. The superpotential is given by 

Wmess = X2ZII + \3Zqq . (175) 

The singlet acquires a non-zero vev and a non-zero F-term (F^). This can be ac- 
complished either substituting z into an O'Raifeartaigh type model [295, 304], or by a 
dynamical mechanism [265, 296, 297]. One may parameterize supcrsymmctry breaking 
in a superpotential Wbreak by {dWi^reak/ dz) = —{F*). As a consequence, the messenger 
fermions acquire masses 

£ = - {X2{z)Mj+Xs{z)i^,^Pg + c.c) , (176) 

while the scalar messenger partners have a scalar potential given by 

V = \\2{z)\' + + \\s{z)\' + \q\') - {\2{F,)ll + \s{F,)qq + c.c.) 

-|-quartic terms , (177) 

where we have used (dW^ess/dz) — 0, and we have replaced z and F^ by their vevs. It 
is easy to read off the eigenvalues of the squared scalar masses and the fermionic and 
bosonic spectrum of the messenger sector; for {I, I), rri^fermions — l'^2(-z) P, and ml^g^i^^g — 
\\2{z)\'^ ± |A2(i^^)|; for (g,g), m^^^i^^^ = \X3{z)\'^, and m^^^;^^^ = \X3{z)\'^ ± |A3(F^)|. 

Supersymmetry breaking is then mediated to the observable fields by one-loop 
corrections, which generate masses for the MSSM gauginos [265] . The q, q messenger 
loop diagrams provide masses to the gluino and the bino, while I, I messenger loop 
diagrams provide masses to the wino and the bino, i.e., Ma^i 2,3 = {aa/4:Tr)A, where 

For squarks and slcptons the leading term comes from two-loop diagrams, e.g. 
oc a^. The 74-terms get negligible contribution at two-loop order compared to 
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the gaugino masses, they come with an extra suppression of a/Ai: compared with the 
gaugino mass, therefore a^, = = Oe = is a good approximation. The Yukawa 
couphngs at the electrowcak scale are generated by evolving the RG equations. 
One can estimate [265] the soft supersymmetry breaking masses to be of order 



If Ms ~ {z) and J (F) are comparable mass scales, then the supersymmetry breaking 



4.5 Supersymmetry breaking in the early Universe 

Non-zero inflationary potential gives rise to supersymmetry breaking, the scale of 
which is given by the Hubble parameter. At early times this breaking is dominant over 
breaking from the hidden sector. After the end of inflation, in most models the inflaton 
oscillates and its finite energy density still dominates and breaks supersymmetry in the 
visible sector. Supersymmetry is broken also by quantum mechanical effects but these 
are negligible compared to the classical supersymmetry breaking from the non-zero 
energy density of the Universe. 

4.5.1 Inflaton-induced terms 

The early Universe supersymmetry breaking can be transmitted to the MSSM flat 
directions either by renormalizable or non-renormalizable interactions [46]. However 
at least for a single flat direction, renormalizable interactions do not lift the MSSM 
fiat directions. In contrast, the effective potential generated by non-renormalizable 
interactions can induce a mass for the ffat direction which is independent of the field 
values as long as they are below the Planck scale. 

At tree level N — 1 SUGRA potential in four dimensions is given by the sum of F 
and D-terms [37] 



aa{F)_ 
An Ms 



(178) 




can take place at about ^/{F) ~ 10^ - 10^ GeV. 




(179) 
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where 

.W 

F'^W' + K'j^, D'' ^ -K'{Tyi(f)j + C . (180) 

where we have added the Fayet-Ihopoulos contribution ^" to the D-term. Here K is 
the Kahler potential, which is a function of the fields 0i, and = dK/d(j)i, and Re fab 
is the inverse of the real part of the gauge kinetic function matrix. 

A particular class of non-renormalizable interaction terms induced by the inflaton 
arise if the Kahler potential has a form [305, 46, 306] 

K = jd'e^ym^), (181) 

where / is the inflaton whose energy density p ^ {J d^OPl) dominates during inflation, 
and (j) is the flat direction. The interaction Eq. (181) will generate an effective mass 
term in the Lagrangian in the global supersymmetric limit, given by 

>C = = 3i/,V0 , (182) 

where Hi is the Hubble parameter during inflation. 

4.5.2 Supergravity corrections 

In addition, there are also inflaton-induced supergravity corrections to the flat direc- 
tion. By inspecting the supergravity potential, one flnds the following terms 

\W{I)\^\ 



Ml, 



and (K,K''Dr'''^P'U^.c). (183) 

Above Di = d/dl + KiW/Mp. All these terms provide a general contribution to the 
flat direction potential which is of the form[46] 

1^(0) = H'M^f , (184) 

where / is some function. Note that this contribution exists also when the flat direction 
is lifted by non-renormalizable superpotential terms. 
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For a minimal choice of flat direction Kahler potential K{(f)\ 0) = (f)^(f), during 
inflation the effective mass for the flat direction is found to be [46] 

Here it has been assumed that the main contribution to the inflaton potential comes 
from the F-term. If there were D-term contributions Vd{I) to the inflationary poten- 
tial, then a correction of order Vf{I) / {Vf{I) + Vd{I)) must be taken into account. 
In purely D-term inflation there is no Hubble induced mass correction to the flat di- 
rection during inflation because Fj = 0. However, when D-term inflation ends, the 
energy density stored in the D-term is converted to an F-term and to kinetic energy 
of the inflaton. Thus again a mass term = ±0{1)H'^ appears naturally, however 
the overall sign is undetermined [228]. 

There are additional inflationary contributions to the potential if the flat direction 
is lifted by the non-renormalizable operators discussed earlier in this Section. These 
new terms come explicitly from the superpotential part of the flat direction 

^W^K^'^K^^^ + h.c^ , (w^K^'DjW*{I)+ h.c.) , 

The flrst one comes from the cross term between the derivative of the flat direction 
superpotential and the inflaton superpotential, the second is due to the Kahler poten- 
tial coupling between the flat direction and the inflaton, and the third term is a cross 
term between the two superpotentials. All these terms give a generalized contribution 
equivalent to an A-term of the MSSM: 

V{^) = HMlf , (187) 

where d is the dimensionality at which the flat direction is lifted. The induced A-ierm. 
has an important role to play during the evolution of the flat direction. A possible A- 
term can also be generated from the expansion of the Kahler potential for fleld values 
I,(f> < Mp, which is of the form [46] 



/ d'9I<P^ + h.c. ^ + h-c. , 

Mp y Mp 
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l-Jd^eiW, + b.c.~-^m + h.c. (188) 

Note that the A-terms arise only from terms with a hncar couphng of the inflaton 
superfield to a gauge invariant-operator 0,. If / were a composite field rather than a 
singlet, then such a term will not arise and an A-term will not be generated. Also, in 
the case of D-term inflation, the inflaton cannot induce an A-term because Fj — 0. 
More generally, if there is a symmetry preventing a linear coupling of the inflaton, then 
order H A-terms can be eliminated also in F-term inflation. As long as the thermal 
bath of the inflaton decay products dominates over the low energy supersymmetry 
breaking scale, we should have Hubble induced corrections to rriQ, 1111/2, A. 

If there is a non-minimal dependence of the gauge superfield kinetic terms on the 
inflaton field, a Hubble-induced gaugino mass can also be produced. Generally the 
gauge superfield kinetic terms must depend on the field(s) of the hidden sector in order 
to obtain gaugino masses of roughly the same order as (or larger than) scalar masses, 
as required by phenomenology. Having m.1/2 ~ H thus appears to be quite natural 
unless an it!-symmetry forbids terms which are linear in the inflaton superfield [46] . 

Since the //-term does not break supersymmetry, there is a priori no reason to 
assume that a //-term of order H will be created. (For a discussion, see [279]). In what 
follows we will treat as a free parameter. 

So far we have not discussed the sign of the Hubble induced mass correction. In 
fact with a general Kahler term either sign is possible. Depending on the sign, the 
dynamical behavior of the AD field is completely different and therefore the predictions 
depend crucially upon the sign. There are however certain cases where the Hubble- 
induced terms might not occur at all. An i?-symmetry [45, 46] or a special choice 
of the Kahler potential could forbids the AD field getting the Hubble-induced mass 
correction [226, 305]. 
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4.6 The potential for flat direction 
4.6.1 F-term inflation 

Let us collect all the terms which contribute to the flat direction potential, which in 
the case of F-term inflation can be written as [45, 46] 

(189) 

The first and the third terms are the Hubble-induced and low-energy soft mass terms, 
respectively, while the second and the fourth terms are the Hubble-induced and low- 
energy A terms. The last term is the contribution from the non-renormalizable super- 
potential. The coefficients |C/|, a, ~ 0{1), and the coupling A ~ l/((i — 1)!. Note 
that low-energy term is dimensionful. 

Note here the importance of the relative sign of the coefficient Cj. At large field 
values the first term dictates the dynamics of the AD field. If C/ < , the absolute 
minimum of the potential is (f) — and during inflation the AD fleld will settle down 
to the bottom of the potential roughly in one Hubble time. In such case the AD field 
will not have any interesting classical dynamics. Its presence would nevertheless be 
felt because of quantum fluctuations. These would be chi-squared in nature since then 
the classical energy density of the AD fleld would be due to its own fluctuations. 

If C/ <^ 1, the AD fleld takes some time to reach the bottom of the potential, and 
if it has a non-zero amplitude after the end of inflation, its dynamics is non-trivial. 

The most interesting scenario occurs when Cj > 0. In this case the absolute value 
of the AD fleld settles during inflation to the minimum given by 

' ^ -HjM'^-'^] . (190) 



Here we have ignored the potential term oc a; if C/ > 0, the a-term will not change 
the vev qualitatively. On the other hand, even for C/ < the potential Eq. (189) will 
have a minimum with a non vanishing vev if |ap > 4{d ~ l)Cj. However the origin 
will also be a minimum in this case. The dynamics then depends on which minimum 
the AD fleld will choose during inflation. 
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The a-term in Eq. (189) violates the global U{1) symmetry carried by 0. If \a\ 
is C(l), the phase 9 of (0) is related to the phase of a through nO + 6a = vr; oth- 
erwise 6 will take some random value, which will generally be of 0{1). This is the 
initial CP- violation which is required for baryogenesis/leptogenesis. In practice, the 
superpotential term lifting the flat direction is also the B and CP violating opera- 
tor responsible for AD baryogenesis, inducing a baryon asymmetry in the coherently 
oscillating 4> condensate. 

4.6.2 D-term inflation 

In D-term inflation one does not get the Hubble induced mass correction to the flat 
direction so that C/ = 0. Also the Hubble induced a-term is absent. However the 
Hubble induced mass correction eventually dominates once D-tcrm induced inflation 
comes to an end. The potential for a generic flat direction during D-term inflation is 
given by 

V{<l>) = mlld^f + [a,X,^ + h.c.) + \Xf\^ , (191) 
and after the end of inflation the flat direction potential is given by [228] 

y(0) = {ml - CH^) + (^.A,^ + h.c.) + \X\^\^ , (192) 

where C ~ For C positive, the flat direction settles down to one of its minima 

given by Eq. (190) provided > y^m^M/A, otherwise 

Note that in this case that the 74-term is also responsible for B and/or L, and CP 
violation. Another generic point to remember is that in i?-parity conserving models 
the B and/or L violating operators must have even dimensions, so that d — A yields 
the minimal operator for AD baryogenesis. 
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5 Dynamics of flat directions 

After the end of inflation (0) continues first to track the instantaneous local minimum 
of the scalar potential, obtained by replacing Hj with H{t) in Eq.(190) or by follow- 
ing Eq. (193) in the D-term inflation case. Once H ~ itlq ~ ^^1^3/2, the low-energy 
soft terms take over. Then m| becomes positive and (0) starts to move in a non- 
adiabatic way (the phase of (0) differs from the phase of A-term during inflation). 
As a result (0) begins a spiral motion in a complex plane, which charges up the fiat 
direction condensate, and eventually leads to generation of a net baryon and/or lepton 
asymmetry [46]. 

For baryogenesis purposes it is essential that the AD condensate obtains a non-zero 
vev during the inflationary epoch. In Sect. 4, we pointed out that a non-zero vev of 
the flat direction condensate is acquired only when the negative (mass)^ contribution 
dominates the potential. The MSSM flat directions which arc made up of squarks 
and sleptons have Yukawa and gauge interactions. The couplings render the evolution 
of a particular flat direction non-trivially, especially when the flat direction has a 
time varying mass due to the Hubble expansion [305, 46, 307, 308]. Moreover, if 
thermalization is not instantaneous, thermal effects from reheating can be substantial 
and might trigger the motion of the flat direction at an earlier time, there by changing 
the evolution of the flat direction condensate in a significant way [309, 310]. 

5.1 Running of the couplings 
5.1.1 Running of gravitational coupling 

Any flat direction has two kinds of interactions: renormalizable gauge or Yukawa in- 
teractions, and a non-trivial coupling to the curvature. Both types of interactions 
contribute to the logarithmic running of (mass)^ of the flat direction condensate. The 
coupling to the curvature is generic because in principle any scalar fleld in an expand- 
ing background receives a contribution from the curvature by virtue of the Lagrangian 
term ^R(f)^, where is a coupling constant. Note that R oc in an expanding 

background. Any scalar fleld always gets an additional positive Hubble induced mass 
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correction, provided ^ is positive. The fundamental theory might have a conformal 
invariance, in which case the couphng strength ^ = 1/6 [311], but it is known that 
conformal invariance is not protected by any symmetry, and that quantum corrections 
always break conformal invariance. Especially for the flat direction condensate, spon- 
taneously broken supersymmetry induces soft supersymmetry breaking terms which 
break conformal invariance, and the value of ^ remains undetermined. 

It is of course possible to simply set ^ = 0. If initially ^ = at some high scale, 
renormalization effects due to scalar field self-interaction will nevertheless generate a 
non-zero ^ at lower scales. In an expanding Universe the value of ^ also changes under 
the influence of a varying curvature (see [312], and references therein). In the simplest 
case of a single scalar field with a quartic self-interaction strength A leads (at one-loop 
level) to a logarithmically running ^ [313] 




(194) 



It is obvious that ^ = 1/6 is a fixed point of the RG equation. If the theory has 
fermions and gauge fields, then obviously the coefficient in front of the logarithmic 
term in Eq. (194) will be modified [312]. 

As we have seen in Sect. 4.2, when supersymmetry is promoted to a local theory, 
a supergravity correction is induced to the fiat direction which is proportional to 
the curvature, and supergravity theories also allow for ^R4>^ (e.g. super conformal 
supergravity [270]). 

In the context of MSSM fiat directions we have implicitly assumed ^ = 0. This is 
justified from the very definition of F- and D-fiat directions. The only leading order 
self coupling term in the fiat direction potential is the Hubble induced A-term in 
Eq. (189). The overall self coupling constant is relatively large when the flat direction 
is lifted at 0? = 4, i.e. the suppression is proportional to 0{1){H/Mp), where we have 
replaced M by Mp in Eq. (189). In any inflation model the ratio Hj/Mp <^ 1, which 
in conjunction with Eq. (194), suggests that the effect of running on ^ is minimal. 
For a running ^ the curvature term in Eq. (194) dominates over the mass term. This 
might not be the case with the flat direction condensate because the condensate also 
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receives a field dependent mass while it is evolving. As long as the vev dependent 
mass is larger than the curvature induced mass, the running of any parameter in the 
theory will be dictated mainly by the renormalizable quantum efi^ects. We therefore 
conclude that the running of ^ can be neglected. 

When the field dependent mass of the fiat direction field becomes of order ~ 
0{H^), it might be prudent to start worrying about the curvature induced term, 
especially during infiation. A simple inspection of Eq. (194) suggests that ^ is always 
of order ln(0(l)), with virtually no alteration in ^eff- From now onwards we fix the 
non-minimal coupling to be ^ = 0. 

5.1.2 Renormalization group equations in the MSSM 

Let us consider the running of the flat direction (mass)^ below Mqut by assuming 
that it is the scale where supersymmetry breaking is transmitted to the visible sector, 
in order to avoid uncertainties about physics between Mqut and Mp. The running 
of low-energy soft breaking masses has been studied in great detail in the context of 
MSSM phenomenology [314], in particular in connection with radiative electroweak 
symmetry breaking [315]. 

Let us recall some of the salient features of the MSSM one-loop RG equations. The 
ones relevant to fiat directions involve the Higgs doublet which couples to the top 
quark, the right-handed stop us, the left-handed doublet of third generation squarks 
Qs and the A— parameter At associated with the top Yukawa interaction. The RG 
equations read [34] 

^ 2 '^^t / 2 I 2 I 2 , \ A \2\ ^ 2\ |2 , ^ 21 \2\ 

dq'^H. = ^ K« + "^Qa + + 1^*1 J - ^ [l9i\mi\ + -9M j , 

2 ^t/2 I 2 I 2 , \ A \2\ ^ f ^ 2\ |2|32| \2 , ^ 2\ \2\ 

—m~ = — [TnTr +m~ +m~ + \At\ -[ — Qiimil H — Qolmol H — q-timsi , 

d , 3/i2 1 /13 2 3 2 4 2 \ 

Here q denotes the logarithmic scale; this could be an external energy or momentum 
scale, but in the case at hand the relevant scale is set by the vev(s) of the fields 
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themselves, ht is the top Yukawa couphng, while Qi and are respectively the gauge 
couphngs and soft breaking gaugino masses of f/(l)y x SU{2) x SU{?>) . If ht is 
the only large Yukawa coupling (i.e. as long as tan/9 is not very large), the beta 
functions for (mass)^ of squarks of the first and second generations and sleptons only 
receive significant contributions from gauge/gaugino loops. A review of these effects 
can be found in [314]. Here we only mention the main results for the case of universal 
boundary conditions, where at Mqut all the scalar masses are ttIq and the gauginos 
have a common soft breaking mass mi/2- For a low value of tan/3 = 1.65 ^, 

"^k - -^"^0 - 2^2/2 (196) 

at the weak scale, while and m~ remain positive. The soft breaking (mass)^ 
of the first and second generations of squarks is ^ rn^ + (5 — 7)m^^2! while for the 
right-handed and left-handed sleptons one gets ~ ttIq -|- 0.1m^y2 ~ ttIq -|- 0.5m^y2) 
respectively. The important point is that the sum m|^^ -|- m|, which describes the 
mass along the H^L flat direction, is driven to negative values at the weak scale only 
for mi/2 > ruQ. This is intuitively understandable, since Eqs.(195) have a fixed point 
solution [317] m^Hu + "^ug + = = when mi/2 = 0. 

5.2 Hubble induced radiative corrections 

Here we describe radiative corrections in a cosmological set-up relevant for the AD 
mechanism [308]. When the Hubble induced supersymmetry breaking is dominant, 
i.e. for H > O(TeV), the evolution of the soft terms is different from the vacuum 
RG equations given in Eq. (195). For the low-energy supersymmetry breaking case, 
constraints from the weak scale (e.g. realization of electroweak symmetry breaking, 
and experimental limits on the sparticle masses) give information about the soft break- 
ing parameters ml and mi/2- Together with fine tuning arguments, these constraints 
^This value corresponds to the case of maximal top Yukawa coupling, so called fixed point scenario 
[316, 317]. Such a low value of tan /3 is excluded by Higgs searches at LEP [318], unless one allows stop 
masses well above 1 TeV. We nevertheless include this scenario in our discussion since it represents 
an extreme case. 
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imply that ttIq > and that mo, mi/2 are C(TeV). In the Hubble induced supersym- 
metry breaking case ml and mi/2 are determined by the scale of inflation (and the 
form of the Kahler potential). At low scales the Hubble induced terms are completely 
negligible because at temperature T ~ M^, H ~ 0{1) eV; at present the Hubble 
parameter is tiny, Hq ~ O{10~^^) eV. 

There exists an even more fundamental difference between the Hubble induced 
and ordinary radiative corrections. In Minkowski space the loop contributions to 
beta functions freeze at a scale of the order of the mass of the particles in the loop. 
In an expanding Universe the horizon radius oc defines an additional natural 
infrared cut-off for the theory. The masses of particles coupled to the flat direction 
receive contributions from two sources. There is a supersymmetry preserving part 
proportional to the vev (0), and the Hubble induced supersymmetry breaking part. 
The loop contributions to beta functions should thus be frozen at a scale given by 
the largest of and H (recall that ht and gauge couphngs are close to one). In 
particular, if the squared mass of the flat direction condensate is positive at very large 
scales but turns negative at some intermediate scale Qc, the origin of the flat direction 
potential will cease to be a minimum, provided the Hubble parameter is less than 
Qc- On the other hand, if < at the GUT scale, its running should already be 
terminated at the scale \{(f))\ determined by Eq.(190).'^ 

In the following two subsections we discuss separately the cases of positive and 
negative GUT-scale (mass)^ for the flat direction condensate. 

5.2.1 The case with d pa -1 

In this case all scalar fields roll towards the origin very rapidly and settle there during 

inflation, provided radiative corrections to their masses are negligible. A typical flat 
^Here we note that the Hubble cut-off usually plays no role in loop corrections to the inflaton 
potential. In most inflation models the masses of the flelds which may run in the loop arc larger than 
the Hubble expansion rate due to the presence of a flnite coupling to the inflaton. This will happen 
if the (time varying) inflaton vev is large and the couplings are not very small. In those cases, which 
are somewhat similar to our case with C/ > 0, one can trust the usual loop calculation evaluated in 
a flat space time background [218]. 
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At/H 


mi/2/H 


Qc{ht = 2) 


Qc{ht = 0.5) 


+1/3 (-1/3) 


1/3 


X 


X 


+1/3 (-1/3) 


1 


10^ - 10^ 


10^ 


+1/3 (-1/3) 


3 


10^1 


10^ - 10^ 


+1 (-1) 


1/3 


X 


X 


+1 (-1) 


1 


10^ - 10^ 


10^ (x) 


+1 (-1) 


3 


10^1 


108 (106) 


+3 (-3) 


1/3 


X 


10^ 


+3 (-3) 


1 


10^4 (10^) 


109 (103) 


+3 (-3) 


3 


10^5 (IQii) 


10^0 (10^) 



Table 2: The scale Qc (in GeV) where the squared mass of the if„L flat direction 
changes sign, shown for Cj = —1 and several values for the ratios At/H and mi/2/H 
as well as the top Yukawa coupling ht, all taken at scale Mqut = 2 • 10^^ GeV, from 
[308] 

direction condensate ^ is a linear combination = Y,^i o-iVi of the MSSM scalars </7j, 
implying that = Yh=i In [308], it was noticed that with small values for 

the 11 parameter, the running of crucially depends on mi/2- 

Let us consider sample cases with gaugino masses mi/2 = {H\ 3H; H/3), the 
A-term^ ^^(Mgut) = {±H; ± 3H;±H/3 ), top Yukawa /^^(Mgut) = (2, 0.5) and 
couphngs 5'i(Mgut) = 5'2(-^gut) = 5'3(-^gut) = 0.71, and follow the running of scalar 
soft masses from Mqut down to 10^ GeV, where low-energy supersymmetry breaking 
becomes dominant. The main result is that only the LH^ flat direction can acquire 
a negative (mass)^ at low scales. In this case = [mfj^ + m\ + /U^)/2, where the 
last term is from the Hubble induced term. The results are summarized in Table 2, 
where it has been assumed that //(Mqut) ^ H/A so that the the //-term contribution 
to m| is neghgible. In general changes sign at a higher scale for /it (Mqut) = 2. 

^The RG equations (195) for At show that the relative sign between At and mi/2 matters, since it 
affects the running of \At\, and subsequently, scalar soft masses. Without loss of generality we take 
the common gaugino mass mi/2 to be positive. 
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This is expected since a large Yukawa coupling naturally maximizes the running of 
rriff^. Furthermore, the difference between At/nii/2 < and At/mi/2 > becomes 
more apparent as \At/'mi/2\ increases and ht decreases. The quasi fixed-point value of 
At/mi/2 is positive [316, 317]. Positive input values of A^ will thus lead to positive A^ 
at all scales, but a negative ^((Mgut) implies that At ^ for some range of scales, 
which diminishes its effect in the RG equations, see Eq. (195). The sign of ^((Mgut) 
is more important for smaller ht, since then At/mi/2 will evolve less rapidly. 

It was noticed in [308] that the squared mass of the H^L fiat direction does not 
change sign when mi/2 — H/3, except for^ At — ±3H and ht — 0.5. This can be 
explained by the fact that for small mi/2 and small or moderate \At\ we are generally 
close to the fixed point solution 

mi ~ --H^; ~ 0; ~ -H^. (197) 

Nevertheless, even for mi/2 <^ H the squared mass of the LHu flat direction as well as 
are < 0.2H^ above 1 TeV, because of the fixed point behavior. This implies that 
the LHu flat direction can still be viable for baryogenesis, as pointed out by McDonald 
[324]. Flat directions built out of will be marginal at best, since the decrease in 
will be counteracted by other contributions to m^; e.g. for the u^did2 flat direction 
we find > 2i?^/3 at all scales. 

The AD mechanism for baryogenesis should always work if Qc> Hj, since in that 
case the global minimum of the potential during inflation is located at |(0)| 7^ 0. Note 
that in this case the vev \{(f>)\ is usually determined by Qc rather than by Eq. (190). 

For scales close to Qc the mass term in the scalar potential Eq. (189) can be 
written as /S^iif^l^p log(]0|/(5c), where the coefficient /S^ can be obtained from the RG 
equations. If > 0, which is true for the Hy,L flat direction for Cj < 0, this term will 
reach a minimum at logd^l/gj = -1- HQc < (HiM^y^y/'^-^ the non-renormalizable 
contributions to the scalar potential are negligible for |0| ^ Qc, so that the minimum 
of the quadratic term essentially coincides with the minimum of the complete potential 

given by Eq. (189). In models of high scale inflation (e.g. chaotic inflation models), 
^For this choice of parameters, At runs initially very slowly. It will therefore remain large for some 
time and helps to decrease quickly towards lower scales. 
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the Hubble constant during inflation Hi can be as large as 10 GeV. This implies that 

for the H^L flat direction can only become negative during inflation if ^ 
which includes the "no-scale" scenario studied in [305]. 

The region of the parameter space safely allowing AD leptogenesis is much larger 
in models of intermediate and low scale inflation (e.g. some new inflation models [14]) 
where Hi is substantially smaller. In such models one can easily have Hi < Qc at 
least for the HuL fiat direction, unless ^ Lf^ or fi"^ > m\i2 [308]. 

If Qc < Hi, the condensate settles at the origin during inflation and its post- 
inflationary dynamics will depend on the process of thermalization. If the inflaton 
decay products thermalize very slowly, is only subject to zero-temperature radiative 
corrections and {(p) can move away from the origin once H < Q^, a necessary condition 
for this scenario is that inflatons do not directly decay into fields that arc charged under 
SU{3) X SU{2) X U{1)y- If Qc > 1 TeV, will readily settle at the new minimum 
and AD leptogenesis can work. 

The situation will be completely different if inflatons directly decay into some 
matter fields. In such a case the plasma of inflaton decay products has a temperature 
T ~ (TdH Mpi^^^^^y^'^ [22] (Td is the inflaton decay rate). Fields which contribute to 
the running of are in thermal equilibrium (recall that the flat direction fleld is stuclc 
at (f) — 0) and their back reaction results in thermal corrections of order to m^. 
For generic models of inflation T > H, implying that thermal effects exceed radiative 
corrections. Therefore (0) remains at the origin at all times and AD leptogenesis will 
not work. 

5.2.2 The case with Ci ^ +1 

In this case all flat directions are viable for baryogenesis purposes provided the running 
of is negligible. Radiative corrections may change the sign (in this case to positive) 
at small vev(s), resulting in the entrapment of 4> at the origin. 

A quantitative study of the sample cases discussed above can be summarized as 
follows [308]. The squared mass of the LHy, flat direction is always negative at small 
scales, unless //^ > H'^/2. For mi/2 = SH, changes sign twice; it is positive for 
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scales Q between roughly 10^^ and 10^ GeV, the precise values depending on ht and 
At. Slepton masses only receive positive contributions from electroweak gauge/gaugino 
loops. As a result, the squared mass of the LLe fiat direction remains negative down to 
1 TeV, unless mi/2 > 2i?; for mi/2 > 3i?, Qc > 10^ GeV even for this flat direction. 
The squared masses of all squarks (except U3) change sign at Qc > 1 TeV unless 
mi/2 < H/3; we find Qc ~ 10^° (10^^) GeV for mi/2/H = 1 (3). This is due to the 
large positive contribution oc mg to the squared squark masses at scales below Mgut- 
The corresponding values for the u^didj and LQd fiat directions are usually somewhat 
smaller, due to the Yukawa terms in the function and the slower running of the 
slepton masses, respectively; however, the values of Qc hsted in Table 2 are still a fair 
approximation. 

The positive contribution to the scalar potential from the non-rcnormalizablc su- 
pcrpotcntial term now dominates —H"^ (see Eq. (190)). If Qc > (-f/'/MQu^^)^/'^"^, 
is positive for all vev(s) and hence the flat direction will settle at the origin during 
inflation and remain there. In such a case the flat direction is not suitable for AD 
baryogenesis. This can easily happen for flat directions involving squarks in mod- 
els with low scale inflation, but is not likely for high scale inflation models (unless 
'mi/2 > 3H). For Hj <Qc< {HjM^Jj^y/''-''^ , feasible for some flat directions in both 
intermediate/high scale and low scale models. During inflation the potential has two 
minima, at (0) = and at |(0)| ~ (iJ/MQux)^^'^"^- Depending on the initial condi- 
tions, (p can roll down towards either of them and settle there but only the latter one 
will be useful for AD baryogenesis. 

If Qc < Hi, the flat direction condensate will settle at the value determined by 
Eq. (190) (the only minimum during inflation) and remain there. The appearance of 
another minimum at the origin after inflation, which is possible once H < Qc, does 
not change the situation since these minima are separated by a barrier. In this case 
radiative corrections will not change the picture qualitatively; however, they will still 
modify the quantitative analysis, since C/ in Eq.(190) will become scale-dependent. 

In brief, the main conclusion is that among the flat directions LHy^ is the only 
robust one in the sense that it gives rise to AD leptogenesis independently of the sign 
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of Cj. 



5.2.3 Running of the flat direction field in no-scale supergravity 

So far we have dealt with a minimal choice of the Kahler potential. An alternative is 
non-flat Kahler potential; an example of this is provided by e.g. no-scale models, for 
which K ~ la{z + z* + where z belongs to supersymmetry breaking sector, and (pi 
belongs to the matter sector, and both are measured in terms of reduced Planck mass 
(for a review, see [303]). In no-scale models there exists an enhanced symmetry known 
as the Heisenberg symmetry [325], which is defined on the chiral fields as Sz — e*0*, 
Scff — e^, and = 0, where are the hidden sector fields, such that the combinations 
T] — z + z* — and yi — are invariant. For a especial choice 

K = f{ri) + HW{<P)/Mlf + g{y) , (198) 



The A^ = 1 supergravity potential reads [305, 307] 



\w\' _ 1 m\\ . ,\w\'^ 

Ml, "^1 f,2 n/r2 ^ 9a\g )b9 ,^2 



(199) 



Note that there is no cross term in the potential such as As a consequence 

any tree level flat direction remains flat even during inflation [305] (in fact it is the 
Heisenberg symmetry which protects the flat directions from obtaining Hubble in- 
duced masses [325]). The symmetry is broken by gauge interactions or by coupling in 
the renormalizablc part of the Kahler potential. Then the mass of the flat direction 
condensate arises from the running of the gauge couplings. 

For f{r]) — —3 In 7?, the one-loop corrected supergravity induced mass term has 
been calculated in [326, 327, 305], which gives an effective mass for the flat direction 
field in the presence of finite energy density stored in the infiaton sector. The typical 
mass of the fiat direction has been computed and comes out to be m| ~ 10~^H^ during 
inflation. The only constraint is that flat direction must not involve stops [305]. 

5.3 Post-inflationary running of the flat direction 

Now we focus on the running of the flat direction after inflation. Here we must take into 
account the low energy supersymmetry breaking effects. In particular, the running of 
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the condensate mass will depend on how supersymmetry is transmitted to the visible 
sector. 



5.3.1 Gravity mediated supersymmetry breaking 

For gravity-mediated supersymmetry breaking the scalar potential along a flat direc- 
tion has been evaluated as [48, 49] 

l + «og(Ll))|*r + -LI_+(^ + A.e.) , (200) 

where is the conventional gravity-mediated soft supersymmetry breaking scalar 
mass term (m^ ~ 100 GeV), X is a parameter which depends on the flat direction, and 
the logarithmic contribution parameterizes the running of the flat direction potential 
with M = (Mp~^m3/2/|A|)*^^/'^~^). In the gravity mediated case \Ax\ < dm^/2-i for 
d = 4,6. 

K can be computed from the RG equations, which to one loop have the form 



dt 



Y: a,,ml + Yhl{Y: h,,m] + A^) , (201) 



where Ojg and &y are constants, mg is the gaugino mass, ha the Yukawa coupling, A 
is the 74-term, and t — In Mx/q- The full RG equations have been listed in [34]. The 
potential along the flat direction is then characterized by the amount of stop mixture 
(where appropriate), the values of gluino mass and A, and in the special case of the 
d — A HuL -direction, on the H^Hd -mixing mass parameter ^h- 

The mass of the AD scalar (f) is the sum of the masses of the squark and slepton 
flelds 0j constituting the flat direction, m| = Z^aPf?7T.f , where Pi is the projection of 
(p along 0j, and = 1. The parameter K is then given simply by 

1 dmi 



''S 



dt 



(202) 

t=\ogq ^ ' 



To compute one has to choose the scale q. The appropriate scale is given by the 
value of the AD condensate amplitude when it flrst begins to oscillate at if or 



Q = \<f>o\ 
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Figure 1: Contours of K for two d=4 flat directions in the {A,^ = mg/m{t = 0))- 
plane: (a) X = (b) K = -0.01; (c) K = -0.05; (d) K = -0.1. The directions 
are (i) Q3Q3QL; (ii) QQQL, no stop; (iii) Usude; (iv) uude with equal weight for all 
■u-squarks, from [329]. 

The RG running of the flat directions in the case of gravity mediated supersymmetry 
breaking was studied in [329, 330], where unification at t = was assumed and all the 
other Yukawa couplings except the top Yukawa were neglected. 

The contours of K for the d — 4 uude and QQQL directions are shown in Fig. (1) 
in the (A, ^)-plane, where ^ = mg/m{t — 0) (for tan/3 = 1 and A = 1). These are 
representative of all the other directions, too, except for Hy^L. For ^ ~ typical 
value for K is found to be about —0.05. Similar contours can be obtained for the d = 6 
{uddY and {QLdY directions, see [329]. For all the squark directions with no stop, as 
long as hb and hy, can be neglected, K is always negative, and the contours of equal 
K do not depend on A. In the presence of stop mixing X < is no longer automatic 
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even in the purely squark directions. The more there is stop, the larger value of ^ is 
required for K < 0. Even for pure stop directions, positive K is typically obtained 
only for relatively light gaugino masses with ^ < 0.5. 

In contrast to the squark directions, K was found [329] to be always positive in 
the if„L-direction. This is due to the fact HuL does not involve strong interactions 
which in other directions are mainly responsible for the decrease of the running scalar 
masses. Very roughly, instability is found when m~ > mj, although the exact condition 
should be checked case by case. In general, the sign of K could be deduced from the 
observation of SUSY parameters such as tan the gluino mass and the supersymmetry 
breaking parameter m<^ [329]. 

5.3.2 Gauge mediated supersymmetry breeiking 

A similar analysis can be made for the gauge mediated case, where supersymmetry 
breaking is transmitted to the observable sector below some relatively low messenger 
sector scale A5, above which the potential is completely flat (see Sect. 4.4.3.). In the 
gravity mediated scenario the soft masses stay intact, modulo RG running, up to the 
Planck scale; in gauge mediation the masses simply disappear above A^. For a large 
condensate vev, one can integrate out the gauge and chiral fields coupled to the flat 
direction in order to obtain an effective low energy theory. In such a case, as was first 
pointed out by Kusenko and Shaposhnikov [151], the potential along the flat direction 
obtains a logarithmic correction of the form [151, 331] 

= mllog 1 + + ' + ^^^^ + h.c. . (204) 

where ~ 1 — 100 TeV. Because of the differences in the potential, the dynamical 
evolution of the condensate field will be markedly different from the gravity mediated 
case. Because the messenger sector is not constrained by experiments, one cannot 
provide a detailed description of the mass parameters. Here one should note that in 
order to have an AD condensate, the A-term is actually constrained. In the gauge 
mediated case I^a] < (10~^ — 10^'^)m<^, for = 4, 6. 

Eqs. (200) and (204) are two book-keeping equations which are useful for the rest 
of this review. 
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5.4 Density perturbations from the flat direction condensate 

The role of MSSM flat directions is not just limited to generating the lepton and/or 
baryon asymmetry in the Universe, but they also play an interesting role in the dy- 
namics of density perturbations. 

5.4.1 Energetics of flat direction and the inflaton field 

Once the flatness of the flat direction potential is lifted by non-renormalizable terms, 
for large fleld values the condensate energy density can dominate over the inflaton 
potential. This could be disastrous: either inflationary expansion would come to a 
halt, or the flat direction condensate fluctuations might ruin the successful predictions 
for the angular power spectrum [52, 332, 53, 54]. 

In [53], the generation of adiabatic density perturbations was studied for both D- 
and F-term inflation models. Note that in the former case there is no Hubble induced 
mass correction to the flat direction condensate. The scalar potential for F-and D-flat 
direction of dimension d is given by (see Eq. (191) in Sect. 4.6.2.) 

where only the dominant term from Eq. (191) corresponding to superpotential term 
of the form W = X^"'/nM'^~^ has been kept. Throughout this discussion i?-parity is 
conserved and therefore we deal only with even dimensions d = 4, 6, 8. 

For illustrative purposes let us assume that the D-term inflationary potential is 
given by (see Eq. (146) in Sect. 3.5.3.) 

where S is the inflaton component; Q is here the renormalization scale. For a large 
initial vev for (f),S ~ 0{Mp), the dynamics is flrst dominated by V{(f)). For a suf- 
flciently large vev of the effective condensate mass squared ^"(0), becomes larger 
than H'^. This occurs ii (j) > (j)H, where [53] 

2W=^ fq\— 2 #4 , ^ 

<1>H = ^ T i—M^--. (207) 



(6(2d-2)(2d-3))5P=:^ 
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If initially 0i > (pH, then (f) will first oscillate in its potential with a decreasing ampli- 
tude: (p{t) oc a~^/'^{t) [230]. This period ends before the onset of infiaton domination. 
The system then enters a regime where both and S are slowly rolling. 
The slow rolling dynamics of the scalar fields is given by the solution of 

where \l'a = S, (p. By taking the ratio of the equations for (p and S, one obtains 

dS ' ^ ' 

which has a general solution of the form 

= 0. [l + (5^ - S')] ; a. = ^,_,^2(.-3)^,^, ' (210) 

where (pi and Si are the initial values at the onset of inflation. There are two features 
about this solution. First, since 5'^ is large compared with the value of 5" at A?^ = 50 
e-foldings before inflation, we see that for sufficiently large 0^ the value of at late 
times is fixed by Si, 

1 \ 2d-4 1 

(211) 



This is true if 0j > 0*, otherwise, (j) simply remains at 0j. Second, we can relate Si 
to the total number of e-foldings during the V{S) dominated period of inflation. In 
general, for sufficiently large 0^, we could have an initial period of V{(j)) dominated 
inflation. During this period S does not signiflcantly change from Si. The potential is 
dominated by T^(0) once > 05, where [53] 

^,= ^(ifY^ . (212) 

\d-i \ ^ J 

(f)s is generally less than 0^ (see Eq. (207)), therefore will be slow rolling during 
V{S) domination. 

From Eq. (210), we flnd that the condition for 5* to change signiflcantly from Si at 
a given value of is given by 

S. < -L (i)"' . (213) 
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The condition for S to change significantly during dominated infiation is given 
by Eq. (213) with = 0^; one finds 

d-2 d 2 

Si<S,,^ ^ ^ ^ . (214) 

Since Si c is small compared with Mp, whereas the value of S required to generate 50 e- 
foldings of inflation, 5*50 = gy/50Mp / {2n) is close to Mp, it follows that Si (> 5*50) will 
generally be larger than Si c, and so the inflaton will remain at Si until the Universe 
becomes inflaton dominated. 

In this case the total number of e-foldings during inflaton domination is given by 

1 /2 

Ns where Si = {g/2Tx)Ng' Mp. If (pi > (p^, then at pa 50 e-foldings of inflation 
will be given by [53] 

1 \ 2db / 2% ' 

^OidJ \gMpNs' 



Note that the dependence on Ns is quite weak; for the case oi d — A (d — 6) AD 
baryogenesis, 0, oc ATg^/^ (Ng^^^). If there is no large number of inflationary e-foldings 
one can essentially fix the value of 0*. In this case one can predict the magnitude of 
the baryonic isocurvature perturbation. 

Imposing a chaotic- type initial condition V{(f)i) ^ Mp yields 

V2Mp , , 

(t>i ~ . (216) 

Ad-i 

By directly solving the slow roll equations for (f) and 5", we obtain the total number of 
e-foldings of infiation: 

where A^^ is the number of e-foldings during V{(j)) domination, provided (pi > (ps- y{S) 
will dominate the total number of e-foldings only if 

Since Ns > 50, the above condition will be satisfied so long as A is not very small 
(for example, ifARil/(d — 1)!). In this case the value of at the time when the 
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CMB perturbations are generated will be determined mainly by the total number of 
e-foldings of inflation, i.e. Nt ~ Ng. 



5.4.2 Adiabatic perturbations during D-term inflation 

The potential for the flat direction condensate is far from flat, and so if the magnitude 
of the flat direction condensate is large, it will cause a large deviation from scale- 
invariance to the adiabatic perturbation. This will impose an upper limit on the 
amplitude of the flat direction condensate at 50 e-foldings before the end of inflation. If 
we assume that the flat direction follows a late time attr actor trajectory together with 
the inflaton, then following the analyses in [53], and Kawasaki and Takahashi [54], the 
flat direction induced adiabatic density perturbation can be estimated from Eq. (125). 

For a potential of the form V = V{S) + V{(j)), one obtains (with the help of 
Eqs. (84,85,115, 126)) [53] 



V 



Ml 



{Vs + V;)V 
and 



VsVs + V.V, - 



^{Vs + V',){V'sVs' + V;v',^) 



Ml. 



(219) 



2\^ 

(220) 



iy's + V:;)V [ 2V 

Particularly, for the case of D-term inflation, if V"^ < and V'^ < Vg , we obtain the 
conventional result, see Eq. (115). Here the isocurvature contribution to the spectral 
index has been neglected; we will discuss it in the next subsection. 

Since the main contribution to the scale-dependence of the perturbations comes 
from ?7, let us estimate the deviation from scale-invariance due to the presence of the 
flat direction condensate. Note that when V'^ > Vg , with <^ Vg and <^ Vs still 
satisfied, we can expand rj in order to obtain corrections to the conventional D-term 
inflation model [53] 

^«.M||^-M|^. (221) 

The condensate scalar induced deviation from the scale invariance in the spectral index 
is given by 

" VV-^ ■ ^222) 

VSVs 
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Requiring that |An<^| < the present CMB observations imply that n — 1.2 ± 0.3; 
adopting K < 0.2 imposes an upper bound on 0, 



1 

\ 4d-7 c ^ o 4d-15 



< 0, = fc, ( -= 1 g^X^^^Mp''-' , (223) 



where 

(l287r(ci-l)2(2rf-3)J " ^^^^^ 
For the case of minimal d — 4 AD baryogenesis, one obtains [53] 

= 0.53 (^-^ j ' (/A-^e'Mp) ^ ~ 10^6 GeV , (225) 

while for d = 6 AD baryogenesis scenario, one gets 

0, = 0.77 (-^) " (g'X-'eMl) ^ ~ 10^^ GeV . (226) 



5.4.3 Adiabatic perturbations during F-term inflation 

During F-term inflation, the dominant part of the flat direction potential is given by 
(see Eq. (189), Sect. 4.6.1.) 

Vtotai(<l>)^^^^ + V(<l>) , (227) 

where V{(f)) is the usual part from the non-renormalizable superpotential term. Here 
we assume that C/ ~ — (9(1). In such a case the local minimum of the flat direction 
condensate is given by Eq. (190), denoted here by (pm- 

Note that if (p is close to (pm = |0 — 0m| ^ (pm), then inflation will damp S(f) to 
be close to zero. The equation of motion for perturbations around this local minimum 
is given by 

S4> + 3H6^ = -kH'^dcp ; k = {2d - 4)Ci > 1 , (228) 
which has a solution of the form: 

5(f) = 50oe"^* ; a = ^(-3 + V9^^) . (229) 
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As long a,s Ht ^ 1, i.e. there are a significant number of e-foldings, the amphtude of 
the flat direction condensate will be damped to be exponentially close to the minimum 
of its local minimum. 

In general, it is likely that the initial value of 4> will not be close to 0^- It has been 
shown that the deviation of the adiabatic perturbation from scale-invariance implies 
that the value of the potential at A?" 50 cannot be very much larger than (f)jn [53]. 
The deviation from scale-invariance due to the flat direction condensate is then 

For ^ 0m, the field will be rapidly oscillating in its potential and the change in the 
amplitude of over an e- folding due to damping by expansion will be dcp/dN ~ — 0. 
Requiring that |An^| < K imposes an upper bound on 

^ / Kd 1 d-2 ^ ^ 

<^ ~ V2H—^M^ . (231) 



For d = 4, one finds [53] 



while for d = 6 



6{d-l)\\ 

< 0.9 fXMp\^ 



(232) 



(233) 



where we have used K = 0.2. For typical values of Hj, the scale-invariance of the 
density perturbations implies that at ~ 50 e-foldings cannot be much more than 
an order of magnitude greater than 0^. Since there is no reason for to be close to 
this upper limit when A^ 50, it is most likely that will be close to 0^ when the 
primordial perturbations responsible for large scale structure formation have left the 
horizon during inflation. 

5.4.4 Isocurvature fluctuations in D-term inflation 

The isocurvature perturbation in the baryon number arises from the AD scalar if the 
angular direction is effectively massless, i.e. mass is small compared with H during 
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and after inflation [52, 53, 54]. The resulting perturbations will be unsuppressed until 
the baryon number of the Universe is generated. This in turn requires that there are 

no order H corrections to the super symmetry-breaking A-tcrms. 

The baryon number from AD baryogenesis is generated at ~ uisusy ~ 

100 GeV 

when the A-term can introduce B and CP violation into the coherently oscillating 
AD scalar. If the phase of the AD scalar relative to the real direction (defined by the 
A-term) is 9, then the baryon number density given by (see Eq. (253), in Sect. 5.5.) 

ub ~ msusy4>l sin 29 , (234) 

where 0o is the amplitude of the coherent oscillations at « irisusy One can then 
obtain fluctuation in the baryon number or an isocurvature perturbation as 

5nB 259 H 



He tan(2e) 7r0tan(2e) ' 



(235) 



where 69 ~ (if/27r0) is generated by quantum fluctuations of the AD field at the time 
when the perturbations cross the horizon. The magnitude of the CMB isocurvature 
perturbation relative to the adiabatic perturbation can be written as [53, 203, 208], 
(see Eq. (129), in Sect. 3.4.3.) 



a — 



3 \V{S)tsin{29)(l)J ' ^ ' 



where V{S) is the infiaton potential (see Sect. 5.4.1). For purely baryonic isocurvature 
perturbations ui = QB/^m, where fi^ is the baryon density and flm is the total matter 
density. For the case of D-term infiation one obtains [53] 

" " 6^0Ari/2tan(2^) ' ^ ^ 

where N 50. 

Requiring that the deviations from the spectral index due to the AD scalar are 
acceptably small, for d — A, one finds 

3.3o;((7A)^/9 



irV9tan(2^) 



(238) 
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and ior d — 6 

The range of allowed by nucleosynthesis is 0.006 < < 0.036 (for 0.6 < /i < 0.87) 
[3]. For = 0.4, K = 0.2, and for d = 4, we obtain 



and for ci = 6 



= (3.0X10- -0.018)^1^. (241) 



(The lower limits should be multiplied by 0.4 for the case fi^ = 1.) If, for example, 
5( ~ A ~ 0.1 and tan(26') ^ 1, one would obtain a lower bound a ^ 10~^ for d = 4 
and a ^ 10~^ for o? = 6. Such small isocurvature contamination could be detectable in 
future CMB experiments. 

Present CMB and large-scale structure observations require that a ^ 0.1 [203, 208]. 
COBE normalization combined with the value of (the rms of the density field on a 
scale of 8 Mpc) as obtained from X-ray observations of the local cluster together with 
the shape parameter F Q^/i = 0.25 ± 0.05 [333] from galaxy surveys, which is also 
consistent with the recent observations of high-redshift supernovae [20, 21]) yields the 
limit a ^ 0.07. The limit may however rely too much on COBE normalization, which 
is just one experimental result among many. 

Future CMB observations by MAP will be able to probe down to a ~ 0.1, while 
PLANCK (with CMB polarization measurements) should be able to see isocurvature 
perturbations as small as 0.04 [204] (see also [205]). For the case of minimal id — 4) AD 
baryogenesis, there is a good chance that PLANCK will be able to observe isocurvature 
perturbations at least if infiation is driven by D-term. For higher dimension AD 
baryogenesis [d > 6) the situation is less certain. 

All this assumes that can take any value. This is true if 0j < 0*, in which case 
(f) remains at its initial value 0^. We have seen that the dynamics of the flat direction 
during D-term infiation implies that if 0j > (j).^ then will equal 0* at A?^ ~ 50. In this 
case we can fix the magnitude of the isocurvature perturbation. For d = 4, N 50 
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and flrn = 0.4 it is given by [53] 



. = (0.17 -1.03)1^) (242) 



(For — 1 this should be multiphed by 0.4.) For d — 6 and — 0.4, 
« = c...(4.4xl0-3-2.6xl0-)y) 1-^. 



(243) 



If g', A ^ 0.1 then for the d — 4: case one expects a* 0.01 — 0.1. For the d — 6 case 
the isocurvature perturbation might just be at the observable level. 

It is important that one can fix the isocurvature perturbation to be not much larger 
than the lower bound coming from adiabatic perturbations. This is because there is 
typically a very small range of values of over which the isocurvature perturbation 
is less than the present observational limit, a ^ 0.1, but larger than the adiabatic 
perturbation lower bound, a ^ 0.01 for d — A. 

5.4.5 Isocurvature fluctuations in F-term inflation 

If the flat direction condensate is stuck in a local minimum (p 0^ given by Eq. (190), 
the isocurvature perturbation is given by [52, 53] 

Q, ~ f244l 

3 tan{2e)6p(j)m ' 

where 5p — ?,5T/T Ri 3 x 10~^ is the density perturbation. Given H, d, and the value 

of (pmi the magnitude of the isocurvature perturbation is essentially fixed. For d = 4 
and flm — 0.4, the isocurvature perturbation has been found to be 

\ 1/2 / ^ X 1/2 

a = (3.1-18.6) X 10^ —— (-/) (245) 

CySan(2^) \MpJ 

while for d = 6 

\ 1/4 / ^ X 3/4 

a = (2.9 - 17.4) x 10^ — — ( — ^ . (246) 

Cy'tan(2^) \MpJ 

If we require that a ^ 0.1 we find the upper bounds Hj/Mp < 10~^/A (for d — 4) and 
Hj/Mp < 10~^/A^/3 (for d — 6). For typical values of H the isocurvature perturbation 
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in the F-term inflation can be close to the present observational limits. In [54] , however, 
it was pointed out that there would be negligible isocurvature perturbations produced 
from the MSSM fiat directions, and the present observations would not be able place 
any independent constraint upon the initial amplitude of the flat directions. 

5.4.6 MSSM flat directions as a source for curvature perturbations 

In a very recent development, a new paradigm has been laid, where adiabatic density 
perturbations were generated from the decay of the pure isocurvature perturbations. 
Though it was first suggested in [319], but it was implemented recently in a pre-Big- 
Bang scenario [320] , where the axion field which generates isocurvature perturbations 
decays late in the Universe. 

In order to create pure adiabatic density perturbations it is important that the 
field, known as curvaton [321], is subdominant during inflation, but becomes domi- 
nant during the late phase of the Universe especially when it is decaying. The curvaton 
field a generates isocurvature perturbations during inflation, assuming that the per- 
turbations generated from the infiaton field can be negligible, and the curvaton mass 

~ Vo-o- ^ ^inf- this limiting case the power spectrum for the curvaton will be 
equivalent to a massless scalar field during inflation; Va — -f^*/47r^, where * denotes 
the epoch when perturbations are crossing the horizon k — a^H^.. The curvaton field 
follows it trajectory after the end of inflation, and when H ~ m^, the curvaton os- 
cillates and eventually decays through its coupling to the SM relativistic degrees of 
freedom. While oscillating it produces density contrast 5cr = 25a /a, assuming that 

< cr* [321], and the perturbation spectrum is given by 

pV2 ^ 2^ = ^ . (247) 

When the curvaton decays it converts all its isocurvature perturbations to the adiabatic 
ones by following that, before decaying the relativistic degrees of freedom due to 
the infiaton decay products gives rise to density perturbations in the radiation as 
Cr = (l/4)5r, and Co- = (1/3)(5ct- With these results the curvature perturbations is 
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given by [321] 



Apr + Ap„ 3 



supposing that before decaying C,^ is negligible. If the curvaton decay products do not 
dominate the Universe, then there will adiabatic and isocurvature perturbations both 
[321, 322]. 

In [323], the authors have pointed out that the MSSM flat directions which are 
lifted by the non-renormalizable operators such as o? = 7, 9 are the best candidate 
for the curvaton. The flat directions d — 7, LLddd (hfted by H^LLLddd) , and d — 
9, QuQuQue (lifted by QuQuQuHdee), which are lifted by superpotential: W ~ 
a'^~^'4)/M'^~^. Note that is the superfield other than the curvaton, does not produce 
any A-term in the potential, since {ip) = 0, and therefore does not give rise to any 
U{1) violating terms in the flat direction potential. It was shown in [323], that these 
flat directions can dominate the energy density of the Universe and while decaying the 
squarks and sleptons can directly decay in the MSSM relativistic degrees of freedom. 
Therefore the virtue of this scenario is that the MSSM flat direction is solely responsible 
for reheating the Universe, barring any need for speculation from the inflaton coupling 
to the SM fields. Infiation was supposed to happen in the hidden sector of the theory, 
which does not necessarily couple to the SM sector. 

5.5 Baryon number asymmetry 

In both D-and F-term inflation the inflaton and other scalar flelds begin to oscillate 
coherently about the minimum of their respective potential after the end of inflation, 
and the post-inflationary evolution of the flat direction condensate is no exception. If 
Cj and C are positive in Eqs. (189) and (192), the corresponding vevs arc cither given 
by Eqs. (190) or (193). In fact, in D-term inflation models for |C| less than about 0.5, 
it is possible to have a positive correction and still generate the observed baryon 
asymmetry as shown by McDonald [324]. In addition, there has been attempts for 
AD baryogenesis in F-term inflation, basing on the low energy effective action of the 
heterotic string theory, where inflation is driven by T-moduli [334]. Flat directions 
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beyond MSSM involving a triplet Higgs has also been considered [335] . Here we will 
mainly concentrate on the negative correction only. 

After inflation, (0) initially continues to track the instantaneous local minimum 
of the scalar potential, which can be derived by replacing Hi with H{t) in Eq.(190). 
Once if ~ mo, the low-energy soft terms take over. The condensate mass squared 
turns positive, and since the phase of (0) differs from the phase of A, {(f)) starts to 
change non-adiabatically. 

In an absence of H corrections to the A-terms, the initial phase 6 of the AD field 
relative to the real direction is random and so typically ~ 1. As a result (0) starts 
a spiral motion in the complex plane (see Figs. (2), and forthcoming discussion on 
the condensate trajectory), which leads to a generation of a net baryon and/or lepton 
asymmetry [45, 46]^. 

The baryon number density is related to the AD field as 

nB,L = /3^(0V-0V), (249) 

where P is corresponding baryon and/or lepton charge of the AD field. The equations 
of motion for the AD field are given by 

4> + 3H^+^^^ 0. (250) 

The above two equations lead to 

\dV((f)) 

hB,L + SHnB,L = ^pim 

By integrating Eq. (251), we obtain the baryon and/or lepton number as 

a'{t)nB,L{t) = 2/3|a|^ f a' {t'Mit')]" MO) dt' , (252) 

Note that a introduces an extra CP phase which we may parameterize as sin(5). After 
a few expansion times, the amplitude of the oscillations will become damped by the 
expansion of the Universe and the A-term, which is proportional to a large power of 
^There have been attempts for AD baryogenesis in local domains, see [336] 
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(f), will become gradually negligible. The net baryon and/or lepton asymmetry is then 
given by [46] 

nB,L{tosc) = g| rw00o sin 26' sin 5, 

~ ^ 3(^-3) ""^ (m^M'^-^)'^^''"'^ sin 29 sin 5 , (253) 

where sin 5 ~ sin 2^ ^ ^(1)- 

When the inflaton decay products have completely thermalized with a reheat tem- 
perature Trh, the baryon and/or lepton asymmetry is given by [337] 

where we have used H{tosc) ~ fn<j), and s is the entropy density of the Universe at the 
time of reheating. For d = i, the baryon-to-entropy ratio turns out to be [337] 

and for d = 6 

!!M « 5 X 10-0 X /, (^] fi^V' (^f , (256) 

s ' \ m4, j\ j \MpJ VlOO GeV; ' ^ ^ 

where we have taken the net CP phase to be ~ The asymmetry remains frozen 

unless there is additional entropy production afterwards. Note that for d — 4, the 
required reheat temperature of the Universe is below the gravitino overproduction 
bound (see Sect. 3.6.2). For higher dimensional non-renormalizable operators, a low 
reheat temperature is favorable, which is indeed a good news. 

In this regard low scale inflation, which guarantees a low reheat temperature, has 
been given some consideration [338] (see also [339] where AD baryogenesis after a brief 
period of thermal inflation, required to solving the cosmological moduli problem, has 
been discussed). Although, in [340], it was pointed out that in gauge mediated super- 
symmetry breaking it is hard to reconcile AD baryogenesis with a moduli problem. 

Among the host of MSSM flat directions which are lifted by non-renormalizable 
operator and listed in Table 1, the LH^ flat direction carrying the lepton number 
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is the candidate for producing lepton asymmetry in the Universe (there has been 
some earher attempts of direct baryogenesis via udd directions, see e.g. [341]). The 
lepton asymmetry calculated above in Eqs. (255,256) can be transformed into baryon 
number asymmetry via sphalerons ub/s — (8/23)ni/s. AD leptogenesis has important 
imphcations in neutrino physics also, because in the MSSM, the LH^ direction is lifted 
by the d — A non-renormalizable operator which also gives rise to neutrino masses [46] : 



where we have assumed the see-saw relation nii, i — {Hy)'^/Mi with diagonal entries for 
the neutrinos Ui, i — 1, 2, 3. The final ub/s can be related to the lightest neutrino mass 
since the flat direction moves furthest along the eigenvector of LiLj which corresponds 
to the smallest eigenvalue of the neutrino mass matrix [46] . 



where m„ i denotes the lightest neutrino. 
5.6 Thermal effects 

In our discussion on the baryon/lepton asymmetry we have tacitly assumed that the 
asymmetry has been generated before the Universe has thermalized and reheated. This 
might not be the case if there were light degrees of freedom which have thermalized 
with an instantaneous plasma temperature Ti^st < {HT^Mpy/^ before the inflaton 
has decayed. We remind that the bulk of energy density is still in the form of inflaton 
oscillations, and only a fraction of the energy density has gone into these light Standard 
Model degrees of freedom. If the MSSM flat direction couples with this thermal bath, 
there arises a modiflcation in the flat direction potential [309, 310, 342]. 

5.6.1 Thermal corrections to the flat direction potential 

Besides the D-tcrm couplings of the form g'^cfP'a^^ where a is some field with gauge 
interactions, there are also F-term Yukawa couplings to fields x which result in a term 
/i^l^PlxP in the flat direction potential, x ^^id o; obtain large masses due to the 




(257) 




(258) 
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flat direction vev, and therefore do not feel the effect of temperature if the condensate 
amphtude is large. As pointed out by AUahverdi, Campbell and Ellis [309] , the back- 
reaction effect induces a mass-squared term h^T'^ for the ffat direction. If this exceeds 
the negative Hubble induced mass squared term, the ffat direction oscillations starts 
earlier than otherwise expected. In order for thermal correction to play a significant 
role the couplings h, g must have intermediate strength. Otherwise, a large coupling 
would induce a large vev dependent mass for a, which would prevent its thermal 
excitations, and a very small coupling would not have significant thermal backreaction 
at all. 

For the infiationary scale Hj ~ 10^^ GeV and M — Mqut, it has been found [309] 
that a generic MSSM fiat direction with a Yukawa coupling h ~ 10~^ starts oscillating 

H ^ 10^ GeV for 4 < d < 8. Since thermal effects induce early oscillations, baryon 
asymmetry is also produced much earlier, which could have interesting consequences. 

Another thermal effect has been discussed by Anisimov and Dine [310, 342]. All 
the ffat directions which are lifted at large d give rise to a large mass for a, and 
consequently one should account for their effect by integrating out the heavy modes. 
This would result in terms like 

^r^ln^. (259) 



167r2 '^'^ M2 

In particular for the ffat direction LHy_ the effective potential thus obtained has the 
form [310] 

12 



K// = «KT)a,TMn M_ , (260) 

where ctg = gl/'i'K, and ag — |^ + "^y) includes leading order contribution of 
the gluons, gluinos, and quarks to the free energy for a non-abelian group SU (N). The 
oscillations in the flat direction are induced when H^^^ = d'^V^ff /d\(f)\'^ = a^CgT'^/l^p 
and one can check that [309] for d — 4, the thermal mass correction ~ h^T^ wins 
over the logarithmic counterpart Eq. (260), but ior d — 5 and/or 6, the logarithmic 
correction dominates and the oscillations start earlier than otherwise one would have 
expected. 

There could also be a thermal enhancement of the A-term [309], which can arise 
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from the cross terms 

Wz,h<t>aa + Xaj^^. (261) 

However a symmetry forbids such enhancement [342], although the situation might 
change if one adds more terms in the superpotential such as [342] . 

W ^ — [al + h— \ . (262) 

5.6.2 Thermal evaporation of the flat direction 

It has been argued that in general the flat direction condensate decays as a result of 
scattering with the thermahzed decay products of the inflaton [46] . Usually the scat- 
tering interactions preserve B and L, and therefore the previously produced baryon 
and/or lepton asymmetry remains unchanged. In [310] it was assumed that thermal- 
ized fermions scatter with the condensate with a rate 

Tscatt-yg^T. (263) 

where yT corresponds to the mass of the condensate. A complete evaporation was 
found to be avoided only after reheating if T^^ < {yg'^)'^/^Mp^H^/^, which is usually 
satisfied for a reasonable range of reheat temperatures and Yukawa couplings. 



5.7 Baryosynthesis and neutrino mass 

As discussed in Sect. 5.4., the lepton asymmetry via LH^ direction leads to a predic- 
tion on the lightest neutrino mass. It is however pertinent to include also the finite 
temperature effects [343]. At finite T, the flat direction potential potential for LH^ 
direction can be written as [343] 



^a,alT'^.[^^U§-^, (264) 



V„., = [ml-C,H'+ cA^r^j |0p + |22(„,„^4 + h.c.)+^(a„0'' + h.c) 

fk\'l>\<T 



|2\ |J,|6 



where = are real positive constants and couplings = 1 — 10 ^ in MSSM [337]. The 
mismatch in phases between and an leads to the helical motion of the flat direction. 
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Once the inflaton decay products generate a thermal plasma with a temperature T — 
{TrhMpHiY^^, thermal corrections take over the Hubble induced term 

H'<ml+ c,f,T' + a,al(T)^ . (265) 

fk\<l>\<T 1*^1 

The flat direction starts to oscillate when [343] 



Hn^r ~ max 



where Hi is given by [337, 343] 



M 



(266) 



H, - 1;^^' ^^ift^^^ruY"] ■ (267) 
The lepton asymmetry is then given by [343] 

AtWit) ^ 1^ / dt'a^{t')Xuv (am(t>') , (268) 

The right hand side of Eq. (268) initially increases until H ^ Hose, after which the 
integrand is rapidly damped because a^0^ ~ forn > 1. The final lepton asymmetry 
is determined approximately by the configuration at the time when the oscillations 

commence [343] 

-Jm (a™0^) t = lm3/2MHoscSeff , (269) 

^ ^ H=H, 



riL = 



2M 



osc 



where Seff — sin(4arg(0) + arg(aOT)) is the net CP phase. The final baryon to entropy 
ratio turns out to be [343] 



Note that the final expression obtained does not depend upon the reheat temperature 
Trhi mainly due to the fact that the Hose is determined by thermal correction ~ 
T^\a.{\(j)\'^). This is however true only for 10^ GeV < T^h < 10^^ GeV as pointed out 
by Fujii, Hamaguchi, and Yanagida in [343]. For 10^ GeV < T^u < 10^ GeV, the 
dependence on reheat temperature appears as ub/s oc T^l^, because then the thermal 
mass term ~ T^|0p dominates. 

It is possible [343] to obtain the right amount of baryon asymmetry with the lightest 
neutrino mass rriui — (0.1 — 3) x 10~^ eV and with a CP phase Sgff — (0.1 — 1) for a 
fairly wide range of reheat temperature 10^ GeV < Trh < 10^^ GeV. 
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Figure 2: Affleck-Dine condensate formation with x = (pi and y = 4>2i for (a) gravity 
mediated case with d = 4 (sohd) and d = 6 (dashed), and (b) gauge mediated case with 
d = 4, nifj) = 1 TeV (sohd) and = 10 TeV (dashed) with the initial condition 9i = — tt/IO, 
from [344]. 

5.8 Trajectory of a flat direction 

Let us now turn our attention to the dynamical evolution of the flat direction after 
the end of inflation. Here we assume that the flat direction is tracking its minimum 
which is determined by Eq. (190). The trajectory of the flat direction depends upon 
the potentials Eqs. (200) and (204). Here we sketch the main differences between the 
gravity and gauge mediated cases. 

Jokinen [344] has studied numerically the trajectories of the flat direction con- 
densate in gravity and gauge mediated cases, following Eqs. (200) and (204). The 
rotation of the condensate depends on the low energy supersymmetry breaking mass 
terms. The classical motion for the condensate = (0i + ^^2) / is illustrated in 
Fig. (2). In the gravity mediated case. Fig. (2a), we see that the orbit is a spiraling 
ellipse and in the gauge mediated case. Fig. (2b), a precessing trefoil. From Fig. (2), 
one can see that there is a twist on the orbit much before the rotation starts properly. 
This is the time of the phase transition, when the condensate starts to rotate in 
the pit of the symmetry breaking minimum. The rotation begins when the symmetry 
breaking minimum is the vacuum, and ends when it has become a false vacuum, and 
twists when the false vacuum has completely vanished forming a kink on the orbit. It 
is possible to produce a condensate through a second order phase transition, but the 
charge in that case will be small. It should also be noted that in the gravity mediated 
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case condensate formation starts when CjH'^ ~ "^3/2 values of d, A and a. In 

the gauge mediated case the condensate formation starts at CiH^ ~ m^/|(/)p, so that 
the formation happens earlier if the condensate mass is increased, as can be seen 
from the different positions of the kink in Fig. (2b). 

5.9 Instability of the coherent condensate 
5.9.1 Negative pressure 

The effective equation of state of a coherent scalar field oscillating in a potential 
U{(j)) with a frequency which is large compared with H is obtained by averaging 
vl p— (|0P/p) — 1 over one oscillation cycle T. The result is [230] 

p=(7-l)p, (271) 

where 

For the case IJ ~ m^(/)^, one finds 7 = 1, so that one effectively obtains the usual case 
of pressureless, non-relativistic cold matter. 

When the motion of the condensate field is not simply oscillatory, such as in the 
case for the condensate trajectory, one can generalize Eq. (272) by integrating over 
the orbit c of the AD field. In that case 

^ ^ 21A4>\ {\-U{4>)Ip)^ .273) 
/,d|0|(l-t/(0)/p)-^ ■ 

In practice the orbits are nearly eUiptical. Then the arc length is given by 



(274) 



^\ A\\-4>\IA^) ' 

where A and B < A are respectively the semi-major and the semi-minor axis of the 
ellipse, and 4>\ = Re4>/\/2. For a circular orbit B — A, whereas for pure oscillation 
(no charge in the condensate) 5 = 0. 

It is therefore obvious that small corrections to a harmonic potential of a coherent 
condensate can easily generate a pressure. As we have seen, in the gravity mediated 
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Figure 3: Prcssure-to-energy density ratio, p/p = ttj = 7 — 1, in the gravity mediated case 
vs. (a) time in logarithmic units for cZ = 4, (b) different initial conditions for d = 4, 6; (c) 
ellipticity e = B/A vs. initial conditions for cZ = 4 (thin lines), d = 6 (thick lines), D-term 
(solid), F-term (dashed) with K = —0.01 and t = lOOm^l^^- (b) w is shown att = SOOm^l^^ 
with dotted lines for the d = 4 D-term case, from [344]. 

case quantum corrections typically modify the flat direction mass terms by 

Ui'P) = \ml<t>' + Kmy log + ■ ■ ■ ^ (275) 



where K is some constant. If one writes 



1 /^2\3; 
,2 i2 



U{<l>) = 2^^r [j^ j (276) 
one finds that 

In the case of the logarithmic potential x ~ 2K. There arises a negative pressure 
p — Kp whenever K < or whenever x is small and negative. This is a sign of an 
instability of the condensate under arbitrarily small perturbations. 

This is exactly the situation one finds in the MSSM with flat directions. The 
effective mass ^^^(0) = dU/d(j)'^ decreases for a range in 4>, albeit for different reasons, 
both for gravity mediated and gauge mediated supersymmetry breaking. According 
to Eq. (273) this results in a negative pressure, which has been computed numerically 
by Jokinen [344]. The results are shown in Figs. (3) and (4). 

The pressure-to-energy density ratio; w = 7 — 1 for the gravity mediated case is 
plotted in Fig. (3a), which shows that there are time-dependent oscillations in pressure 
[344] . The average pressure is slightly on the negative side. The average value of w is 
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Figure 4: Pressure-to-energy density ratio, w, in the gauge mediated D-term case (without 
Hubble induced ^-term) vs. time in logarithmic units for (a) d = 4 and (b) d = Q; (c) 
ellipticity of the orbit where d = A (thin lines) and d = 6 (thick lines). The scalar masses 
= 1, 10, 100 TeV are denoted respectively with solid, dotted and dashed lines, from 
[344]. 

shown in Fig. (3b) at t ~ lOOmg^g ^ ^^'^ different initial conditions. In Fig. (3c), 
the ellipticity of the orbit, e — B/A, is plotted to show that w is more negative if e 
is small. It should be noted that w achieves values which are more negative than the 
absolute lower bound coming from pure oscillation. 

A similar analysis has been made for the gauge mediated supersymmetry breaking 
case. The results have striking similarities [344]. In Figs. (4a) and (4b), the time 
development of the pressure-to-energy density ratio, w, for = 4 and d — Q has been 
depicted. One can see that the pressure is always negative. The calculation of average 
pressure is even more involved than in the gravity mediated case, since the oscillation 
frequency becomes very large. In Fig. (4c), the ellipticity of the orbit is shown as 
a function of different initial conditions. Jokinen [344] has pointed it out that quite 
generically £ < 0.1. 

5.9.2 Growth of perturbations in the AD condensate 

As a result of internal negative condensate pressure the quantum fluctuations in the 
scalar condensate grow according to [9] 



5k = -Xk^^k 



(278) 
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Figure 5: Energy-to-charge ratio, x, in the gravity mediated case vs. (a) time in logarithmic 
units for d = 4, 6; (b) the D-term case; (c) the F-term (with Hubble induced ^-term) case 
with d = 4, 5,6, 7 (solid, dash-dot, dashed and dotted lines), K = —0.01 and t = lOOm^^^^ 
from [344]. 



U K < 0, quantum fluctuations of the condensate fleld at the scale A = 27r/|k| will 
grow exponentially in time as 

(5(/)k(t) = (5(/)(0)exp (^-KkH) . (279) 

In reahty the onset of non-hnearity sets the scale at which the spatial coherence of the 
condensate can no longer be maintained and the condensate fragments. For the AD 
condensate the initial perturbation originates from inflation. Note that since the AD 
condensate carries a global charge, due to charge conservation the energy-to-charge 
ratio changes as the the condensate fragments. 

The energy-to-charge ratio has been estimated numerically for both gravity and 
gauge mediated cases by Jokinen [344]. The time evolution of the energy-to-charge 
ratio X is shown in Fig. 5, where x is also plotted for various initial phases in F-and 
D-term inflation models. For the gauge mediation case the plots are quite different 
from the gravity mediated case, see Fig. (6). 

5.9.3 The true ground state 

Under the negative pressure the homogeneous AD condensate fragments and forms 
lumps. The question then is, what is the true ground state? The answer is, a non- 
topological soliton with a flxed charge, called the Q-ball [345, 346, 47], which in general 
is made up of a complex scalar fleld with a global U(l) symmetry, for which the 
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Figure 6: Encrgy-to-chargc ratio, x, in the gauge mediated case vs. (a) time in logarithmic 
units of time c/ = 4, (b) o? = 4 and (c) d = 6 with D-term (thin hues, without Hubble induced 
A-term)and F-term (thick lines, with Hubble induced A-term) and = 1, 10, 100 TeV 
(solid, dashed, dotted lines) at t = 4 • 10^m^\ 10^"z^\ 4 • lO^m^^ (d = 4) and i = 4 • 
10^ 10^ 4 • 10^ m^^ (d = 6), from [344]. 

Lagrangian is 

jC = di,(f)d^(f)* - U{(f)(f)*) . (280) 

When supplemented by the CP violating terms, this is the Lagrangian for the MSSM 
flat directions. 

The conserved current is — i{(t)d'^(t)* — (f)*d^(f)), and the conserved charge, and 
energy are given by 

Q = jd'xf, (281) 

E = I d^x[\^\^ + \V(j)\^ + U{(l)(l)*)] (282) 

If the charge is kept fixed, the state of lowest energy is found by minimizing [47] 

E^^E- cu{Q - i J d^xf) , (283) 

with respect to variations in cf) and the Lagrange multiplier uj. From 5^E^ = 0, it 
follows that (j) ~ i^4' = 0, so that we may write 

(/)(i,x) = e^'^V(x) (284) 

where </? may be chosen real by virtue of U{1) invariance. The charge and energy of 
such a configuration read 



Q ^ 2u J d^x ip^ , 



(285) 
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E ^ j d?x [{y^f + U{^'') + , (286) 

and one has to minimize = J (i^a;[(Vv?)^ + Uui{f^)] + uiQ where 

UUip'')^Uiip^)-u;V ■ (287) 

To find a locahzed configuration that vanishes at spatial infinity one may make use 
of the spherical rearrangement theorem, which implies that E^j is minimized by (p{'x) 
which is spherically symmetric and monotonically decreasing. This is equivalent to 
solving the equation of motion 

(f(p 2 dip dU^ ^ f288) 

If by convention we set the globally symmetric minimum to (f = with U{0) = 0, one 
can then show that a non-trivial solution to Eq. (288) is obtained whenever U{ip'^)/ip^ 
has a minimum at </? 7^ 0, i.e. U {(f"^) grows more slowly than m^</7^ over some range. 
We will discuss Q-balls in detail in the next Sects. 6 and 7. 



5.10 Numerical studies of fragmentation 

Although the homogeneous AD condensate is not the ground state, it is not obvious 
that the ground state should always be reached within cosmic time scales. It is then 
essential to study the dynamical evolution of the AD condensate. Since the Q-ball 
formation is inherently a non-linear phenomenon, analyzing small perturbations is 
not sufficient to determine the full dynamical evolution of the AD condensate. One 
can nevertheless gather some information about the gross features of the condensate 
fragmentation by perturbative considerations alone. 

5.10.1 Perturbation theory 

Negative pressure is equivalent to an attractive force between the condensate quanta 
which induces a growing mode in spatial perturbations. A linearized description of 

the evolution of perturbations has been given by Kusenko and Shaposhnikov in [151], 
and by Enqvist and McDonald in [48, 347, 49, 328]. For a maximally charged AD 
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condensate {B — in Eq. (274)), the linearized perturbation takes the form = 
(p{t) + 6(j){x, t), and 9 — 9{t) + 59{x, t), where the homogeneous condensate is described 

by 

$ = ^e^^W , (289) 
v2 

with (j){t) — {ao/a)^^'^(j)o and 9{ty ^ m^. Initially S(f){x,t) and 69{x,t) should satisfy 
the relationship [151, 48, 49] 

The solution of the linear perturbation equations then has the form [48] 



a. 



^^«^-j''%0„exp(/d<(ii!^j je- (291) 
and 

f r /l \K\m'i\ ^^'^\ 
^»«^».exp^/*(^--,^j Je-. (292) 

For the gravity mediation case, the above condition applies if |k^/a^| ^ \2Km'^\, and 
is small compared with and ^ 1. If the first condition is not satisfied, then 

the gradient energy of the perturbations produces a positive pressure larger than the 

negative pressure due to the attractive force from the logarithmic term, preventing the 

growth of the perturbations. 

For the case of a matter dominated Universe, the exponential growth factor is then 

[48, 49] 

[,j\%mf.im%f. (293) 



2\ 1/2^ 



2 a' 9{ty 

The largest growth factor will correspond to the largest value of for which growth 
can occur, 

\ ^2\K\ml. (294) 

max 

The value of H at which the first perturbation goes non-linear is [48, 49] 



H. « ^ . (295) 
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with 



„(A) = -logf^l . (296) 



where (po is the value of when the condensate oscillations begin at if ~ m<^. A 
typical value of a (A) (e.g. with d = 6) is a (A) pa 30. The initial non-linear region has 
a radius Aj at H^, which is given by [48, 49] 

For the case of a non-maximally charged condensate the situation is slightly differ- 
ent. It is likely that the initial radius and the time at which the spatial perturbations 
initially go non-linear will roughly be the same [49] as for the maximally charged con- 
densate. In general, the charge density of the initial non-linear lumps will essentially 
be the same as that of the original homogeneous condensate. 

The perturbative evolution of a single condensate lump was considered in [49]. In 
terms of = (0i -|- i02)/\/2, the initial lumps are described by 

0i(r,t) = Acos{m^t){l + cos{7ir /ro)) (298) 

(p2{r,t) = Bsm{m^t){l + cos(7rr/ro)) , (299) 

for r < To and by 0i^2 = otherwise. The initial radius of the lump is 2ro, where 
^0 = T^/iV^l^l^^"^^,/))- The maximally charged condensate lump corresponds to A = 
B, while the non- maximal lump has A > B. The total energy and charge in a fixed 
volume are given by [328, 329] 

E^Att j drr^p , Q ^ An j drr^q - AB, (300) 

with Qmax = A'^- 

In [49], the behavior of the solutions was found in a perturbative analysis to de- 
pend on K, and to a greater extent on Q/Qmax- The condensate lump was found to 
pulsate while charge is flowing out until the lump reaches a (quasi-) equilibrium pseudo- 
breather configuration, also called (J-axiton, with the lump pulsating with only a small 
difference between the maximum and minimum field amplitudes. For the Q-axiton, 
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in which the attractive force between the scalars is balanced by the gradient pressure 
of the scalar field, the energy per unit charge is much larger than m<^; indeed, the Q- 
axiton exists even if Q = 0. Only for a maximally charged Q-axiton are the properties 
similar to that of the corresponding Q-ball. It is however unclear whether Q-axitons 
are just an artifact of perturbation expansion. 

5.10.2 Lattice simulations 

The features of the fragmentation of the AD condensate cannot be fully captured 
by studying various mean field theory approaches, such as in large N-approximation 
and Hartree-approximation [239, 238, 243]. The formation of a Q-ball is a non-hnear 
process for which various mode-mode interactions become important. This can be 
seen by expanding the perturbed and as shown by Kasuya and Kawasaki in [348, 
349, 350] 

. . 

5(t) + ?,H8(i)- 29(1)5(1) - —5(t) + U"U)5(t) = 0, 

(f)9 + 3H(f)S9 + 2((j)S9) -2'i^eS(p- (p^59 = 0. (301) 

Although the potentials differ in the gauge and gravity mediated cases, it is never- 
theless always possible to identify the fastest growing amplified mode. In the gravity 
mediated case we have already obtained that by inspecting Eq. (293). A similar anal- 
ysis can be performed for the gauge mediated case by noting that C/"(0) Ri — 2m^/0^. 
Taking into account the conservation of charge OcfP'a^ = const., along with the approx- 
imation of a circular orbit, one may simplify Eq. (301) by seeking a solution of the 
form S(f) — 6(f)Q exp(Q;i -|- ikx) and 89 — 59q exp(Q;i -|- ikx). In order to further simplify 
the analysis, one can also assume a — const, and (p — (po — const., so that the phase 
velocity 9 — (C/'/0)^/^ ^ v^m^/^o- H « is real and positive, the fluctuations grow 
exponentially and become non-linear. Solving for 6(j)Q, 69q, Kasuya and Kawasaki finds 
for gauge mediated case [351] 

Note that in order for a to be positive, one must require the expression in the second 
parenthesis to be negative. This means that the instability band for the fluctuations 
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is given by 

< - < — ^ . (303) 
a 00 

The most amplified mode appears at {yia)max ~ (3/2)^/^m^/0o in the gauge mediated 
case. 

Various groups have studied the fragmentation of the AD condensate and the 
formation of Q-balls numerically. In [330], condensate fragmentation was simulated 
numerically on a 2 + 1 dimensional 100 x 100 lattice, starting with a uniform AD- 
condensate with 0o = 10^ GeV and with an arbitrary phase cu, 4> — 0oe*'^* + S4> with 
uniformly distributed random noise S(j) ~ O(lO~^^)|0o| added to the amplitude and 
phase. The parameter values chosen for the simulations were = 10^ GeV, K — 
—0.1, and A = 1/2. The results indicate that first the charge density of the condensate 
decreases uniformly due to the expansion of the Universe. As time progresses a growing 
mode can be seen to develop. White noise is still present but the growing mode soon 
starts to dominate. This process continues until lumps of positive charge develop. 
These are dynamically arranged in string-like features but the filament texture is a 
transient feature which disappears soon, see forthcoming Figs. (9,10). 

The further evolution of the fiat direction depends on the initial energy-to-charge 
ratio of the condensate, defined as 

x=^, (304) 
mQ 

and hence on the value of a;. If a; = 1, i.e. the energy-to-charge ratio of the condensate 
is equal to that of a Q-ball, no negatively charged Q-balls, anti-Q-balls, are formed. 
After the modes grow non-linear, the lumps just evolve into Q-balls and finally freeze 
due to the expansion of the Universe. 

If a; < 1 so that x ^ 1, the fragmentation process has a much more complicated 
history [330]. After the positively charged lumps have formed, expanded linearly and 
then developed non-linearly, the extra energy stored in them causes the lumps to frag- 
ment as they evolve into Q-balls. In this process a large number of negatively charged 
Q-balls forms. The total charge in the negative and positive Q-balls is approximately 
equal so that the initial charge in the condensate is in fact negligible compared to the 
amount of charge and anti-charge created. 
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Figure 7: Q-ball formation in gravity mediated supersymmetry breaking. Prom the left, the 
first plot shows the amplitude of the condensate after Q-balls have formed at 2; = 6.3. The 
second plot shows around 40 Q-balls with a largest charge Q ~ 5.16 x 10-^^, and the third 
plot shows the power spectra of the condensate fluctuations {k\6ipk\'^) when the amplitude 
of the fluctuations has become as large as the homogeneous mode: {Sif"^) ~ ip'^. The top 
panels (a) and (c) show the full fluctuations calculated on one dimensional lattices, while 
the bottom panels (b) and (d) show the linearized fluctuations without mode mixing, from 
[350] 

Full 3 + 1 dimensional simulations have been presented by Kasuya and Kawasaki 
[348, 349, 350] (for both gravity mediated and gauge mediated cases), see Figs. (7,8), 
and by Multamaki and Vilja [352] (for a wide range of the energy-to-charge ratio in 
the gravity mediated case with), see Figs. (9,10). 

In [348], the authors simulated Q-ball formation in the gravity mediated scenario 
on a (64)^ lattice with a lattice spacing = 0.1. The results are shown in Fig. (7). 
The initial fluctuations in the real and the imaginary direction were taken to he 5 (pi — 
54)2 ~ 0(10-^), together with initial field values ^1(0) = ^2(0) ~ O(IO^) GeV. In 
order to obtain the spectrum the authors relied on a Id lattice with N = 1024 and 

= 0.1. The result is depicted in the third plot of Fig. (7) for two different comoving 
time scales. One can see the marked difference between the linearized perturbations 
and the lattice simulated ones. In the latter case the spectrum does not fall sharply, 
which can be attributed to mode-mode interactions or rescattering effects which kick 
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Figure 8: Configuration of Q-balls on a three dimensional lattice in the gauge mediated case. 
More than 30 Q-balls have formed, and the largest one has the charge with Q 2± 1.96 x 10^^. 
The second plot shows the power spectrum of the condensate fluctuations when the amplitude 
of the fluctuations has become as large as the homogeneous mode. The top panels (a) and 
(c) show the full fluctuations calculated on one dimensional lattices, while the bottom panels 
(b) and (d) show the linearized fluctuations without mode mixing, from [349]. 

the lower momentum modes higher, leading to a broadening and smoothening of the 
spectrum. Note that in linearized fluctuations the instability band is almost the same 
as Eq. (294). For example k/m^ = V2a{T)\K\^^^ ^ 2 for \K\ = 0.01 and r = 5.5 x 10^. 

In [349], Kasuya and Kawasaki repeated their simulation for the gauge mediated 
case, which is shown in Fig. (8). Note that the size of the Q-ball is bigger than in the 
gravity mediated case. 

Multamaki and Vilja [352] had typical lattice sizes of 120^. They verifled that the 2 
dimensional simulations [330] capture all of the essential features of the AD condensate 
fragmentation: transient filament structures resulting in a large number of Q-balls and 
anti-Q-balls, which can also be seen in Fig. (10). Note that when the condensate has 
the exact energy-to-charge ratio of a Q-ball so that x — 1, no anti-Q-balls form in 
Fig. (9), whereas for x ^ 1, the number of Q-balls and anti-Q-balls are practically 
equal. 



118 




Figure 9: Formation and evolution of charged lumps on 3D lattice from comoving time 
r = 875 — 3000, when x = 1. Only positively charged Q-balls have formed, from [352]. 




Figure 10: Formation and evolution of charged lumps from comoving time r = 1500 — 7500, 
when X = 10^. Here positive and negatively charged Q-balls have formed, depicted in 
different shades, from [352]. 
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5.11 Equilibrium ensembles 

Some of the lattice results can be understood by analytical arguments. In particular, 
the distributions and number densities of Q-balls and anti-Q-balls may be obtained 
simply by maximizing the entropy. Such approach appears justified, since after frag- 
mentation, the AD lumps arc expected to interact vigorously and the field fragments 
will settle to the state of lowest energy by emitting and exchanging smaller fragments. 
If the interaction is fast enough compared with the expansion rate of the Universe, i.e. 
r = ritotO'v > H — (2/3i), where ritot is the total number of Q-balls and anti-Q-balls, 
a PS ttRq is the geometric cross-section of a Q-ball coUision, and v is the average 
velocity of a Q-ball, then it is naturally expected that the final state should consist of 
an equilibrium distribution of Q-balls and anti-Q-balls in a state of maximum entropy. 
Considerations supporting this has been given in [330] for the case of gravity mediated 
supersymmetry breaking. 

The Q-ball (anti-Q-ball) distributions N^{Q,p) {N_{Q,p)) are subject to the fol- 
lowing constrains: 

Etot^E+ + E_, E±^JdQdpE{Q,p)N±{Q,p) 

Qtot = Q+-Q-, Q± = JdQdpQ N±{Q,p) , (305) 

where E{Q,p) y/p"^ + w?Q'^ is the energy of a single Q-ball, £"+ {E_) and (5+ {Q-) 
are the energy and charge of Q-balls (anti-Q-balls), and Etot and Qtot are respectively 
the total energy and charge of Q-balls and anti-Q-balls, which are equal to the energy 
and charge of the initial AD-condensate (unless significant amounts of energy and/or 
charge are transformed into radiation). It then follows from Eq. (305) that 

_ m (g+ + g_) 

X = Etot > 7^ > 1- (306) 

mQtot 

This condition is independent of the precise form of the Q-ball distributions. 

If all the baryon asymmetry resides in Q- and anti-Q-balls, then at times earlier 
than 10-^, Qtot/Q+ ~ AS ~ 10"^. (Since B - L is conserved in the MSSM, this 
holds also for the purely leptonic fiat directions.) Prom Eq. (306) it follows 

X ~ 10^ . (307) 
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Even if all of the baryon asymmetry were not carried by Q-balls, a natural expectation 
is x » 1 so that the number of Q-balls, iV+, and the number of anti-Q-balls, iV_, are 
approximately equal and the total number of Q-balls is Ntot ~ 2iV_|_. The main bulk of 
the Q-balls may be expected to be relativistic [330] , and the collision rate much larger 
than the Hubble rate. 

Indeed, one can verify that the equilibrium assumption is self-consistent, and it 
is also supported by numerical studies in 2 dimensions [330] as well as by the 3 di- 
mensional simulations of Multamaki and Vilja [352], who actually observe the Q-ball 
distribution to relax into equilibrium as shown in Fig. (10). 
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6 Q-balls 



Q-balls with a global charge have many interesting astrophysical and cosmological 
consequences. They may reheat the Universe [353, 354] and serve as candidates for 
dark matter candidate [151, 355, 356, 351, 357, 358]. They may provide baryogenesis 
and leptogenesis [48, 49, 347] and while decaying, produce LSP dark matter [359, 
360, 351, 361, 362]. They could be responsible for the generation of cosmic magnetic 
fields [363]. Q-balls could stabihze neutron stars [364] or even form solitonic q-stars 
[365, 366, 367]. They could act as a laboratory for physics beyond the electroweak 
scale [331]. 

6.1 Q-ball as a non-topological soliton 
6.1.1 Proofs of existence 

The Q-ball is a generic ground state in a broad class of theories with interacting scalar 
fields carrying some conserved global charge [345, 346, 368, 369, 47, 365]. The Q-ball 
is an example of a non-topological soliton whose boundary condition at infinity is the 
same as that for the vacuum state, unlike in the case of topological solitons such as 
magnetic monopoles [370, 371] (a detailed review on non-topological solitons can be 
found in, e.g. [372, 365]). 

Q-balls are quite generic solitons in 3 -|- 1 dimensions, which can be associated 
with many scalar fields with various U{1) charges [50], with a non-Abelian symmetries 
[373, 374, 375], and also with local gauge symmetries [376, 377, 378, 356, 379, 380, 381]. 
The main difference which distinguishes global Q-balls from a local Q-balls is that in 
the latter case the charge of the stable Q-ha\\ is bounded from above. 

Recall that in order to find a Q-ball solution one must minimize the energy E^^ 
(see Eqs. (285,286)) for a fixed charge with respect to the variations of lu and (p{x) 
independently. Obviously one could try finding directly a Q-ball solution by solving 
Eq. (288). A more effective way is to seek a bounce solution (p{x) for tunneling in three 
Euchdean dimensions in the potential given by Eq. (287) [382, 384, 158]. Note that 
the first term in Eq. (286) is then the three dimensional Euclidean action S^lifujix)] of 



122 



this bounce solution which satisfies Eq. (288) in radial coordinates with the boundary 
conditions (p{0) = (p'{oo) = 0. 

A theorem [47, 365] (see also [382]) states that if there exists a range of ip however 
small, and the potential U {(p^) contains an attractive interaction however weak, then 
a non-topological soliton solution exists for 

i/^ < tu^ < , (308) 

where the mass parameter is defined as U{ipP) — > m^(/?^ when — > 0. The value of 
uj = ^Jk^ + m| determines the frequency of the (p quanta in the field space. For a plane 
wave solution it is always true that a;^ > m| while for a solitonic solution a;^ < m^. 
This suggests that there exists a parabola v'^ipP' tangent to U{ip'^) at (p — ±(/?o, with 
< m^. Another useful way of expressing this is through 

= min , for = (^0 > . (309) 

For a sufficiently large Q given by Eq. (285), the energy of a sohton is then given 

by 

E = \pQ\ < m^\Q\ , (310) 

which ensures its stability against decay into plane wave solutions with (/? ~ (/7o inside 
and ~ outside the sohton. Note that the global U{1) symmetry is thus broken 
inside the soliton by the vev, however, remains unbroken outside the soliton. The 
most crucial piece is the presence of a global U{1) charge of the Q-ball which actually 
prevents it from decaying and makes the soliton stable. (For analytical results on 
Q-ball properties, see [383]). The above discussion can be repeated in the presence 
of several charges gj, and an analogue of Eq. (309) has been estabhshed by Kusenko 
in [50]. The main difference between a single and multi-charged Q-ball is that there 
exists different vevs corresponding to different charges, which modifies the appropriate 
bounce solution Eq. (309). 

6.1.2 Beyond thin wall solution 

The above discussion tacitly assumes a thin wall approximation where the edge of a 
soliton is sharply defined, which means that the gradient energy is small compared 
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to the volume energy. As (/7o increases, the thin wall limit breaks down for a fixed 

charge, and one must seek other methods in order to guarantee the existence of a 

Q-ball solution as pointed out by Kusenko [382]. For a flat potential, which mimics 

the large field value situation, the equations of motion Eq. (288) can be solved near 

the origin r — 0, and for large r. The Q-ball profile can be found to be 

, , sin(a;r) „ 

^{r) = ^0 , r<R, 

cur 

= (^le-"*^^ r>R, (311) 

where the values of ipo^ipi^cu and R are such that they minimize E^, while ip{r) in- 
terpolates r — R, the size of the Q-ball, continuously. In a thick wall limit one can 
also write down Ei_j in terms of dimensionless variable ^ = cux and — (p/cu, while 
neglecting all the terms in except the constant term and the cu'^^p'^ term [382, 158], 
one obtains E^^ ^ acu + b/uj^ + ujQ, where a, b are constants independent of uj. The size 
of the Q-ball turns out to be R (x 1/cu, and the vev of the field in the Q-ball interior 

IS (fio oc 

from Eq. (311). Note that there is no classical hmit on the charge of a 
Q-ball. In fact no matter how small Q is there always exists a value cu close to m^, 
for which E^^ is minimized. Quantum stability requires that Q > 1. When Q ^ 1, 
the quantum corrections will indeed become important. For Q > 7 the Q-balls are 
quantum mechanically stable configurations [386]. 

The Q-ball has been shown to be classically stable by Coleman [47] (see also [365, 
382]). The semiclassical approach obviously breaks down when quantum fiuctuations 
around the Q-ball are comparable to the energy of the system itself. This happens 
when Q ^ 1. Even though the charge becomes small, the Q-ball size remains large 
in comparison to the Compton length m^^ [382, 385, 386]. One can also establish a 
Virial theorem, which holds for any Q and does not require any approximation [50]. 

6.2 Varieties of Q- balls 
6.2.1 Thin wall Q-balls 

There exist thin wall Q-hall solutions, where the boundary is a well defined edge, as 
well as thick wall Q-balls, where the boundary is not localized in a narrow region and 
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the soliton is typically described by a Gaussian profile. Both Q-ball types may appear 
within the same theory. 

The thin wall Q-ball is the simplest and arises naturally in any suitable scalar 
potential that allows for the existence of a Q-ball. As mentioned already, thin wall 
solution has the profile (f{r) ~ (fioS{r — R) in the radial direction, where R is the size 
of the radius of a spherically symmetric Q-ball. This obviously neglects the surface 
energy contribution and yields [47] 



Energy is thus growing hnearly with charge. Note that the radius of such a Q-ball can 
be very large. 



These are useful relationships for the purposes of this Section. 
6.2.2 Thick wall Q-balls in the gauge mediated CEise 

Thick wall Q-balls have been widely considered in the literature within gauge and 
gravity mediated supersymmetry breaking scenarios. In both cases the mass of the 
Q-ball grows more slowly compared to a thin wall case, i.e. the scalar potential grows 
slower than 9?^. In this case Q-ball never reaches a thin wall regime, even if Q is large. 
The value of ip inside a Q-hall extends as far as the gradient terms allow, and the mass 
of a Q-ball is proportional to Q^, where p < 1. 

In the context of gauge mediated supersymmetry breaking the AD potential takes 
the form (without the non-renormahzable contributions) [151, 331] (see Eq. (204), in 



where m,^ ~ 0(TeV) represents the supersymmetry breaking scale. In [387], the 
authors have considered the effective potential of the form 




(312) 




(313) 



Sect. 5.3.2.) 




(314) 




(315) 
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where F^/^ :$> m^. Despite the difference between the forms of Eqs. (314) and (315), 
the dynamics of the flat direction is similar to the one given in Eq. (314), which yields 
the equation of motion 

if H — (p — —ijj </? . (316) 
r 

For large r, (/^(r) ~ exp(— m^r), where is the mass of near the origin. The 
interpolating solution was already presented in Eq. (311). The energy of such a Q-ball 
grows as [151] 

The profile of the Q-ball is given by (^(r) ~ exp(— m^r). The radius and the value of 
the vev inside the Q-ball is roughly given by [151] 

R ^ -j^Q^/\ (318) 

'"'^gV4. (319) 



6.2.3 Thick wall Q-balls in the gravity mediated case 

If the potential grows only slightly slower than cp'^ as in the case of gravity mediated 
supersymmetry breaking scenarios, the potential may be approximately written as 
[48, 347, 49] (see Eq. (200), in Sect. 5.3.1.) 

,2 



U{(p) ftimi l + iriog 



M2 



ip' , (320) 



where K < 0, and M is the largest mass scale. Note that at small vevs we have again 
neglected the non-renormalizable contributions in the above potential. The mass scale 
is given by m<^ ~ 1^3/2 ~ 0{TeY). 

The Q-ball equation of motion is written as 

/ + V^-^oV + ^Jv'^logJ^^j , (321) 

where cuq is defined by 

ujI = uj^- ml (l + K) . (322) 
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For thin wall Q-balls, the initial value of </? is very close to (pd, the value of (p for 
which the right-hand side of Eq. (321) vanishes. In this case (p will remain close to (pci 
up to a radius of order luq^ \og{ipci/S(p{0)), where 6ip{0) = {(pci — V^(0)). It will then 
decrease to zero over a distance Sr ~ coq^, corresponding to the width of the wall of 
a Q-ball. The radius of a thin wall Q-ball can be made arbitrarily large by choosing 
50(0) small enough. 

For a thick wall Q-ball, the initial value of ip can be much smaller than (pci- In 
this case the non-renormalizable terms may be neglected. In general, the right-hand 
side of the Q-ball equation of motion vanishes for three values of (p, which correspond 
to (pci, (pc2 and zero. (pc2 corresponds to the point at which, assuming that the non- 
renormalizable terms can be neglected, the first two terms on the right-hand side of 
Eq. (321) cancel, and we obtain [49] 



ipc2{r) = Me' 

In deriving the last correspondence we have used Eq. (322). Note that (pc2 is an 
attractor, in a sense that if (p{0) is close to (pc2 it will tend towards (pc2 as r increases. 
The radius of a thick wall Q-ball in the gravity mediated case is given by [48, 347, 

49] 

^ ^ 1^11/2 ' (324) 

where R is defined as the radius within which 90% of the Q-ball energy is found, and 
[48, 347, 49] 

uo ^ \K\^/^m^ . (325) 

Since typically is small compared with 1, we find uj ^ m^. In the gravity mediated 
supersymmetry case the size of a Q-ball does not depend on charge, unlike in the 
gauge mediated case, see Figs. (7,9,10), where the sizes of the Q-balls are all equal. 
One may take a Gaussian ansatz for the profile of a thick wall Q-ball [49] 

cp{r) = cp{0)e~^ , (326) 
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provided one identifies uj^ ^ 3\K\m1 and R'^ ^ 2(\K\m1)~^. 

The total charge of the Gaussian thick wall Q-ball is given by [49] 

3/2 



dr ATrr'^ujiple'l^ = (-] uj<^IR^ (327) 
while the total energy is given by [49] 

E^l{lY\lR+{lY'rnyoR\ (328) 

where the second term in the above equation is the combined contribution from the 
potential energy and the charge term, where we have used o;^ ~ m^. Since R is large 
compared with m^^ for small \K\, the potential plus charge term dominates the energy. 
The radius within which 90% of the energy lies is then given by Rc = 1.25R. The 
energy per unit charge is given by [49] 

where wc have used the Gaussian result cjg = Sli^lm^. For all practical purposes we 
can take E m^Q. Although the energy per unit charge is larger than m^, the mass 
of the scalar at small values of will have the form m^{l ^ ol\K\) (with a ^ 1) once 
the logarithmic correction to the potential is included, so that the binding energy per 
unit charge will be positive and of order \K\m^. 

The last two examples Eqs. (314,320) exhibit two extremes of any thick wall type 
Q-ball. Any thick wall Q-ball should belong somewhere in between, such as in hybrid 
case. 

6.2.4 Hybrid case: gauge and gravity mediated Q-ball 

In Sect. 6.2.2., we discussed the Q-ball potential in the gauge mediated case, but it 
is true that any generic flat direction is also lifted by gravity mediation as well. In 
the gauge mediated case the full flat direction potential relevant for Q-ball formation 
should read as [349, 351] 

12' 



^W) = "^^l0g 1 + ^ +rn3/2|¥' 
Tflj 



l + |is:|loe 



m_ 

Ml 



(330) 
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where 7713/2 takes values between 100 KeV and 1 GeV. The gaugino loops lead to 
K < as discussed before (see Sect. 5.3.1.), but Yukawa couplings give rise to K > 0. 
On the other hand if the Yukawas dominate, the AD condensate can be stabilized and 
Q-balls can only form for ip < ipgg. The second term in the above potential dominates 
when 



<p><Peq = . (331) 

7773/2 

There are two distinct regimes. When (p > (peq, the Q-ball properties resemble the 
gravity mediated thick wall case. Otherwise, when (p < (peg, the Q-ball properties are 
similar to the gauge mediated thick wall case. The energy per unit charge can be 
written as [349, 351] 



Q 



(332) 

7773/2 <P><Peq- 

Obviously in between there should be a hybrid regime which interpolates smoothly 
between the gauge and gravity mediated cases. 

6.2.5 Effect of gravity on Q-balls 

In principle gravity can give a significant contribution to the Q-hall energy, as shown 
in a study by Multamaki and Vilja [388]. In a thin wall case, the interesting result 
is that gravity limits the maximal size of a Q-ball. The reason is that besides the 
rotational motion in the complex field space which generates outward pressure, there 
exists a gravitational attraction. Since the gravitational contribution to the Q-ball 
energy is negative, it is possible that gravity can render an otherwise unstable Q-hall 
stable. The effect of gravity on Q-balls remains small provided the soliton is much 
larger than the Schwarzschild radius and the charge smaller than the gravitational 
charge Qg ~ (Mp/m^)'^, which is quite large Qg ~ 10^^ for 777,^ ~ 100 GeV. 

6.2.6 Q-balls and local gauge inveiriance 

So far we have considered Q-ball solutions in theories with a global U{1) symmetry. 
The symmetry group can however be extended to include global non-Abelian symme- 
tries [373, 374, 375, 380]. The existence of Q-balls in a supersymmetric Wess-Zumino 
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model has been demonstrated by Axenides, Floratos and Kehagias [375], who showed 
that Q-balls form domains of manifestly broken supersymmetry. 

A gauged Q-ball has some additional interesting properties [376, 378, 356, 379, 381]. 
Taking a complex scalar field 0(r, t) — /(r, t) exp(— i^(r, i))/-\/2 coupled to an Abelian 
gauge field A^, the charge of a given field configuration is given by [376] 



where e is the gauge charge. Since the gauge field inside a Q-ball is broken by the non- 
zero vev of (p which couples to the gauge field canonically, the gauge field is massive 
and acts as a U{1) superconductor, provided the Compton wavelength of the gauge 
field is smaller than the size of a Q-ball. As a result of broken gauge symmetry, 
there is an extra source of electrostatic self-energy contribution which comes from 
the electrostatic repulsive force due to the presence of a gauge charge. The gauged 
charges are also repelled to reside on the boundary of a Q-ball. In a gauged Q-ball, for 
a fixed charge Q, both radius and energy are relatively large. The remarkable feature 
is that there is a maximum allowed charge and correspondingly a maximum radius. 
This is due to the repulsive electrostatic potential Aq which tends to destabilize the 
gauged Q-ball. The maximum charge-to-radius ratio is governed by the gauge charge: 

Qmax/Rmax = 47r(m0 - UJc)/e'^. 

In [379] , it has been argued that the presence of fermions could stabihze the gauged 
Q-ball. During Q-ball formation, fermions could be trapped inside the Q-ball. Those 
with charges equal to the charge of the scalar quanta would be repelled from the inside, 
whereas fermions with opposite charges would remain and render the Q-ball neutral. 
Large electric fields inside the Q-ball could also lead to a pair production and to a 
subsequent screening of the charge. 

An interesting application of a gauged Q-ball could be the hadronic structure of 
QCD, as in the Friedberg-Lcc model [389], where hadrons are modeled by phenomeno- 
logical non-topological solitons. Gauged Q-balls in theories with a Chern-Simons terms 
has also been considered in [377]. 




(333) 
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6.3 Q-haW decay 

6.3.1 Surface evaporation to fermions 

In the MSSM, the scalar field forming a Q-ball can interact with fermions. Then a 
Q-ball can decay into a pair of fermionic quanta. Q-hall decay has been considered in 
[390, 391], where fermion production was studied in a classical background of a Q-ball. 
As it was first pointed out in [390], for a large Q the Q-ball evaporates through its 
physical surface and there exists an upper bound on the evaporation rate per unit 
area. This behavior has been verified by numerical studies [391]. 

In order to understand this, let us imagine that a region inside a Q-matter with 
a vanishing 0, which might appear due to fermion pair production, forms a cavity. 
Suppose L is the linear size of the cavity. The energy of Q-matter formerly inside the 
cavity is given by [390] 

Eq - u;l4>lL^ , (334) 
while the charge Q within that region is given by [390] 

AQ ~ uocplL^ . (335) 

N massless fermions inside the cavity will have a free Fermi gas distribution, and 
therefore, the energy of the fermions is given by [390] 

[-J-) ■ (336) 



The ratio of energies is given by [390] 



En I huj, 



2\ 1/3 
1 



(337) 



E^ \ (t)l j 

In the semiclassical limit ^ ^ 0, cavitation is energetically forbidden irrespective of 
the size of a cavity. In a sense, it is the Pauli exclusion principle which keeps the 
Q-ball stable. Fermions are produced but there exists a Fermi pressure which prevents 
further production. Inside Q-matter fermions gain mass of order g(f)Q and saturate the 
Fermi energy so that fermions can be produced only from the surface. 

Following [390], let us assume that each fermion carries the energy ~ huJo/2. A 
simple bound on the average fermionic pair production from the surface can then 
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be obtained by noting that maximum current density can be reached only when the 
outward moving fermions from the Q-wall are occupied while every level in the inward 
moving fermions is empty. Then, by assuming massless Wcyl fermions with a single 
helicity state for each k, one obtains the limit on the outward moving current [390] 

(n • j),<,„/2 < k^dk cos{e)d{cos{e)) d4> = ■ (338) 

Integrating this expression over the enclosing surface yields the maximum evaporation 
rate per unit area in the limit R oo [390] 

"^"^ <7^- (339) 



dtdA - 1927r2 

A detailed calculation relies on estimating the reflection and transmission coefficients 
of the ingoing and outgoing waves outside the Q-ball and matching these with the 
solutions obtained inside the Q-ball. In order to calculate the transmission coefficient, 
one has to sum over infinite sequence of scatterings. In [390], it was shown that 
the maximum transmission coefficient is 1 when ooq < g(f)Q. The weak coupling limit 
"^0 > 5*00 leads to a different transmission coeflicient and the rate of evaporation is 
given by [390] 

dtdA " Uo ; V1927r2 ) ~ 32n^ ' ^'^^^^ 
The factor g(f)Q determines the penetration width of the fermions. In other words, in 
the weak coupling limit fermions can penetrate deep inside a Q-ball without completely 
filling the Fermi sphere. Within Pauli exclusion principle, Q-matter could decay into 
very weakly coupled fermions within the whole Q-ball volume. This will be discussed in 
the context of L-balls in Sect. 7.1, which can decay into massless neutrinos throughout 
the interior. 

In [390, 391], it was shown numerically that the evaporation rate is strongly de- 
pendent on R but approaches the limiting profile given in [390]. In a realistic case a 
step-function is not always a good approximation, in particular for a thick wall Q-ball. 
In the thick wall case the problem has been investigated numerically by Multamaki 
and Vilja [391], who found that for a sufficiently large Q the evaporation rate decreases 
with increasing Q. As a result the Q-ball evaporates faster when its size decreases. 
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6.3.2 The decay temperature 

For a thin wall Q-ball the surface area is related to the charge via A = (367r)^/3g2/3^^^2/3^^ 
where k — (2(/7qC/((/?o))^'^^. The lifetime of a Q-ball in this case is given by [49] 

V 3 y 



1447r ( - ) . (341) 



For a thick wall case, the area of a Q-ball is independent of its charge, being fixed by 
its radius R {\K\^/'^m^)~^ . The Q-ball lifetime in this case is then given by [49] 

From the above expression, one can estimate the temperature at which the Q-balls 
decay By assuming radiation domination, the decay temperature is defined as [49] 

T.^{^fm'\ (343, 



where kx = (^^^^^^) • By taking uj ~ a; 100 GeV, thick wall Q-balls will decay 
at a temperature [49] 

lb f u \i/2/ioi5\i/2 



/ uj 10'° \ ' 



where we have set kx ~ 17. The Q-balls will decay at a temperature less than 100 GeV 
if g > 2 X 10^^|X|-^ For thin wall Q-balls, the right-hand side of Eq. (344) has an 
additional factor (Q/Qc)^^^, where Qc is the value at which the thin wall limit becomes 
valid. 

6.3.3 Q-ball decay into a pair of bosons 

If the Q-ball forming scalar is not the hghtest, as might be the case in the MSSM, it 
is also possible that the Q-ball decays into hghter bosons. The decay into scalar fields 
within the Q-ball volume is not blocked by the Pauli exclusion principle. In the case 
of bosons one should replace N^^^ by N in Eq. (336), because the bosons can condense 
in the lowest mode of the cavity. The ratio Eq. (337) becomes 

l^-cuoL. (345) 
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Cavitation in the case of bosons is always energetically favorable for sufficiently large 
L. Q-ball decaying into light scalars could be significantly enhanced relative to the 
decay into fermions. 

The decay into light scalars will only be possible near the edge of a thick wall 
Q-ball. This is because particles coupling directly to the condensate scalars will gain 
a large effective mass from (0) inside the Q-ball. As a result, decay into light scalars 
will only occur via loop diagrams with rates suppressed by the large effective mass. 
This tends to make MSSM 5-balls long-lived, as will be discussed later. 

Not all the flat directions would have a scalar decay mode. An example is udd 
direction, which is lifted at d = 6 and contains the right-handed squarks. For the 
universal boundary condition for squarks at large scales, RG flow analysis suggests 
that the left-handed squark masses are typically heavier than the right-handed ones. 
The Higgs scalar mass could also be heavy compared to the right-handed squark mass 
and the slepton masses usually come out to be lighter than right-handed squarks. 
Hence the decay of such baryonic direction would be kinematically forbidden. Even if 
the right-handed squark decays to Higgses and to sleptons, it would certainly involve 
either a pair of (light) quarks, gauginos, Higgsinos or leptons in the final state. As a 
result the core of a udd Q-ball would be Fermi suppressed. Any final state involving 
a pair of fermions can arise only from the surface. 

Generically Q-balls made up of left-handed squarks and sleptons, such as the d — & 
MSSM flat directions dQL or eLL are expected to decay into a pair of light bosons. 
Note that these particular directions are good candidates for carrying B charge. One 
may parameterize by fs the possible enhancement factor of the scalar decay rate over 
the fermion decay rate, so that [49] 

\ / boson \ / jermion 

This gives for the decay temperature [49] 

\ A&TikrQ ) \\K\) VlOOGeV; \Q ) ^ ' 

For dQL or eLL one may estimate the largest possible enhancement factor using 
the Gaussian thick wall ansatz as shown in [49]. Within the Q-ball, for » m^, the 
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lowest possible dimension operator which could allow the condensate scalars to decay 
at one loop level to hght particles is the one lifted by the d — 5 operator 



^ / d'e^X^rj , (348) 



M 

where x ^^'i V represent the light particles and M fa g(f, where g is the coupling of 
the heavy particles to (p. The decay rate of the condensate scalars to hght scalars will 
then be given by [49] 

dQ 



dt 



-J uf'^{r)r{r)A7rr^dr , (349) 



where uj^p^{r) is the charge density within the Q-ball, and [49] 



r(r) ; g(l)>m^ 

~ am^j, ; < , (350) 

with a — gi^/(47r) (for simplicity we consider a single coupling constant g). Let r* be 
the radius at which ifi{r) — m/g. Then the largest contribution to the decay rate will 
come from a region of width 5r ~ E?/{Ar^) around r*, over which (/? has a roughly 
constant value ip ~ m/g, where 

r* = 7i? ; 7 = In^/^ i '-^] • (351) 

Prom Eq. (349), the rate can be deduced as [49] 
dQ 



dt 



(rrt roo \ 

m^a / drr^ + 0^(0) / drr^e~^ j (352) 

I I \ / \ / Jermton 

where we have used u ~ m^. For a thick wall Q-ball and for typical values of the 
parameters, g(p{0)/m (0.1 — 0.01)Q^/'^, so that 7 ft! 4.5 in Eq. (353). One finds that 
the enhancement factor is typically [49] 

1 70 

if g is less than 1 then for \K\ ft 0.01 — 0.1, wc expect an enhancement factor not 
much larger than about 10^. For g less than 1, most of the enhancement factor comes 
from the unsuppressed tree level decays occurring at r > r*. 
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6.4 Cosmological formation of (^-balls 

Cosmological formation of Q-balls was initially proposed by Kusenko in [50], by 
Kusenko and Shaposhnikov in [151], and by Enqvist and McDonald who studied Q- 
balls from the fragmentation of the fiat direction condensates in the gravity mediated 
case [48, 49]. Subsequently, the problem has been attacked numerically by Kasuya 
and Kawasaki in case of gravity mediation [350, 351] as well as in gauge mediation 
[349, 351], and by Multamaki and Vilja in the gravity mediated case [352]. 

6.4.1 In gravity mediated case 

In the gravity mediated case the size of the Q-ball depends on the charge; the larger 
the charge, the larger is the size of a Q-ball, see Fig. (8). In a cosmological context it 
is natural to think that the Q-balls form when the the most amplified AD condensate 
mode is as large as the horizon size just after the AD field starts rotation, such that 
~ uj^^ ~ ipo/m^. The charge of a Q-ball should be given by [351] 

l+D 

where D — 1,2,3 stands for the number of spatial dimensions, and ui ~ m^. This 
naive expectation has been verified numerically on lattice by Kasuya and Kawasaki 
[351], although the formation turns out to be slightly delayed [350]. The maximum 
charge of a thick wall Q-ball can be written as 

Qma. = /9d<^o^-° , (356) 

where [3ds are some numerical factors with /3i ~ 0.1, /32 ~ 0.02, and /Ja ~ 6 x 10~^ 
[351]. 

Kasuya and Kawasaki [351] also noticed that the charge of a Q-ball depends on 
the helical motion of the AD condensate, and it is proportional to e 

where e = 1 corresponds to a circular motion, and £ = for the radial motion. The 
numerical calculation in [350, 351] indicates that Qmax is constant for small e where 
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Q ^ H-^n^ ^ uj-^uj^l ^ml'^l^] , (355) 



both positive and negative Q-balls with charges of the same order of magnitude are 
produced, while hnearly dependent on £, around £ ~ 1. It was noticed that dominantly 
positive Q balls were formed. Numerical simulations also reveal the presence of small 
negatively charged Q balls [350, 351]. 

6.4.2 In gauge mediated case 

The potential for the condensate $ forming a Q-ball is given by Eq. (204). If the 
AD condensate couples directly to particles in a thermal bath, the potential receives 
a thermal mass correction of order ~ T^l^l^ while integrating out the heavy modes 
yields a contribution ~ iT^ logd^l^/T^). In the latter case the actual sign depends 
upon the integrated modes. If the integrated heavy modes mainly belong to the matter 
multiplet then the sign comes out to be positive, otherwise if it were dominated by the 
heavy gauge degrees of freedom then it turns out to be negative. In a particular case 
of LHu flat direction this sign is positive as shown by Fujii, Hamaguchi and Yanagida 
in [343]. Combining these gives the relevant part of the effective potential (without 
non-renormalizable terms) as [351] 



Here we have assumed that the flat direction obtains a positive contribution to the 
thermal potential. The effective mass can be written as [351] 



Note that T oc t'^l^ during inflaton oscillations dominated phase. 

At very large amplitudes of the AD condensate, the gravity mediation effects for 
supersymmetry breaking dominates and a stable Q-ball of hybrid type, as discussed 
in Sect. 6.2, can form [349, 351]. In this case the AD potential is dominated by the 
terms 



Since here the curvature of the potential depends weakly on the amplitude, the AD 




(358) 




(359) 




(360) 



condensate starts rotating when if ~ 7713/2 <S 1 TeV. In [351], the authors have 
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simulated the dynamics of the AD condensate on 1,2, 3-d lattices, and verified the 
formation of Q balls for X = —0.01. 

In a cosmological context one may estimate the maximum charge of a Q-ball from 
the fragmentation of the AD condensate. The analysis is similar to the gravity medi- 
ated case, except for the fact that the times when the AD condensate starts oscillating 
and the most amplified mode (or instability band) enters the horizon do not coincide. 
The charge of a Q-ball is then given by [349, 351] 

Q ~ Hf^m^^j ~ i\K\'/'m^)-''ms/2\K\^l ~ lirr^/^^^ f— V ' (361) 

\'>TT'(t>J 

where the subscript '/' denotes the time when the Q-ball forms, and we have assumed 

7713/2 ~ rricj) in our final expression. Note that Q oc (/^q, with the proportionality 
constant to be determined numerically; one finds Qmax ^ Ps^'^ [351], where /Ss ~ 
6 X 10"^. 

6.5 (5-ball collisions 

The dynamics of any extended object carrying charge is quite different from the dynam- 
ics of charged point-like objects. In this respect studying the Q-ball collision is impor- 
tant in order to understand the charge distribution of Q-matter in the Universe. Q-ball 
collisions have been studied by a number of authors [392, 151, 350, 393, 394, 395, 396]. 
In [351, 395], collisions have been considered in the context of gravity mediated super- 
symmetry breaking, while in [396], Q-ball collisions have been tackled within gauge 
mediated supersymmetry breaking. An animation of Q-ball collisions can be found at 
www. utu.fi/r^tuomul/ collision, rnpg. 

The usual ansatz is to take initially two spatially well separated Q-balls so that 
initial field configuration (in Id) is [394] 

$(t, X) = e*-i*+^V,,(|a; + a\) + e^^^Vc^d^ - «!) • (362) 

The two Q-balls have the profiles ip^^-^ and (/p^^j) correspondingly to the frequencies uoi 
and UJ2-I and they are separated by a distance 2a. The total charge of the configuration 
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(in Id) is given by [394] 

/+00 
ipu,i{\x + a\)iptj2{\x - ci-Ddx . (363) 
-oo 

Note that the last term is exponentially small in the separation parameter a, because 
the profiles of the individual Q-balls die off exponentially outside Q-matter. Usually, 
for a single Q-ball the phase is unimportant because of the global U (1) symmetry, but 
for multi-Q-ball case the relative phase a plays an important role and affects the total 
charge of the configuration. 

There are two extreme cases: the two charges are equal; or they have opposite 
sign. The key parameters are the relative phase, the incident velocity, and the charge 
[393, 394]. The generic interaction for two Q-balls of equal charges is attractive if the 
relative phase a — 0, and repulsive when a — tt . In case of attraction the two Q-balls 
coalesce to form one larger Q-ball with a resultant charge less than the sum total of 
the individual charges. A loss in charge also occurs when the Q-ball suffers a large 
distortion. 

If the initial phase a 7^ 0, tt or if the charges of the Q-balls are not equal, then the 
dynamics of the Q-ball coUisions result in charge transfer. Q-balls tend to repel each 
other, which happens even after charge transfer. If the incident velocity is extremely 
high (relativistic), then Q-balls simply pass through each other without losing much 
charge [394, 393]. 

A collision of a Q-ball and an anti-Q-ball exhibits several interesting features. A 
naive expectation would have them annihilating. Instead, they bounce back or pass 
through each other. Charge is partially annihilated, though. The main reason is the 
fact that generically Q-balls can transfer their charges only very slowly. The charge 
transfer is very seldom complete [394]. 

These conclusions hold mainly for thin-wall Q-balls. There is not much difference 
between gauge and gravity mediated cases. Because in gravity mediation the Q-ball 
size is smaller, in a fixed volume the AD condensate tends to break into larger number 
of Q-balls [350]. As a result the Q-balls can have larger peculiar velocities than in 
the gauge mediated case. Multamaki and Vilja have studied the gravity mediated 
case in [395] for 2-d lattice for a range of velocities between v — 10~^ and v — 10~^. 
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The authors allowed all possible values of the relative phase and found that the fusion 
cross section of two Q-balls appears to be smaller than the geometric cross section, 
whereas the cross section for the charge exchange is larger than the geometric cross 
section. The probability for charge exchange processes increases with increasing cu. 
The actual velocity of the Q-ball is an open issue, which however, to some extent, can 
be determined in a cosmological context [350] . In the gauge mediated case the peculiar 
velocity of the Q-balls is small and the main interactions are either elastic scattering 
or partial charge exchange [396]. 

6.6 (5-balls in a thermal bath 

There are several effects which one must take into account when Q-balls are immersed 
in a thermal bath at a temperature higher than the Q-ball formation scale. The 
temperature dependent effective potential might not even allow for a Q-ball solution, 
but if it does, then for a large Q thermal corrections are neghgible as argued by Kusenko 
and Shaposhnikov [151], and Lainc and Shaposhnikov [397]. Finite temperature effects 
always lead to an erosion of the condensate and therefore to a loss of charge from the 
Q-ball into the ambient plasma. As a result a chemical potential /ipiasma ~ ^Q/iVT"^) 
arises in the surrounding plasma, where V is the Hubble volume and T is the ambient 
temperature. As pointed out in [397] , the process of Q-ball evaporation will stop when 
the chemical potential fiQ ~ m^{T){Q — 5Q)^^^^ associated with the Q-ball becomes 
equal to fipiasma- The conclusion is strictly valid for the flat potential studied in [397], 
but should hold quahtatively for all kinds of Q-balls. 

In most cases of interest Q-balls are produced non-adiabatically with a formation 
time scale much shorter than the evaporation scale. Then it is natural to ask how 
Q-balls come in thermal and chemical equilibrium. Since the Q-ball does not lose its 
charge significantly, chemical equilibrium between Q-balls and thermal plasma may 
never be reached, in which case /iq always dominates over /ipiasma- It is possible to 
obtain a thermal equihbrium at least between the soft edge of the Q-ball and the hot 
plasma. 

In general there are three different thermally induced effects: dissociation, diffusion. 
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and evaporation of the Q-ball. 



6.6.1 Dissociation 

If Q-balls never reach thermal equihbrium with the plasma, then Q-ball dissociation 
by the bombardment of thermal particles as discussed in [151, 49] is important. In 
case of charge dissociation, the thermal particles in the plasma collide with the Q- 
ball and may even penetrate inside. The penetration width depends on the kinetic 
energy of the particles. For definiteness, let us focus on the gauge mediated case with 
thick wall Q-balls. A particle which interacts with (/? receives a mass contribution 
~ g{(fi{x)), while outside ~ gT. At an ambient temperature T, the particle 
cannot penetrate the Q-ball beyond Xgt, known as stopping radius and determined by 
gip{xst) pa 3T. If %Ij imparts sufficient energy to the Q-ball in order to overcome its 
binding energy within the dynamical time scale then the Q-ball may simply break up. 
On the other hand, if the energy is dehvered to the Q-ball is below the dissociation 
limit, the Q-ball will be able to radiate the excess energy away adiabatically and will 
not dissociate. 

The rate of dissociation depends on the flux of the incoming thermal particles 
/ = (5f*(T)/7r^)47rx^^T^, and the energy per thermal particle transferred to the Q-ball 
~ ^tT, where 7t < 3. Then the rate of energy imparted to the Q-ball is given by [49] 

dE ^ Ag,{T)^TT'(3^R^ ^^^^^ 

dt TT ' 

where Xgt = PR is defined by the Gaussian thick wall profile ip{r) — ip{0) exp(— r^/i?^). 



and therefore /3 = ^J\og{gip{Q)/2>T). In order to evade complete dissociation /S.E <^ 
H^m^Q, where Am^ ~ \K\'m^. The dissociation will not be completed provided the 
temperature of the thermal bath is given by [49] 

A J^, a2 ^^Q'"^ (365) 

where the dynamical time scale is assumed to be ~ kr/m^, with kr > 1- 

In a realistic case dissociation alone cannot erode the Q-ball completely. The Q-ball 
will rather come into thermal equilibrium with the ambient plasma. In an expanding 
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Universe the minimum energy is configured in such a way that the Q-charges are 
always present in the Universe along with the Q-balls. In this case the energy of the 
Q-balls decreases as the temperature of the Universe decreases. 



6.6.2 Diffusion 

Diffusion takes place only through the soft edge of the Q-ball at a distance over which ip 
does not change much. There are two factors which determine the diffusion rate; firstly, 
how efficiently the edge of the Q-ball diffuses, and secondly how fast the core of the Q- 
ball readjusts itself in order to compensate for the loss of charge. At large temperatures 
the diffusion rate is large as the Q-ball tries to relax into chemical equilibrium with a 
thermal plasma. The net diffusion rate is given by [398, 349, 351] 



where D — A/T is the diffusion coefficient with A — 0{1)^ and //q ~ a; is the chemical 
potential of a Q-ball. When the temperature of a thermal bath drops due to expansion 
of the Universe, surface evaporation rate takes over. 

If there is a thermal bath already prior to reheating, the instantaneous temperature 
of the plasma would be large: T ~ {MlViHyl^ ~ (MpT^i/) where Vi is the decay 
rate of the inflaton field. During this period, before reheating, t oc T"^, while after 
reheating t oc T"^. The diffusion rate can be obtained from Eq. (366) [351] 



6.6.3 Evaporation at finite T 

In a thermal bath the surface evaporation rate is no longer given by Eq. (339), since 
one must take into account of thermal corrections. At finite T the evaporation rate of 
a Q-ball has been found to be [397, 349, 351] 



dt 



AttDRq/iqT^ ~ -AttAT , 



(366) 




(367) 
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where jipiasma <^ //Q ~ ~ m{T)Q^f^^ while C, and m(T) are given by 

2 



m{T) 



T 



T 
1 



(T < m^), 



(369) 



The evaporation rate can be calculated as [351] 



r 



2^4 
3 



(370) 



In order to determine which rate is dominating, let us consider the ratio Rdiff = 
^difff^evap- For T > m^, the ratio is given by Rdiff = AQ'^^^. If R^iff < 1, the 
diffusion rate is the bottle-neck for the charge transfer. This condition is met when 
the Q-ball charge is large enough: 

Q>W'{^)\ (371) 

On the other hand when T < m^, the condition Rdiff < 1 corresponds to 

T>% = ^l/3m<^Q-^/'^ (372) 

and the transition temperature is lower than for large enough Q-ball charge. 

If there exists a thermal bath prior to reheating then there arc additional compli- 
cations regarding the evaporation rate, which now exhibits four different possibilities 
depending on how the temperature compares with the reheating temperature and 
the mass of the AD particle [351, 49]. The time-temperature relationship changes at 
T = Trh, while the rate dQ/dt changes at T = m^. Combining all the effects, one 
obtains [351] 

io^^^gV4 (r>r,,,m^). 



'dQ^ 
dT. 



evap 
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{Trh <T < m^), 
{m^<T<Trh), 
{T < m^,Trh). 



(373) 
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The estimate for the loss charge turns out to be of similar magnitude in all possible 
regimes with [351] 

for any case. 

In the gravity mediated hybrid case the evaporation and diffusion rates have the 
same forms in terms of Q-ball parameters Rq ~ \K\~^/'^m-^i2, and u ~ ■^3/2- The 
transition temperature at which Vf.^ap — ^diffi reads = A^/^\K\^/^{rn2,/2'ni'^^Y^'^ . As 
in the 'usual' type of Q balls where the potential is dominated by the logarithmic 
term, the charge evaporation near is dominant and the total evaporated charge is 
found to be [349, 351] 

i^Qr^l^-^(^y'"(^Y'\ (375) 
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7 Cosmological consequences of Q-balls 



In this section we discuss various issues concerning Q-ball cosmology. If the MSSM 
fiat direction carries some combination of B and/or L, the charge will be stored in 
Q-balls created by the fragmentation of the initial condensate. If Q-balls are stable, 
like in the case of gauge mediated supersymmetry breaking, then they will be a good 
candidate for dark matter which can be searched directly. If Q-balls eventually decay, 
the charge will be released in a form of baryonic quanta, providing an interesting alter- 
native mechanism for baryogenesis which does not necessarily depends on sphalerons 
transitions. The evaporation of Q-ball also gives rise to supersymmetric dark matter. 
This is an added advantage of Q-ball cosmology: it provides a physical mechanism for 
relating the dark matter and the baryon densities. 

7.1 L-ball cosmology 

Q-balls only carrying a leptonic charge are known as L-balls. They emerge from 
the LHu fiat direction which is quite different from the rest of the MSSM fiat direc- 
tions. Even though LH^ flat direction might not fragment as already mentioned in 
Sect. 5.3.1., there are choices for the initial conditions at Mqut which make it possible 
to obtain a decreasing mass with decreasing /i in the RG equations given in Eq. (201). 
It has been noticed that becomes negative for scales typically smaller than 10^ 

GeV or so [48, 49, 347]. Depending on the choice of parameters there can be a "hill" 
in the plot of Tn\jj^ versus \ip\ (representing LHu), such that starts decreasing with 
increasing \(p\ for sufficiently large values of \ip\. The effect of negative at small 
enough \(p\ will generate a minimum for C/(|0|)/|0p as required for L-ball formation, 
typically at |</?o| ~ 1 TeV. Effectively such a potential can be given by [49] 

2 

C/(0) ^ ^(2e-^ - 1)^' , (376) 

where s ^ 1 TeV~^. This gives rise to thick wall L-balls with radius it! fa rn'^^. The 
charge of a L-ball is bounded, i.e. L < m^s~'^V, which tends to become zero once 
^ ^ Since L < (sm,/,)"^, the L-balls will have a maximum charge, which for 
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typical values of s and m<^ cannot be larger than 10^ (with an essentially fixed radius). 
This stands in contrast to other thin wall Q-balls for which the charge is proportional 
to the volume. 

Inside L-balls the field strength is of order TeV, which is much smaller than 
the initial amphtude ior d — A AD condensate at if m^. This suggests that L- 
balls cannot form by the collapse of an unstable condensate at H m^. Therefore 
the AD baryogcncsis along the LHu direction will be essentially unaltered from the 
conventional scenario. Even if there were a primordial formation of L-balls, these 
objects would decay at ^ 10^ GeV. It is therefore unhkely that such L-balls could 
have any cosmological consequences. Though it is possible that L-balls, which have a 
field strength of order 1 TeV or less, could play a role in the physics of the electroweak 
phase transition (we will discuss about phase transition aided by solitons in Sect. 8.3.1). 

7.2 B-ball cosmology 

Apart from LH^ flat direction, there are purely baryonic directions such as udd. Q- 
ball forming along this direction carries only baryonic charge and is dubbed as B-hall. 
There are also the dQL and eLL directions. One may expect these directions to 
be phenomenologically similar to udd direction. i?-parity conservation allows d — A 
non-renormalizable superpotential term (HuL)'^ and d = 6 term (udd)'^. In addition, 
there is also the d = A, B — L conserving uude direction (and phenomenologically 
similar QQQL direction). Although this will not produce a B asymmetry via a direct 
decay of the AD condensate in case there is subsequent anomalous B + L violation 
(i.e. sphalerons, see Sect. 2.3.4.), but it can generate a baryon asymmetry via Q-ball 
decays occurring after the electroweak phase transition [48, 49, 347]. 

For large enough B, a 5-ball cannot decay into the lightest 5-carrying fermions 
(the nucleons) , and so it is completely stable. Stable S-balls could have a wide ranging 
astrophysical [151, 364] and experimental implications [331], as will be discussed in 
Sects. 7.4 and 7.6. 

In the gravity mediated case the 5-balls may decay at temperatures below Tg^, 
whence the observed baryon number will be a combination of baryon number originat- 
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ing from the decay of the S-balls and baryon number from free squarks left over after 
the break-up of the squark condensate. Since the S-balls are composed of squarks, 
when they decay they will naturally produce a number density of neutralinos which is 
of the same order of magnitude as the number density of baryons, as was first pointed 
out in [48]. If the S-balls decay sufficiently below the freeze-out temperature of LSPs, 
and if the number density of thermal relic neutralinos is less than that from the S-ball 
decay, then the dark matter density and baryon number densities in the Universe will 
be naturally related. 

The actual ratio of baryons to dark matter will essentially be determined by two 
variables: (i) the mass of the neutralino LSP, and (ii) the proportion of baryon number 
trapped in B-balls as compared to baryon number in free squarks; this is usually 
referred to as the efficiency of S-ball formation. 

When a B-hall decays, for each unit of B produced, corresponding to the decay of 3 
squarks to quarks, there will be at least three units of i?-parity produced, correspond- 
ing to at least 3 neutralino LSPs (depending on the nature of the cascade produced by 
the squark decay and the LSP mass, more LSP pairs could be produced) . Let N^^^S 
be the number of LSPs produced per baryon number and /s be the fraction of the 
total B asymmetry contained in 5-balls. Then the baryon to dark matter ratio is 
given by [48, 49], see also [399], 

^B = = , 377) 

PDM N^fBm^ 

where m„ is the nucleon mass and is the neutralino LSP mass. It is rather natural 
to have < 1. The present LEP lower bound on the neutrahno mass in the MSSM 
(assuming no constraints on the scalar masses) is 17 GeV [400]. If we were to assume ra- 
diative electroweak symmetry breaking and universal masses for the squarks and Higgs 
scalars at the unification scale, then the lower bound would become ^ 40 GeV for 
tan/? < 3 [400]. For > 3, and with > 17 (40) GeV, we find that tb < 1 occurs 
for /b ^ 0.02 (0.008). As long as more than 2% of the baryon asymmetry is trapped 
in S-balls, the observed dominance of dark matter in the Universe can be naturally 
explained. 
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Primordial nucleosynthesis [3] bounds the density of baryons in the Universe to sat- 
isfy 0.0048 < l^ij/i^ < 0.013, where OA < h < 1 (we adopt the bound based on " reason- 
able" limits on primordial clement abundances [3, 57]). The observed baryon to dark 
matter ratio, ~ 0^/(1 — O^) (assuming a flat Universe), satisfles 0.005 ^ ^ 0.09. 
This can be accounted from S-ball baryogenesis, provided [48, 49] 



For example, if the LSP mass satisfies 17 (40) GeV ^ ^ 500 GeV, then the 
observed baryon to dark matter ratio can be achieved by a wide range of fs, i.e. 
0.007 < fsiN^/S) < 3.9 (1.7). Note that if > 67 GeV, then we must have Jb < 1, 
implying that the observed baryon asymmetry must come from a mixture of decaying 
S-balls and free baryons. 

This all assumes that the asymmetry not trapped in the 5-balls can survive down 
to temperatures below Tg^j. If we were to consider & B — L conserving condensate or 
additional L violating interactions in thermal equilibrium above Tg^j,, then the only B 
asymmetry which could survive anomalous B + L violation is the one associated with 
the B-halls. In this case Jb would be effectively equal to 1 (we refer to this case as 
"pure" S-ball baryogenesis (BBB)), so rriy^ would have to be less than 67 GeV. 

A crucial assumption in all these is that there is effectively no subsequent annihi- 
lation of LSPs coming from i?-ball decays. 

7.3 B-balls in gravity mediated supersymmetry breaking 

The attractive force due to logarithmic radiative correction term in the condensate 
scalar potential is given by Eq. (200). Of particular interest is the d — 6 udd squark 
direction with a non-renormalizable superpotential term of the form [uddY and the 
d — A uude direction, which conserves B — L. The magnitude of K is important for 
numerical estimates. Prom the 1-loop effective potential [34], for the udd direction, 
the correction due to gauginos with supersymmetry breaking masses Mgusy is given by 




(378) 
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(see earlier discussion in Sect. 5.3.1.) [48, 49] 

1 a 

K^-ii E (379) 

a, gaugmos cp 

where the sum is over those gauginos which gain a mass from the condensate scalar ip. 
The main contribution will come from the three gluinos which gain masses from the 
squark expectation values. For ^ 0.1 we obtain \K\ f» 0.004(M3/m)^. Depending 
on the ratio of the supersymmetry breaking gluino mass to the squark mass, we expect 
\K\ to be typically in the range 0.01 to 0.1 (see discussion in Sect. 5.3.1.). 

7.3.1 S-ball Bctryogenesis 

Recall that the perturbations of the AD condensate at some scale A go non-linear 
once t ^ mj^ ~ m^^ (see Sect. 5.10.1), causing the AD condensate to collapse into 
fragments of size A and trapping inside a baryon density. The fragments then relax 
into the state of lowest energy S-balls with a charge of order B. The charge of S-ball 
prevents the soliton from further collapsing and hence the perturbations do not grow 
on length scales smaller than the S-ball radius. Once the length scale going non-linear 
is larger than the final S-ball radius, one expects S-balls to form quite efficiently. 
The time at which a perturbation of scale A goes non-hnear is then given by [49] 



where 



One finds that ~ 34 (44), for d = A{d = 6) directions. In practice S-balls will 
typically turn out to have thick walls with radius R ~ {\K\^/'^m^)^^ . Perturbations 
on this scale which have the largest possible growth in time will go non-linear at 

corresponding to ~ Q.l\K\m^. Assuming that the charge asymmetry corresponds 
to the presently observed baryon asymmetry; r]B ~ 10~^° prior to reheating, then the 
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baryon density is given by 

The charge contained inside a region within a radius it! ^ (|ir|^/^m<^)~^ is given by 
[48, 347, 49] 

Hence with \K\ ^ 0.01, we see that 5-balls of charge larger than 10^^ are hkely to 
form. 

It is quite natural for the MSSM B-balls to decay after the electroweak phase tran- 
sition, see Eq. (344). It is therefore possible to generate the observed B asymmetry 
from the decay of S-balls occurring at relatively low temperatures. This is true even 
if there are rapid L violating interactions or B — L conservation. In a thermal bath 
S-balls could lose their charges by the processes discussed in Sect. 6.5. Therefore we 
naturally have some kind of constraint on the ambient thermal temperature. Combin- 
ing Eqs. (384,375), one obtains an upper bound on the final reheat temperature of the 
Universe [49], which is given by 

T.Mr'i^)(^foe.. (3S5, 

The reheating temperature therefore cannot be larger than 10^ GeV, in order that the 
S-balls have enough charge left over for the purposes of baryogenesis. 

7.3.2 LSP dark matter from i?-ball decay 

S-ball formation from even dimensional d — 4,6 operators is a good candidate for 
generating dark matter in the Universe through S-ball decay into LSPs [48, 49, 361, 
362, 360]. When the 5-balls decay there will be ^ 3 LSPs produced per baryon 
number or equivalently 

n^^3fj,(^!Ihc^nB, (386) 

where fs is the efficiency parameter which denotes the fraction of the total charge 
stored in the 5-ball, which ought to be less than one in a realistic situation, and 
rrin is the nucleon mass. Considering the bounds on baryon number density from 
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nucleosynthesis 0.004 < Qsh"^ < 0.013 [3], and the conservative bound from CMB 
0.004 < Qs/i^ < 0.023, one obtains the upper hmits on the LSP mass [49] 

m, < (17.6 - 20.8)/^ 1 (^) ( A j ' QeV . (387) 

The direct experimental lower limit on the neutralino as a LSP comes from ALEPH 
> 32.2 GeV [401], which requires fs < 0.64. LEP and Tevatron constraints for 
the universal A-teira and gaugino masses lead to > 46 GcV, and for the universal 
scalar masses > 51 GeV, which implies that Jb < 0.45 and Jb < 0.41, respectively. 
This discourages the maximally charged condensate hypothesis for a S-ball formation. 

If one adopts Eq. (386) as it stands, without taking into account of possible LSP 
annihilations, then the amount of dark matter density would be given by [49] 



3 / \m. 



where m„ ~ 1 GeV is the nucleon mass, and we have used the bound Ofi/i^ > 0.004. 
In case of a bino-like LSP, the S-ball formation could be a serious obstacle for the AD 
baryogenesis. 

In this regard a Q-ball with a smaller charge would be beneficial. One option was 
considered in [380], where the authors proposed a gauged B — L symmetry to make 
Q-ball small enough (for a discussion on Q-balls with a local gauge symmetry, see 
6.2.6). The gauged U{1)b-l was assumed to be broken at a scale ~ 10^^ GeV. It was 
argued that in the gravity- mediated models the D-term from U{1)b-l helps forming 
smaller 5-balls from the oscillations of the fiat directions at weak scale. As noticed 
by the authors [380], their mechanism fails to ameliorate the problem for the gauge- 
mediated models. Another solution has been invoked in [402], where it was argued 
that if one takes into account of the Hubble induced radiative corrections to the flat 
direction then for a range of gaugino masses 3H < mi/2 < 5H, the amplitude of 
the AD condensate oscillations is redshiftcd and leads to a formation of considerably 
smaller Q-balls at low scales. The advantage is that the mechanism works for both 
gravity and gauge mediated type AD potentials. 
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Otherwise, because of a large annihilation cross section, Wino and Higgsino hke 
LSP can be rather more promising candidates for CDM if ultimately originating from 
an AD condensate, as pointed out by Fujii and Hamaguchi in [361, 362]. If the LSP 
is not very heavy, say ^ 200 GeV, so that their freeze-out temperature is about 
~ 10 GeV, then S-balls should decay at temperatures below 10 GeV for all 
the produced LSPs to survive. It is quite possible that S-balls are not the only source 
of LSPs and that some of the LSPs produced in B-ha\\ decay will be annihilated 
with the LSPs in the background. Using Eq. (384), one finds that if the B-ha\\ decay 
temperature Td < 10 GeV, there is an upper limit on the reheating temperature [49] 

T., ; 5 X 10« ^ fl^^V (-IS.) GaV . (389) 

fs \ rn^ ) \oikJ VlOGeV; KlO^^V ^ ) 

With \K\ ^ (0.01 - 0.1), and ak ~ 40, Trh < (5 x 10^ - 2 x 10^)/"^ GeV. 
7.3.3 The LSP abundance 

The LSPs produced in 5-ball decays will collide with themselves and with other weakly 
interacting particles in the background and settle locally into a kinetic equilibrium. 
Thermal contact can be maintained until 7/ ~ m^/2Q [285, 18], and a rough freeze-out 
condition for LSPs (if they were initially in thermal equilibrium) will be given by [18] 

nLSp{<J.nnv) ^ Hf-^ , (390) 

^/ 

where (Jann is the LSP annihilation cross-section and the subscript / refers to the 
freeze-out values. The thermally averaged cross section can be written as {(Tannv) — 
a+bT/rriy., where a and b depend on the couplings and the masses of the light fermions 
[18]. For a light neutralino with < mw, neglecting the final state fermion masses 
and assuming an efficient LSP production, so that — 1, one finds for the LSP 
density from Eq. (390), that b Hm^^Tj'^nj^ , and 

^/ = Jj^2i^xTff"e-^-l''f ^ 1.46 X 10-^X ■ (391) 
The LSPs produced in i?-ball decays will spread out by a random walk with a rate 



V determined by the collision frequency divided by a thermal velocity Vth ~ \jT I'^x- 
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Since the decay is spherically symmetric, it is very likely that the LSPs have a Gaussian 
distribution around the central region of radius as pointed out in [48] 

Within this central region annihilation is significant if nLSp(cann'yrei) ~ H. In [48], it 
was shown that the annihilation of LSPs is insignificant provided the Q-ball decay 
temperature is given by 

/ ™ \ 3/16 / 1020 100 \^^^^ 

-'--(loaev) (4) (S ''''' 
The main conclusion is that typically most of the LSPs will survive if the S-ball decay 
temperature is less than a few GeVs. 

The final abundance of the LSPs can be approximately expressed by a simple 
analytical form [361, 362]. Solving the Boltzmann equation analytically for the LSP, 
one finds [362] 



+ \ {av)Mpi{Td - T) 



(394) 



YxiTci) V 45 

where = n^/s. In Fig. (11), we present the numerical solution of Fujii and Ham- 
aguchi for two cases: with and without large entropy production from the i?-balls. In 
the latter case Uq need not be directly related to the present baryon asymmetry. This 
case can be simulated by choosing a very small e and including the radiation energy 
density generated by the S-ball decay in the Boltzmann equation [362] . 

If initial abundance Y^JJI'd) is large enough, the final abundance Y^o for T <^ is 
given by [361, 362] 



/87rV(rd), ^ 
■{av)MpiTd 



(395) 



45 

In this case the final abundance Y^q is determined only by the Q-ball decay temperature 
Ta and the annihilation cross section of the LSP (av) (independently of the initial 
value Y^(Td) as long as Y^{Td) ^ y^^^^P^"^). In terms of the density parameter ^2^, the 
neutralino abundance can be rewritten as [362] 

/0.7\2 f \ /lO-^GeV2\ /100MeV\ / 10 , 
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= 1 GeV (without entropy production) 



e=0.1, N^=10, m^=lTeV, mj,= 100GeV 




time ( GeV"' ) 



= 1 GeV (Q-ball dominant) 



e<e=0.01, Nj,= 100. m^=lTeV, m.^=100GeV 




time ( GeV"' ) 



Figure 11: The evolution of the neutrahno dark matter abundance generated from the Q- 
ball decay for = 1 GeV and {av) = 10~^, 10~^, and lO^^GeV"^, represented by thick solid 
lines. The abundances estimated from the analytic formula in Eq. (395) are shown in dashed 
lines. In the first figure the energy density of the Q-balls has been assumed to be small 
enough with respect to radiation; the parameter values are = 1 TeV, = 100 GeV, 
e = 0.1 and = 10. In the second plot Q-balls have been assumed to dominate the energy 
density before their decay; here the parameters are = 1 TeV, = 100 GeV, e < 0.01 
and A^;^ = 100, from [362]. 

In case of Y^(Td) < y^ppi'o^^ the final abundance is given by 

Y^o^Y^{T,)>e-'(^) , (397) 
\ s /o 

This is the case where the LSP annihilation cross section is small enough, which holds 
for a bino-like LSP, where the relic abundance of the LSPs is directly related to the ob- 
served baryon asymmetry [48] . Unfortunately, a bino-like neutralino will overdose the 
Universe unless we assume an extremely light bino (which is experimentally excluded). 

7.3.4 Which direction? 

Consider the udd and uudde directions, which are both lifted by d = 4 and d = Q 
non-renormalizable operators. For d = Q case the initial value of the field when the 
condensate oscillations begin is given by = 5.8 X lO^^A-i/^ (m^/100 GeVf'^ GeV, 
and for = 4 case = 3.2 x (m^/100 GeVf'^ GeV. Let us compare the 
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lower bound on the reheating temperature for the two cases. For d — 6, it has been 
estimated [49] to be Trh > 0.23 A^^ (m^/100 GeV)^/^ GeV (for d = 6, A ~ 0.003 if 
the strength of the non-renormahzable interactions is set by Mp). Such bound on 
Tj-h fulfills all the requirements for the survival of 5-balls from thermal dissociation, 
diffusion, etc. Repeating the same analysis for the d — A case leads to an upper 
bound Trh t 8 X lO^A GeV (with A ~ 0.1). This is rather hard to satisfy. The 
robust conclusion appears to be that for an efficient B-hal\ production which also 
gives rise to dark matter, should involve the d = 6 udd direction. This particular 
direction is also favored in the sense that fs — 1 can be possible. As long as the 
reheating temperature does not exceed 10^~^ GeV, the baryon to dark matter ratio 
can be naturally accounted by B-hall decays in case of the MSSM fiat directions which 
are lifted by d = 6 operators [48] . 

In the case of the d = A directions B-balls do not form efficiently [49] . For d = 4 
and along the uude direction B — L is actually conserved. In this case only the baryon 
number trapped in the S-balls will survive. Even though it might appear that the 
observed B asymmetry could be obtained if the initial condensate had a large B, this 
would be a rather ambitious program because uude-halh disappear rather quickly in 
a thermal bath as argued in [48]. 

7.3.5 Direct LSP seetrches and i^-balls 

If dark matter is a neutralino, then it should be possible to confirm its existence by 
direct and indirect searches. The direct detection involves the interaction of neutralino 
with matter, which is usually dominated by scalar couplings of relatively heavy nuclei 
A > 20 [403, 18]. The counting rate of the elastic neutrahno-nucleon scattering is 
given by [404] 

R = (^] f l-8xlO"GeVn ( \ ( \ ^^^^^^ 

{m^mNjy Kg -days y VO-SGeV/cm^y \^320 Km/s/ ' ^ ^ 

where p-^ is the density, the average velocity of the neutralino, m„ the mass of the 
target nucleon, and ^ a nuclear form factor. The cross section a is the neutralino- 
nucleon cross-section at zero momentum transfer. These interactions are mediated 
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by heavy Higgs exchanges, or by a sfermion exchange. The forthcoming experi- 
ments, such as CDMS [405], CRESST [406], EDELWEISS II [407], GENIUS [408] 
and ZEPLIN [409], have high hopes in reaching cross-sections up to lO^^'' pb. 

Indirect detection rehes on astrophysics. The annihilation of neutrahno LSPs in 
any astrophysical sources can give rise to fluxes of anti-protons and positrons which are 
usually not seen in the cosmic rays. Other possible ways are neutralino annihilation 
into 27 final states or Z7 final states [410]. The monoenergetic gamma rays with 
energy ~ might not have any competitive background from other astrophysical 
sources, and has a possibility of being detected in the next generation air Cherenkov 
telescopes observing the galactic center, such as VERITAS [411], HESS [412] and 
MAGIC [413]. 

Fujii and Hamaguchi [362], have studied in detail the parameter space where neu- 
tralino production from the decays of 5-balls gives the correct dark matter abundance. 
The allowed region is found where the dominant contribution to the LSP is provided 
by H. A large H content of the LSP enhances the neutralino annihilation cross section 
into W bosons via chargino exchange. There is a small difference in thermally produced 
H and non-thermally produced ones. In the latter case the mass of the neutralino is 
much smaller than thermally produced ones. The annihilation rate of neutralinos 
into 27 final states is enhanced for H-like neutralino (this was previously shown in 
the context of thermally distributed neutrahnos in [410, 414]), and seems to hold true 
even when the neutrahnos are created non-thermally [361, 362]. A typical B-hall decay 
temperature which leads to the desired CDM density is 100 MeV < < (a few) GeV. 
In the anomaly mediation and in the no-scale supersymmetry breaking models, it was 
found that W is the most promising candidate for the LSP in a wide region of the 
parameter space [362]. The conclusion is that if AD baryogenesis is successful within 
minimal SUGRA, then Higgsino and Wino like neutralinos are perhaps the likely dark 
matter candidates [362]. 
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7.4 (5-balls and gauge mediated supersymmetry breaking 

As pointed out by Kusenko and Shaposhnikov, in gauge mediated supersymmetry 
breaking scenarios the salient feature is that the Q-balls are stable against decaying 
into nucleons [151], because the Q-ball energy per unit charge is given by Eq/Q ~ 
f^<t)Q~^^'^ < 1 GeV for ~ 1 TeV. For sufficiently large Q, the Q-ball itself could 
then be a candidate for CDM [331, 151, 356, 351, 415]. 

7.4.1 Baryogenesis and gauge mediation 

A sufficiently large Q-ball in the gauge mediated supersymmetry breaking case will be 
absolutely stable against decaying into nucleons. Although baryogenesis by evapora- 
tion of a Q-ball is not effective, Q-balls nevertheless loose some charge to the ambient 
plasma by thermal dissociation and diffusion, which might create asymmetry in the 
nucleons. An ambitious possibility, but one that would be hard to realize, is to work 
out whether Q-ball evaporation could also lead to dark matter as in the case of gravity 
mediated supersymmetry breaking. 

In the gauge mediated case the baryon asymmetry can be related to the dark 
matter density via 

where £ <S 1 displays the departure from a circular orbit, AQ is the evaporated charge, 
Qq is the density parameter for the Q-balls, pc,o ~ 8/iq x lO^^^GeV^ is the present 
critical density, and n-y^ ~ 3.3 x 10~^^GeV^ is the present photon number density, 
with hp, ~ 0.7. The evaporated charge should yield the baryons while the remaining 
charge is in Q-ball dark matter. The baryon-to-dark matter ratio is given by [351] 

rs^^- VB^Q-'" ~ 10".-',teSJ- (^) Q-^'\ (400) 
Q epc,Q\lQ ^ VTeV/ 

The total evaporated charge from Q-balls in a thermal bath is given by Eq. (374). 
Requiring that Q > AQ, the charge of the Q-ball should be sufficiently large [351]: 

g>1.2xl0V/V/^<(;g|)". (401) 
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The value of e is not an independent quantity and can be related to the charge of a 
Q-ball. At the beginning of the condensate rotation H ~ m'^{T)/ipo ~ T'^/(po, where 
T ~ (MT^f^Hy^^ is the instantaneous temperature prior to reheating. Combining 
these two pieces of information, one obtains T ~ T^^^ M / ipQ. The net charge of the 
Q-ball is then related to the baryon number density through ub ~ tbSu^ ~ tbSoj^Pq, 
where u ~ m'^{T)/ipo ~ T'^/<Po is the rotation frequency of the condensate. With the 
help of Eq. (400), one finds [351] 

UB ~ IO'^Vb^q' (;^) TVog-^/^ (402) 

Kasuya and Kawasaki also estimated the net charge of the Q-ball, given by the initial 
amplitude of the condensate and the reheating temperature, as [351] 
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From Eq. (356) one obtains Q ^ (3 {(po/T) for £ ~ 1 (for £ -C 1, one replaces (3 by 
/3' = 7/3 with 7 ~ 0.1). 

Prom all these considerations one can estimate the amplitude of the condensate 
and the charge of the Q-ball. The amplitude is given by [351] 

(R \ ^/^'^ / \ -2/5 

Inserting this in the expression T ~ Trh\jM/ipQ, one obtains [351] 

T ~ 2.3 X 10' 8-'/''^^'/' ( X ( -Ar^] { GeV , (405) 

^ VlO^GeVy V6xlO~V vTeV; ' ^ ' 

so that the charge of the Q-ball reads [351] 

Now, with the help of Eqs. (401) and (406), one obtains [351] 

llO-ioAlOSGevJ ""UxlO-s UevJ " ^^"^^ 



£ ~ 1.5 X 10-^n 

The above equations determine the parameter space for Q-ball baryogenesis to coexist 
with Q-ball dark matter in gauge mediated models. 
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7.4.2 Generic gauge mediated models 

Now the task is to consider the parameter space for a reaUstic MSSM flat direction. 
A detailed analysis can be found in [351], where it was pointed out that the d — 5,6 
directions are favorable. For d — A, the charge does not accumulate enough to survive 
as a dark matter relic while d = 7 requires unnaturally small values of £. It is worth 
pointing out that as long as the ambient temperature of the plasma is sufficiently 
high with T > (m3/2<^o)^^^ [loglv'o/^^)] logarithmic term dominates over the 

gravity- mediation term. For d — 5,6, the requirement on temperature is thus T > 
10^-6 GeV. 

On the other hand, if the Q-ball forming scalar field has a large vev, or if the 
reheating temperature is extremely low, then gravity mediation effects should also be 
taken into account. However, as kasuya and Kawasaki pointed out in [351], in order 
to produce enough baryon asymmetry and simultaneously produce surviving Q-balls 
for dark matter, one requires a considerably large gravitino mass 777.3/2 > 10^ GeV, 
which is an unacceptable value for the gravitino mass within the gauge mediated 
supersymmetry breaking [39]. 

One could also consider a generic model for gauge mediation where the scale of the 
logarithmic potential is larger than m^, such that [349, 351] 

, _-_plog I — „ , „ -.- ,„ 
Vr^{ ^ ^\Ml) (408) 



ml^'' < Ms), 



where Mg is the messenger mass scale. In this particular case the condensate will start 
oscillating at large field amphtudes. The Q-ball will form at a large vev. The mass 
and the size of the Q-ball are now given by [349, 351] 

Mq^MfQ^'^, R^Mp^Q^'\ u (409) 

If Mp > T, then the Q-ball could be stable against decaying into the nucleons, provided 
the Q-ball mass per unit charge is smaller than 1 GeV. This condition holds for 



159 



The ambient temperature is smaller than Mp, but could be larger than m^. In this 
case the Q-ball would rather evaporate. The total evaporated charge can be estimated 
to be [351] 

Ag ^ 10- (^^y^' (^^y^^ X gvi. . (411) 

^ V102 GeVy Vl06 GeVy ^ ^ ^ 

The survival condition: Q > AQ then imposes a bound on an initially accumulated 

charge to be 

Q > 10" ( '"^ ^ ( \ (412) 

VlO^GevJ llO^GevJ ' ^ > 

While combining with Eq. (399), one obtains a relationship between the baryon number 
and the amount of dark matter as [351] 

o-o.'"r^-(^)-''■X^)-(^)-^ («3) 

The consistency condition requires the initial amplitude of the AD condensate to be 
(/?o ~ 10^^ GeV. Note that when all the above conditions are taken into account, then 
there is hardly any region in the parameter space which allows for the required baryon 
asymmetry from evaporation of Q-balls together with enough surviving Q-balls to 
provide the dark matter, as was concluded by Kasuya and Kawasaki in [351]. 

7.4.3 Late formation of gauged Q-balls 

In [351], the authors have also considered the late formation of Q-ball at a scale when 
gravity mediation and thermal logarithmic correction terms dominate the potential at 
large and small scales, respectively. The Q-ball forms when the instabihty band enters 
the horizon with a wavelength and the angular velocity ujeq ~ ^Iql ^eq)~^ ~ ''"3/2 ' 
The number density of the Q-balls is found to be n^q ~ T'^Jm^^j^^ while the charges arc 
Q ^ (<^/^e?)^ [351]. The exact relationship depends on the helicity of the condensate. 
In this late formation scenario the only realistic flat directions turn out to be d = 6 
with Q ~ 10^4 - 10^1, Trh ~ 1.0 X 10^ - 30 GeV, and Mp ~ 10^ - 10^ GeV; and d = 7 
with Q ~ 10^5 - 10^2, ~ 60 GeV - 10 MeV, and Mp ~ 10^ - 10^ GeV [351]. 

Kasuya and Kawasaki also repeated their analysis when temperature effects arc 
negligible. In this case the Q-ball forms when V ~ Mplog (^(pl^/Mg^ ~ f^l/2V^eq-> 
where (peq ~ Mp/m^i2- The AD condensate fragment just after its amplitude becomes 
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smaller than ipeq. Consistent scenarios arise only in d = 6, 7, 8 cases. For the allowed 
parameter space, one should require 777,3/2 ~ 0.1 GeV, Mp ~ 10^ GeV, T^h ~ 5 GeV, 
and Q ~ lO^o. 

These many proposals indicate very clearly that in order to find a consistent cos- 
mological scenario where Q-ball evaporation leads to baryon asymmetry while the 
survived Q-balls act as a dark matter candidate requires a stringent condition on 
Mp < 10^ GeV. Note that Mp ~ 10^ GeV is required to provide the right spectrum 
for sparticle masses, see [297, 39]. 

7.5 (5-balls as self-interacting dark matter 

Recently, Q-ball has been proposed as a candidate for self-interacting dark matter by 
Kusenko and Steinhardt [357]. Such a consideration is motivated by the fact that col- 
lisionless CDM appears to have certain discrepancies between numerical simulations 
and observations. The halo density profiles, and the number density of satellite galax- 
ies, do not match well with the observations [416, 417]. A possible remedy is dark 
matter that has fairly strong self-interactions [417], a situation which however is not 
easily achieved in the standard particle physics models. 

The self-interaction cross section and the mass of the dark matter particles under- 
going elastic scattering should satisfy the relation s — aon/'mDM ~2xl0^ — 3x 
10"^ GeV~^. This might change if one considers other processes such as dark matter 
annihilation. 

Q-balls, being extended objects, can certainly have a large geometric cross-section, 
but it was found that in order to match the required cross-section to mass ratio s, the 
AD particles should have a very low mass scale ~ 0(1) MeV. This requires that the 
charge of the Q-ball should also be very small, i.e. Q < 10^. Such Q-balls would not 
originate within MSSM but could be possible in some extended theories. There are 
however a number of issues concerning production of such a small charged Q-balls, 
their thermal distribution, their evaporation, annihilation and scattering should be 
taken into account consistently. It has been argued that thick wall Q-balls do not 
seem to have much admissible parameter space [358], while thin wall Q-balls have a 



161 



slightly better chance to succeed with a vev cpo > C(MeV). In any case, such a small 
AD condensate amplitude seems hard to obtain, although not completely impossible, 
in order to reconcile with the existing particle physics models. 

Recently, it has been pointed out [418] that future experiments should be able to 
discern the spatial extent of the dark- matter particle. It was noticed that the extended 
objects such as Q-balls leave its distinct imprint on the spectrum which falls off very 
fast with increasing energy. The signal is primarily dominated by the low-energy events 
near the threshold. 

7.6 Direct searches for gauge mediated (^-balls 

If the initial charge of a Q-ball is larger than the evaporated charge, the Q-hall sur- 
vives and contributes to the energy density of the Universe. Let us take as an example 
S-balls, such as U2did2, which in gauge mediated case are unstable because of the pres- 
ence of baryon number violating operators (required in the first place to charge up the 
condensate) with dimensions larger than 5 [415]. Nevertheless for Q > 10^'^, the life- 
time of the Q-ball is in fact greater than the age of the Universe {to ~ 10^° years) [415]. 
Such Q-balls are potential candidates for CDM. 

Prom Eq.(374), we can read the limit on the initial charge of a Q-ball 

QinH>7.4xlO^^(;^)"''^'\ (414) 

On the other hand, in order for the Q-ball to be stable against decaying into nucleons, 
i.e. Eq/Q<1 GeV, one finds [351] 

Q>lO"{^y. (415) 

The condition relating baryon number and dark matter is given by Eq.(401), and can 
be written as 

Q<l0^^e'/^[^y\ (416) 

where we have taken tjb ~ 10"^'^ and Qq < 1. 

If we do not impose the condition that the evaporated charge accounts for the 
baryons in the Universe, the only constraint is that the energy density of Q-balls must 
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not exceed the critical density. As mentioned in the previous Section, this condition 
is given by Eq.(406). 

Let us now assume that the gauge mediated Q-balls indeed make up the dark 
matter in the galactic halo. Then the corresponding number density is given by [355] 



and the Q-ball flux will be determined by ^ ~ (l/47r)nQV ~ lO'^Q^ ' (1 TeV/m^) 
cm^^ s~^ sr"^, where we have taken v ~ 10~^c. The number of events will obviously 
depend on the surface area of the detector. In case of Super-Kamiokande the number 
of events could be estimated to be ~ {10'^'^ /QsY^^i^ TeV/m<^) yr~-^ [355], where 
the surface area of the water tank is ~ 7.5 x 10^ cm^ [419]. 

The Q-balls can be detected via Kusenko-Kuzmin-Shaposhnikov-Tinyakov (KKST) 
process [355]. Note that in the absence of any fundamental singlet a S-ball would 
necessarily break SU (3) while the electroweak symmetry could be restored inside, 
provided the Higgses do not obtain vevs or if the flat direction is already a SU{2) 
singlet such as udd. The instanton mediated baryon number violation will be much 
weaker in these cases because of the limitation on the size of the instanton, which 
cannot be larger than the size of the Q-ball. When nucleons coUide with a Q-ball, they 
enter the surface layer of the Q ball, and dissociate into quarks, which are converted 
into squarks via gluino exchange. In this process, Q-balls release ~ 1 GeV energy per 
collision by emitting soft pions. 

For an electrically neutral Q-ball the absorption cross-section is quite large and is 
determined by the soliton size to be [355] (Tb ~ cm . The esti- 

mated mean free path in matter is given by A ~ 1O^'^A(1O^^/Qb)^''^(^0/1 TeV)(l gjcvn? j p) 
cm [355, 51], where A is the weight of the atomic nucleus with a density distribution 
p. Obviously absorption of quarks takes place at a higher rate than the collisions of 
Q-balls with nuclei. The BAIKAL experiment [420] sets a limit on the monopole flux 
which also gives a lower bound on the charge of a Q-ball as > 10^^ for ~ 1 TeV. 

If the Q-balls are electrically charged, for instance when the selectrons obtain a 
large vev along the QQQLLLe flat direction, the detection prospects worsen, because 
of two reasons. First, there will be a Coulomb barrier which will prevent absorption 




(417) 
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Figure 12: The bounds on the charge of an electrically neutral Q-ball with respect to the 
supersymmetry breaking scale. Various shaded regions of the parameter space are already 
ruled out. The line "B-Ball Stability Limit" marks the region below which i?-balls are not 
stable. The allowed region Q > 10^^ can be found only in the upper part of the figure, from 
[51]. 

of the incoming nuclei. Secondly, the absorption cross section will be determined by 
the Bohr radius cr ~ vrr^ ~ 10^^^ cm~^, and therefore the corresponding mean free 
path length will decrease to A ~ 10~^A(1 g/cm^/p) cm [355]. The present limit on 
electrically charged Q-ball comes from the MACRO search [421] with a flux JF < 
1.1 X 10~^^ cm^^ s^^ sr^^, which constraints the charge of a S-ball to > 10^^. 

A wide range of charges are already ruled out, see Fig. (12) from [51]. The only 
allowed range of charges is Q ~ 10^^ — 10^^ for ~ 1 TeV — 100 GeV. The range 
however depends on £, which one tacitly assumes to be of order one in Eq. (416). 
For smaller values of e the range of allowed charges will further reduce. In a generic 
logarithmic potential a charge Q ~ 10^^ and Mp ~ 10^ GeV is still allowed [351]. 
In many cases such as in delayed Q-ball formation [351], the charge of a Q-ball is 
large ~ 10^^. Such Q-balls could be detectable in the Telescope Array Project or the 
OWL-AIRWATCH detector. 

There are also astrophysical ways for detecting stable Q-balls. Charged Q-balls 
can dissipate a large amount of energy quite efficiently before they can be detected. 
The Q-balls might leave a track behind whose stopping range is roughly 1000 m 
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for a charge Qb < lQ^^{m^/l TeV)"^/^ [355]. Electrically neutral Q-balls would 
hardly have any impact when they pass through a planet like ours. The loss in Q- 
ball kinetic energy would be almost undetectable; the velocity would decrease by 
5v/v ~ 10-2Q-^/^(l TeV/m^f. 

Q-balls could also be captured in sufficiently dense stars such as inside the core 
of neutron stars. This will result in an increase in the temperature of the neutron 
star studied in [355], but an insignificant enhancement ~ 0.01(Qb/10^'^)^^^^^ may 
not result in any observable consequences. There is still a possibility that the gradual 
accumulation of Q-balls may seal the fate of a neutron star by decreasing its mass and 
thereby reaching the critical condition which may lead to a supernovae explosion [355]. 
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8 Flat directions other than MSSM 



Flat directions, fragmentation of scalar condensates and Q-balls are generic features 
that could be encountered in many cosmological models with scalar fields. Examples 
include such cornerstones of modern cosmology as inflation, as well as the currently 
popular particle physics models with extra dimensions. It has also been suggested that 
the motion of a complex flat direction condensate field, spinning in a C/(l)-symmetric 
potential, could give a dynamical explanation for the dark energy [422] , although the 
situation could be complicated by copious Q-ball formation [423]. 

8.1 Fragmentation of the inflaton condensate 

By definition, the inflaton is a homogeneous scalar condensate with a small quantum 
induced spatial fluctuations. Reheating of the Universe, which is a consequence of the 
inflaton condensate break-up and decay, may take place via ordinary smooth pertur- 
bative decay of the condensate [229, 22] or via a non-perturbative process dubbed as 
preheating [236, 237], which typically involves an amplification in some of the fluc- 
tuation modes as well as the fragmentation of the inflaton condensate. Reheating 
dynamics depends very much on the form assumed for the inflaton potential, and for 
some choices, the inflaton condensate may also form Q-balls, see [353, 354] . 

8.1.1 Reheating as a surface effect 

Usually the process of reheating is taken to be entirely a volume effect. This can be 
problematic, however, especially if the scale of inflation is high, i.e. Hinf ~ 10^^ — 
10^^ GeV, as it is sometimes assumed in order to provide the right magnitude for the 
density perturbations, and also for reasons that have to do with non-thermal heavy 
dark matter production or exciting right-handed Majorana neutrinos for leptogenesis, 
etc. [136]. As it is well known, the entropy thus dumped into the Universe may pose a 
problem for big bang nucleosynthesis by overproducing gravitinos from a thermal bath; 
an often quoted bound on the reheat temperature is Trh < 10^ GeV [258]. Obtaining 
such a low reheat temperature is a challenge for high scale inflation models. (One 
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way to solve the gravitino problem is to dilute them via a brief period of late thermal 
inflation [172]). 

A novel way to avoid the gravitino and other moduli problems is reheating via the 
surface evaporation of an inflatonic soliton. Compared with the volume driven inflaton 
decay, the surface evaporation naturally suppresses the decay rate by a factor 



where L is the effective size of an object whose surface is evaporating. The larger 
the size, the smaller is the evaporation rate, and therefore the smaller is the reheat 
temperature. 

Reheating as a surface phenomenon has been considered [353, 354] in a class of 
chaotic inflation models where the inflaton field is not real but complex. As the 
inflaton should have coupling to other fields, the inflaton mass should in general receive 
radiative corrections [9] , resulting in a running inflaton mass and in the simplest case 
in the inflaton potential that can be written as 



where the coefficient K could be negative or positive, and m is the bare mass of the 
inflaton. The logarithmic correction to the mass of the inflaton is something one 
would expect because of the possible Yukawa and/or gauge couplings to other fields. 
Though it is not pertinent, we note that the potential Eq. (419) can be generated 
in a supersymmetric theory if the inflaton has a gauge coupling [49, 348, 330] where 
K ~ —{a/87i){my2/''^)y where mi/2 is the gaugino mass and m;^ denotes the slepton 
mass and a is a gauge coupling constant. It is also possible to obtain the potential 
Eq. (419) in a non-supersymmetric (or in a broken supersymmetry) theory, provided 
the fermions live in a larger representation than the bosons. In this latter situation the 
value of K is determined by the Yukawa couphng h with K — — C(/i^/167r^), where C 
is some number. 

As long as \K\ ^ 1, during inflation the dominant contribution to the potential 
comes from m^|$p term, and inflationary slow roll conditions are satisfled as in the 
case of the standard chaotic model. COBE normahzation then implies m ~ 10^^ GeV. 




(419) 
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If K < 0, the inflaton condensate feels a negative pressure (see, Sect. 5.9.1, for the 
discussion on negative pressure of the AD condensate) and it is bound to fragment into 
lumps of inflatonic matter. Moreover, since the inflation potential Eq. (419) respects a 
global U{1) symmetry and since for a negative K it is shallower than m^|$p, it admits 
a Q-ball solution (see Sect. 5.9.3). Comparing with Q-balls along the MSSM flat 
directions, here the major difference is that the inflatonic condensate has no classical 
motion along the imaginary direction as usually required for a Q-ball solution. 

8.1.2 Q-balls from the inflaton condensate 

As pointed out in [353] , there are quantum fluctuations along both the real and imagi- 
nary directions which may act as the initial seed that triggers on the condensate motion 
in a whole complex plane. The fluctuations in the real direction grow and drag the 
imaginary direction along via mode-mode interactions, as illustrated by 2 dimensional 
lattice simulation in [354], see Figs. (13). The flrst plot shows the hnear fluctuations 
without rescattering effects; scattering effects are accounted for in the second plot. 
The late time formation of inflatonic solitons is shown in Fig. (14). Q-balls were ob- 
served to form with both positive and negative charges, as can be seen in the first 
plot of Fig. (14), while keeping the net global charge conserved. Inflatonic Q-balls are 
of same size because the running mass potential resembles the MSSM flat direction 
potential in the gravity mediated case, where the Q-ball radius is independent of the 
charge. 

Q-balls of size R ~ \K\^^^'^m~^ form when the fluctuations grow nonlinear (see 
Eq. (324) [353, 354]). Since the growth rate of fluctuation is ~ |i^|m, the Hubble 
parameter at the formation time can be estimated as Hf ~ 7|ir|m, where 7 is a 
numerical coefficient less than one. For <^ 1, we can approximate the decrease in 
the amplitude of the oscillations by 0/ ~ (j)i{Hf/Hi) as in the matter dominated era, 
where (pi ~ Mp denotes the amplitude at the end of infiation in the chaotic model, 
and Hi m when the oscillations begin. The total charge of a Q-ball is given by 
Q ~ (A7r/3)R^nq ~ (l/9)(3C^-f^\K\'^R^mM^, where Uq = 0o ~ C0/, and /3 < 1 

and C ^ 1 are numerical factors. 
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Figure 13: The first plot from left shows the instability bands of the homogeneous mode 
of the inflaton along the real (solid) and imaginary (dotted) directions. The second plot 
shows the result of lattice simulation in the real and imaginary directions, together with the 
evolution of the homogeneous mode. The third plot shows the power spectra of fluctuations 
at late times. All plots assume K = —0.02, from [354]. 

Given an inflaton coupling to fermions of the type h4>ipip it has been shown that 
[353, 354] reheating is driven by surface evaporation of inflatonic Q-baUs for relatively 
large Yukawa couplings h < 1. In general K and h are not independent quantities but 
are related to each other by \K\ ~ C(/i^/167r^). If the inflaton sector does not belong 
to the hidden sector, it is very natural that the inflaton coupling to other matter fields 
is relatively large, i.e. h > (m/Mp). In this regime the evaporation rate is saturated 
by Eq. (340) and [354] 

Note that the decay rate is determined by the ratio m/M-p ~ 10~^, which is fixed by the 
anisotropics seen in the cosmic microwave background radiation. Even though we are 
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Figure 14: The first plot shows the charge density distribution in a small sub-lattice at late 
times. The second plot shows inflatonic solitons forming in 3D lattice. Here K = —0.02, 
from [354]. 

in a relatively large coupling limit, the decay rate mimics that of a Planck suppressed 
interaction of the inflatonic Q-ball with matter fields. Fermionic preheating [246] is 
not a problem in this case because the whole infiaton energy is not transferred in this 
process and the energy density stored in the fermions remains small compared to the 
infiaton energy density, as argued in [354] . Fermions cannot scatter infiaton quanta off 
the condensate [246] , unhke in the case of bosonic preheating [237] . Inflatonic Q-ball 
formation is reminiscent of bosonic preheating [424] due to the presence of attractive 
self coupling of the infiaton, which stems from the logarithmic term in the potential. 

8.2 Affleck-Dine baryogenesis without MSSM flat directions 
8.2.1 Leptogenesis with sneutrino 

It is interesting to note that it is possible to mimic the AD baryogenesis even without 
MSSM fiat directions. This can happen in F-term hybrid inflation with a superpoten- 
tial [137] 

W = -A^^ + XS^^ + , (421) 

where A fa 6.5 • lO^^e^/^ GeV, or equivalently A ^ 1.3- lO^^jr;] A-^/^ GeV, where e, < 1 
are the slow roll parameters. The superpotential Eq. (421) extends the usual hybrid 
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case Eq. (135) in order to provide a natural initial condition for which should be 
zero within the accuracy 10~^ [224], while the inflaton should have an initial value 
close to the Planck scale. ^ is identified with the U{1)b-l carrying right-handed 
neutrino which alleviates the problem of initial conditions through the dynamics of 
Eq. (421) [224]. During inflation the system is trapped in a false vacuum with (f) — 0. 
Once (f) decays, sneutrino ^ obtains a vev of order of the Hubble parameter at the end 
of inflation. It then starts oscillating around the origin and due to dynamical breaking 
of U{1)b-l, a net B — L asymmetry is generated [137]. Right after infiation the real 
and imaginary component of ^, -01 and -02 have a relative sign difference which induces 
a helical motion on ^ similar to the MSSM flat direction condensate. In this model 
the neutrino mass is generated by the ordinary see-saw mechanism [32] 

= ^ ^ X 1.2 • 10-^ eV , (422) 

(where we have taken {1^2) ^ 170 GeV). The Hubble parameter during inflation is 

Consideration of the detailed dynamics and decay of the condensate leads to an esti- 
mation of the lepton asymmetry, which is given by [137] 



where the numerical factor is X = m.m{l, x}Ck~^, with C ~ 1 coefficient if the field 
^ gets a supergravity induced mass term during inflation, otherwise C — S/N^. The 
reheat temperature appears to respect the gravitino bound. 

8.2.2 AD baryogenesis in theories with low scale quantum gravity 

In low scale gravity models the weakness of gravity arises because of extra spatial 
dimensions, the scale of which could be as large as mm, [425]. In the simplest models, 
the SM fields live on a three dimensional brane while gravity can also permeate the 
bulk [425, 426]. The fundamental scale is ~ 0(TeV) in 4 -|- 0? dimensions, which is 
related to the the volume suppression and to the effective four dimensional Planck 
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mass Mp by a simple relationship [425] 

= M^+'^Vd . (425) 

This automatically sets a common size for all the extra dimensions at 60 • For two 
extra dimensions with M^, = 1 TeV, the size is 0.2 mm, just below the current ex- 
perimental limits from the search for deviations in Newtonian gravity [427]. Recent 
astrophysical bounds based on neutron stars suggests M* > 500 TeV for two extra 
dimensions [428] . Naturally, low scale gravity models have important implications for 
collider experiments [429] and for cosmology [430, 431] (for a review on large extra 
dimensions, see [432]). 

Although quite attractive from a particle physics point of view, large extra dimen- 
sions bring along a host of cosmological problems. There are dynamical questions 
regarding the stabilization of the size of the extra dimension(s), or equivalently the 
vev of the radion field, whose mass can be as small as 0{eV) (for two large extra 
dimensions) [425]. Cosmologically stabilization should take place very early by some 
trapping mechanism as discussed in [174]. Another challenge is how to realize inflation 
in these models. There have been many proposals [431], such as invoking a SM singlet 
scalar hving in the bulk [433]. There is the problem of the Kaluza Klein (KK) states 
of the graviton and any other fields residing in the bulk, which above a certain tem- 
perature known as the normalcy temperature, should fill the Universe. The normalcy 
temperature is constrained by cosmological considerations to lie in the range 1 MeV to 
100 MeV [425, 430, 434, 435, 436], and the reheat temperature should be lower than 
the normalcy temperature. 

In addition, one must not only forbid dangerous higher order operators which can 
mediate proton decay, but also ensure that such operators are not being reintroduced 
by whatever mechanism is responsible for baryogenesis. This tends to make baryogen- 
esis in low scale gravity models quite difficult. Moreover, it has been argued [434, 435] 
that leptogenesis is not a viable option; obviously at very low reheat temperatures 
required by the normalcy temperature, the sphalerons cannot be activated (a possible 
way out of this problem could be to increase the reheat temperature by increasing the 
number of large extra dimensions to six [437]). 
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One solution to this predicament was provided in [435] (see also [438]). There 
U{1)^ carried by a gauge singlet x was broken dynamically in order to provide a small 
asymmetry in the current density. Baryon asymmetry is produced by the decays of 
X and X into SM quarks and leptons, analogously to the old AD baryogenesis. The 
decay channels are constrained because quarks and leptons must carry a non zero 
global X charge. This prevents x~X asymmetry to be transferred into non-baryon- 
number violating interactions such as interactions involving the Higgses. 

A consistent model for baryogenesis can be constructed [435] along the lines of 
sneutrino leptogenesis discussed in Sect. 8.2.1. It requires a hybrid inflaton sector, 
described by the fields (f) and N, and the flat direction field x, which are promoted to 
the bulk. This ensures the right amplitude for the density perturbations and provides 
enough baryon asymmetry towards the end of reheating [434, 435]. The potential for 
the zero modes in 4 dimensions can be written as [435] 

Vau{<P, N, xu X2) = n\ ' N\xl + xl) + ^ (^) ' ixl + xlT 

+4[^)\n{xI-xI), (426) 

where Ki, 1^2,1^3 are order one constants, and Xi and X2 are the real and imaginary 
components of the complex field x- Note that all the terms are Planck mass suppressed 
because Eq. (426) is an effective 4d potential derived from a higher dimensional La- 
grangian by integrating out the extra spatial dimensions. Since during inflation the 
auxiliary field = 0, the flat direction condensates Xi and X2 are massless and the 
potential is almost flat. 

The final x asymmetry was found to be given by [435] 

s ~ 27P I, Mp ) [hJ - 27 [J \\nJ UJ ' ^^^^^ 

which at the same time provides an upper bound on the baryon to entropy ratio. As 
an example, taking ~ 100 MeV, ~ 100 TeV, the number of e-foldings A^e ~ 100, 
and couplings of order one, yields an asymmetry of order ~ lO""^*^. 
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8.3 Solitosynthesis 

Solitosynthesis is a mechanism of charge accretion with the help of pre-existing small 
Q-balls in a charge asymmetric background [439, 440, 441, 443, 444]. The accumulation 
of charge and forming a large charged Q-ball has been shown to be quite efficient 
especially in a finite temperature thermal bath. 

8.3.1 Accretion of charge by Q-balls 

Consider thermodynamics of Q-balls surrounded by massive non-relativistic fermions 
ip. The number densities of Q-balls and ip particles are governed by the Boltzmann dis- 
tributions nQ(T) = gQ{MQT /2TJ:f/'^ exp[(/iQ-MQ)/T], andn,^ = 5(^(m^T/27r)^/2 exp[(/XQ- 
'm^)/T], where gg is the partition function for the Q-ball and — 2. The respective 
chemical potentials are denoted hy /i. In chemical equilibrium the absorption and 
evaporation of charge is equally possible, i.e. (Q) + t/j -(r^ {Q + 1) , which relates the 
two chemical potentials through fiQ = Qfi^. The number density of Q-balls can be 
expressed in terms of the number density of the fermions [440, 444] 

(^) ' " ' 

where Bq = Qm^ — Mq is the binding energy per charge of a Q-ball. When Bq grows 
with Q, the formation of large charged Q-ball is likely. The interactions between Q-ball 
and ip quanta leads to a chemical equilibrium when [440, 444] 

n^V^UabsiQ) = riQ+iTevapiQ + 1) ■ (429) 

The accretion and evaporation rate from a charge Q-ball is given by Saha equation: 

^a6s(Q) ^evapiff) ; 



<^abs{Q) CTahsiQ - 1) 



(430) 



where — (T/27rmy,)^/^ is the mean velocity of ip particles. The charge of a Q-ball 
grows when TahsiQ) > T^evapiQ)- The necessary charge asymmetry in ip quanta affects 
the Q-ball abundance. Charge conservation requires 

N ^ - n^* + ^ Qnq ^ Q*nQ* = -q^n^ , (431) 
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where r]^ — n^— riii;* jn^ is the charge asymmetry. If 77 is close to zero, annihilations of 
i) quanta will be effective and there will be virtually no solitosynthesis. The number 
density of the stable quanta is = ?7^n-y, and the charge asymmetry is given by 
[440] 

= 2.5 X lO-^Q^/i^^ (432) 

Obviously the over-closure limit Q^/i^ < 1 yields an upper bound on 77^ in terms 
of m^. At high temperatures large Q-balls arc suppressed by a small asymmetry 
factor ~ At lower temperatures the abundance is dominated by Q^^"^, where 

Mq oc m^Q [440]. 

The Q-ball starts growing at a temperature when VabsiQ) > fevap{Q)- This can 
be estimated from the Saha equation [443, 444] 

J. ^ + Mq., - Mq 

where c ~ 0{1) number for large enough Q. 

When the absorption of Q-charge from the surroundings freezes out, sohtosynthesis 
stops. There are two distinct era which one may consider; solitosynthesis freezes out 
during the radiation dominated epoch, or during the matter dominated epoch. The 
demarcation temperature is Teq ~ 5.5(r2o^^)^^C~^ eV. The freeze-out temperature is 
defined by T[{Q) + -0 ^ (Q + 1)] < H{T). For a geometric cross section Uabs ~ t^R^i 
the freeze-out temperature is given by [444] 



10-9 






IceV 


10-13 






ICeV 




(434) 



where ^ = [g^,{Ti:,)—g^,{T)Y^^ ^ the temperature signifies the decoupling temperature 
of a Q-ball with a rest of the plasma, and /3q = {?>m^^/ Ai^ujip^^/^ . One should require 
that Tp <Ts, which is satisfied for a very low ■0-inass. In order to have some feeling 
for the numbers involved, we note that for the smallest charged Q-ball with Q = 2, 
the freeze-out temperature is smaller than T, only when < 0{1) GeV. During 
the matter dominated epoch, due to the gravitational clustering, Q-ball synthesis is 
favorable but for ambient temperature < 0{1) eV, < 0{1) MeV. 
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8.3.2 Phase transition via solitogenesis 

An interesting observation has been made in [442, 443], where it was pointed out that 
sohtosynthesis may also lead to a first order phase transition. The gradual absorption 
of charge may lead to a critical charge Q — Qc for which the false vacuum inside the 
Q-ball becomes unstable and expands as the the whole space is filled with the true 
vacuum. The value of the critical charge is determined by minimizing the energy of 
the Q-ball with respect to the radius when dE/dR — and d^E/dB? — are satisfied 
simultaneously. In a thin wall limit when 1, the critical charge is given by [443] 

^ 81 U{^,fn ' (435) 



where = dip^2{U{ip) - {uj^/2)ip'^). 

One can also enquire how large the critical charge should be in order to facilitate 
the destabihzation of the otherwise cosmologically stable false vacuum. It has been 
found that the decay of a metastablc false vacuum at zero temperature requires a 
small charge Qc ~ 28{ipo/U{(poV/^) for A - (100 GeV)^ [443]. Similar considerations 
at finite temperature result in a larger critical charge Q > Qc ^ 14:6ipoT/U{(po,Ty^^ 
for the phase transition to proceed. 

Solitosynthesis-catalyzed phase transitions in supersymmetric models would re- 
quire a local violation of lepton or baryon number. As an example (see [443]), one 
could consider a lepton number violating local minimum along H^LlLji ^ 0. This 
false vacuum can decay into the standard true vacuum but the sphaleron induced 
transition rate is almost negligible compared to the cosmological time scale because 
SU{2) X U{1) is broken. If the typical mass scale of squarks and sleptons is consid- 
erably heavier than 1 TcV, L-balls will accumulate when > Tp pa and can 
catalyze a phase transition within one Hubble time at temperatures Ts [442, 443]. 
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